- 在线时间
- 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 编辑
! e; X! B/ p& D9 v% G% G) q7 h8 q X# [" y
这是一篇文献里给出的script,但我放到S-plus里跑不出来,有错误,有可能是我不会用S-plus,所以请教大神下面的script是否有错,没有错的话怎么解决错误?能不能改动换到R里运行?
+ p) D; W/ ?7 S; k% \$ k5 g# [( F6 J4 Q2 U2 I0 Z | V$ ]
library(Design)
7 O' g& g# l5 j& T& \reri <- function(datsam) 7 D- n3 Q! ?% N; T- n: Q
{
5 }% g, @, s7 W+ _( T4 dfitlr <- glm(outcome ~ determinant A * determinant B, family=binomial, data=datsam)
1 ^7 @" m+ j) \# k. ireri <- exp(fitlr$coef[2]+fitlr$coef[3]+fitlr$coef[4]) - exp(fitlr$coef[2]) - exp(fitlr$coef[3]) + 1
2 k4 h: ?# b& L* r2 N8 g} 2 K# n* Q/ F% ?( L1 f- s/ E, I
summary.bootstrap(bootstrap(dataset, reri(dataset), B=10000), probs=c(0.025,0.5, 0.975))
' M3 k0 e2 ^5 K1 P1 { v' O+ Q' N8 N9 Z2 i! Z
————————————我是分割线—————————————————————————————— # T* V2 u8 _0 i3 x/ N" p
/ W# D; F7 o h" x
: o- V0 G" b6 X/ ~ n0 ~
8 s) D( ]0 @6 M2 o& t# h library(Design)
' q' S% K1 c! [) S _ p> reri <- function(datsam)
! l2 H' f: U. x% O/ k+ {fitlr <- glm(DM.0813 ~ smoking * rs163184, famliy=binomial, data=datsam)
3 X$ J7 `& I- j0 w& \+ reri <- exp(fitlr$coef[2]+fitlr$coef[3]+fitlr$coef[4]) - exp(fitlr$coef[2]) -exp(fitlr$coef[3]) + 1}, J# p# M8 [0 D; k$ w* |/ k" D
Last expression in function is an assignment
8 M. I$ D. e8 e; x- F: O7 q' ~, j- |5 Z (You probably wanted to return the left-hand side)
I( ^$ ?0 M% A1 j9 l$ l in: reri <- exp(fitlr$coef[2] + fitlr$coef[3] + fitlr$coef[4]) - exp(fitlr$coef[2]) - exp(fitlr$coef[3]) + 1
/ c1 r* o- y* i, i1 {9 F& d7 E. Y> summary.bootstrap(bootstrap(smoking, reri(smoking), B=10000), probs=c(0.025,0.5, 0.975))
" `! e: k. v/ `# rWarning messages:$ U0 F' C; J! b1 C( w3 t4 G
Numerical expression has 7028 elements: only the first used in: model.frame(formula, data, na.action, dots)2 | _6 d2 x3 G5 T+ F
Problem: Object "DM.0813" not found 3 f. H8 s6 L/ Z9 [
Use traceback() to see the call stack
6 q1 ?4 a- x F' w3 Z3 W4 m& W' D( E9 m* n4 }0 x+ U2 `6 u
& `) z: K! ~/ Y! M: p @0 e, y
不明白这个错误什么意思,小白求指导……
( ?8 \& Y& d- t' r m3 N4 }, l7 R7 c# G5 j5 @
& U# v6 }# P( Z& H
|
zan
|