- 在线时间
- 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 编辑 7 S9 @$ \: e/ E; w/ J
0 x0 C% B! {. R! h6 { o1 ^ x
这是一篇文献里给出的script,但我放到S-plus里跑不出来,有错误,有可能是我不会用S-plus,所以请教大神下面的script是否有错,没有错的话怎么解决错误?能不能改动换到R里运行?
3 k3 i+ U, d2 I( v9 g" {
T; S {- j% G# v. G1 z7 ilibrary(Design)
9 \' S! }: _: o- f K# ?+ treri <- function(datsam) - |! X _* O; W9 W: O# X# E
{ , s' ?* D' k1 E0 l) \! p
fitlr <- glm(outcome ~ determinant A * determinant B, family=binomial, data=datsam)
6 M& A- `# ?& y& m$ \0 ]reri <- exp(fitlr$coef[2]+fitlr$coef[3]+fitlr$coef[4]) - exp(fitlr$coef[2]) - exp(fitlr$coef[3]) + 1
$ N F' W6 W0 t} 7 m. O. i+ v: K! T
summary.bootstrap(bootstrap(dataset, reri(dataset), B=10000), probs=c(0.025,0.5, 0.975)) % o" p. i+ V- o+ h& [
9 \5 k" }4 m4 H————————————我是分割线——————————————————————————————
* Q, a1 z2 o. A- H
; W( f, B, ]$ L* I" P& W1 r$ v: S* f2 }) M
8 X! E/ ]. ~# M3 [, \5 K1 g7 M
library(Design)5 f0 Q7 F& h- _" F9 I
> reri <- function(datsam)
1 H" s- X1 |* {4 m+ {fitlr <- glm(DM.0813 ~ smoking * rs163184, famliy=binomial, data=datsam): K A0 O' ?8 G' V% {
+ reri <- exp(fitlr$coef[2]+fitlr$coef[3]+fitlr$coef[4]) - exp(fitlr$coef[2]) -exp(fitlr$coef[3]) + 1}
9 ]- g0 `& e7 ] Last expression in function is an assignment
+ ?8 f3 x& n9 r1 D9 b* c' _ (You probably wanted to return the left-hand side)! s1 Y; }5 e. B9 O, v2 M o0 W/ z
in: reri <- exp(fitlr$coef[2] + fitlr$coef[3] + fitlr$coef[4]) - exp(fitlr$coef[2]) - exp(fitlr$coef[3]) + 1
4 j0 i$ B# q0 B3 t( P> summary.bootstrap(bootstrap(smoking, reri(smoking), B=10000), probs=c(0.025,0.5, 0.975))
1 S' @3 [9 _2 tWarning messages:' V+ x0 g+ }% p9 U2 W! ~6 P
Numerical expression has 7028 elements: only the first used in: model.frame(formula, data, na.action, dots)
0 p% }' Z- I% S# B# H* [0 RProblem: Object "DM.0813" not found
+ n/ W( p$ ` w& p5 tUse traceback() to see the call stack1 y7 p9 K6 L0 v
. O4 Y. S, e9 q/ |% ?& l: A, D# E: A% k% @- j9 ?
不明白这个错误什么意思,小白求指导……; x- \: x8 w- C" c
: s7 r( N' z& d7 y% G3 [7 C. e J, N5 x, X& p. Q! s( u
|
zan
|