- 在线时间
- 0 小时
- 最后登录
- 2018-3-16
- 注册时间
- 2018-3-16
- 听众数
- 1
- 收听数
- 0
- 能力
- 0 分
- 体力
- 3 点
- 威望
- 0 点
- 阅读权限
- 10
- 积分
- 1
- 相册
- 0
- 日志
- 0
- 记录
- 0
- 帖子
- 1
- 主题
- 1
- 精华
- 0
- 分享
- 0
- 好友
- 0
升级   20% 该用户从未签到
 |
本帖最后由 419595891 于 2018-3-16 13:51 编辑
8 A* W0 m7 D# X L" f# k1 ?; i. [ |; g' [! N
这是一篇文献里给出的script,但我放到S-plus里跑不出来,有错误,有可能是我不会用S-plus,所以请教大神下面的script是否有错,没有错的话怎么解决错误?能不能改动换到R里运行?
; S: O1 `6 B; H: P* F7 L$ I/ c$ m
7 r; a' m: a" J' u/ S8 c: vlibrary(Design) 1 \2 w+ H- \. Q1 W0 T/ ]- ^9 {* u4 w) ]5 U
reri <- function(datsam)
5 `" ~4 [# F: @. j& \{ ; D& M) Q/ ]1 G& {5 N! Q
fitlr <- glm(outcome ~ determinant A * determinant B, family=binomial, data=datsam)
1 {, E3 |: b/ Z# S* u6 L0 jreri <- exp(fitlr$coef[2]+fitlr$coef[3]+fitlr$coef[4]) - exp(fitlr$coef[2]) - exp(fitlr$coef[3]) + 1 * t1 k5 v3 j8 E2 p8 m y) Z
}
5 D# A, r$ Z+ H# x6 }summary.bootstrap(bootstrap(dataset, reri(dataset), B=10000), probs=c(0.025,0.5, 0.975)) ( O3 R. U V5 ^6 U0 P
4 P0 Y, h4 g; ], M
————————————我是分割线—————————————————————————————— ! w4 f$ [7 P" a/ f/ \; }& b. b# y
; ?3 s! r, J4 N8 _
# {) T; R7 @' ?7 R
- D! I$ h) x! {) R library(Design)* z8 o9 x+ q' h* O; J$ J
> reri <- function(datsam)
3 p9 V3 R1 {3 |) y9 w9 T$ B+ m+ {fitlr <- glm(DM.0813 ~ smoking * rs163184, famliy=binomial, data=datsam)/ N% k: R, d& g; A& d2 Z
+ reri <- exp(fitlr$coef[2]+fitlr$coef[3]+fitlr$coef[4]) - exp(fitlr$coef[2]) -exp(fitlr$coef[3]) + 1}5 d4 W6 T/ s1 @$ P
Last expression in function is an assignment4 e2 M1 |: y& N. W+ R3 L) M
(You probably wanted to return the left-hand side)
/ g5 J/ j' M% n in: reri <- exp(fitlr$coef[2] + fitlr$coef[3] + fitlr$coef[4]) - exp(fitlr$coef[2]) - exp(fitlr$coef[3]) + 1; [3 l( o1 Z. v% v" X% S* L
> summary.bootstrap(bootstrap(smoking, reri(smoking), B=10000), probs=c(0.025,0.5, 0.975))
! X6 J0 ^+ |+ WWarning messages:# I) x, W) v$ |
Numerical expression has 7028 elements: only the first used in: model.frame(formula, data, na.action, dots)
( {' v, B6 d2 A& N: \, @; R7 D0 RProblem: Object "DM.0813" not found : }9 r. w: G9 }' c
Use traceback() to see the call stack
0 T4 E9 E, w4 H1 p
' N7 y( R; g+ C
o C. }4 H0 r7 ]不明白这个错误什么意思,小白求指导……6 K. G3 ?1 s" d
" V- n* N+ `! g" S4 k
1 i7 b! B' A# S
|
zan
|