- 在线时间
- 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 编辑
1 A$ i% V6 N) F0 @: q1 m. y
3 o1 o+ B2 @. n: }' s& F这是一篇文献里给出的script,但我放到S-plus里跑不出来,有错误,有可能是我不会用S-plus,所以请教大神下面的script是否有错,没有错的话怎么解决错误?能不能改动换到R里运行?
/ k* C9 w6 u+ R& Y7 }1 z1 Z" Q6 ?* U9 n. N
library(Design) ~$ ~- \' h |7 C
reri <- function(datsam) 4 |' y! b- N- R% `+ J; r+ B
{ 0 m# v0 L' `1 [" n6 e6 P1 N
fitlr <- glm(outcome ~ determinant A * determinant B, family=binomial, data=datsam) 7 B' j/ B' K2 R5 X6 L# }) M9 n
reri <- exp(fitlr$coef[2]+fitlr$coef[3]+fitlr$coef[4]) - exp(fitlr$coef[2]) - exp(fitlr$coef[3]) + 1
: q. C0 z" I( P4 O" v l} # d, A$ e+ x6 n
summary.bootstrap(bootstrap(dataset, reri(dataset), B=10000), probs=c(0.025,0.5, 0.975))
1 T# D2 Z6 y4 N. x- W2 H1 R; t
! R$ Y( \) y; B) N( r( `9 }9 L————————————我是分割线——————————————————————————————
4 ?% r' G9 t" \1 D
$ R* y1 \0 `* P. q$ ]' o6 j
7 Y/ V% Q4 {! d0 l! O
# | J) h, @, R I3 X! w library(Design)) a+ U6 \! Y, S7 n
> reri <- function(datsam)
* G6 p, _, G5 T+ {fitlr <- glm(DM.0813 ~ smoking * rs163184, famliy=binomial, data=datsam)
! I$ v. j P0 O) O, m0 v+ reri <- exp(fitlr$coef[2]+fitlr$coef[3]+fitlr$coef[4]) - exp(fitlr$coef[2]) -exp(fitlr$coef[3]) + 1}) B2 L3 i( u0 q5 G0 _) `
Last expression in function is an assignment
1 f. V$ O b8 i9 t: s (You probably wanted to return the left-hand side)
+ F" \/ d5 k9 @9 S5 Z in: reri <- exp(fitlr$coef[2] + fitlr$coef[3] + fitlr$coef[4]) - exp(fitlr$coef[2]) - exp(fitlr$coef[3]) + 1" |, D; E3 X$ w6 Y1 {
> summary.bootstrap(bootstrap(smoking, reri(smoking), B=10000), probs=c(0.025,0.5, 0.975))
, b! j: _1 H# z3 U7 LWarning messages:* t3 P6 D( k* h# S5 e& a
Numerical expression has 7028 elements: only the first used in: model.frame(formula, data, na.action, dots)
! r* _8 T; `! W3 GProblem: Object "DM.0813" not found
! j v( L/ S8 T l. s& oUse traceback() to see the call stack8 r! O9 _1 u% M3 I: P1 D
; @) N% X0 W0 Z" p
8 f% R& l# b" l) {# K# b
不明白这个错误什么意思,小白求指导……, Q8 l* C0 b* W# p, N9 U
- m% i' M" i) q! _) O6 u
/ z7 o. j3 J; z3 R |
zan
|