- 在线时间
- 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 编辑
/ ]! o& ~7 Z" @8 }8 y3 y* E* i6 F; m9 k4 L# J! N/ }1 A4 F/ Q/ v0 H
这是一篇文献里给出的script,但我放到S-plus里跑不出来,有错误,有可能是我不会用S-plus,所以请教大神下面的script是否有错,没有错的话怎么解决错误?能不能改动换到R里运行?
6 H; o7 K+ X8 r6 k' M6 A" \, z( n8 r
library(Design) 6 O; F9 _2 G6 y$ R; K
reri <- function(datsam) % W: N% y( Y+ T$ P6 o+ u
{
. F1 c" V/ ], j5 wfitlr <- glm(outcome ~ determinant A * determinant B, family=binomial, data=datsam) . g9 M# ?. s1 {: V
reri <- exp(fitlr$coef[2]+fitlr$coef[3]+fitlr$coef[4]) - exp(fitlr$coef[2]) - exp(fitlr$coef[3]) + 1
4 w% o t0 X; b/ T2 |9 k}
) n5 K! `' N5 f2 E! {# Qsummary.bootstrap(bootstrap(dataset, reri(dataset), B=10000), probs=c(0.025,0.5, 0.975))
: y0 w' x/ K% S1 f! j: m' J& z' B3 [. A8 R0 r$ \7 G# M
————————————我是分割线——————————————————————————————
( l% G* V, a1 b) R8 @ z$ z
) p' r% D+ }% ~4 M- R( w6 O8 @) n( P7 w. C
* J6 r2 }8 m) ~( z. u _, |/ |
library(Design)0 ~: L! ?7 R/ H. X- I) E
> reri <- function(datsam)7 C/ h3 ^. z" ~
+ {fitlr <- glm(DM.0813 ~ smoking * rs163184, famliy=binomial, data=datsam)
2 X+ B8 x1 b3 \ K; N+ reri <- exp(fitlr$coef[2]+fitlr$coef[3]+fitlr$coef[4]) - exp(fitlr$coef[2]) -exp(fitlr$coef[3]) + 1}/ n( t# I0 U* i6 r
Last expression in function is an assignment7 [+ ]% v( m8 m }+ A7 U
(You probably wanted to return the left-hand side)
$ f$ i$ k( y8 ?$ d7 T in: reri <- exp(fitlr$coef[2] + fitlr$coef[3] + fitlr$coef[4]) - exp(fitlr$coef[2]) - exp(fitlr$coef[3]) + 17 `7 u: `" o4 r3 p1 c( L b7 k6 X6 f. q
> summary.bootstrap(bootstrap(smoking, reri(smoking), B=10000), probs=c(0.025,0.5, 0.975))
/ w3 J; N; Q! }, s. S8 l9 E7 S5 kWarning messages:3 T f, z9 t4 ], [) y) V
Numerical expression has 7028 elements: only the first used in: model.frame(formula, data, na.action, dots)
/ r1 Q/ w! w- w: X: K7 w0 ?Problem: Object "DM.0813" not found
4 t3 p5 m: f' bUse traceback() to see the call stack
% @8 R, w: u* K% ]8 B+ F" k h1 d) \8 `, v2 n
! H& y! L% }# _
不明白这个错误什么意思,小白求指导……- K2 V& @1 f- Y
* Q5 x6 x9 a: M, v7 K7 u
* I _- C+ U8 v/ d
|
zan
|