- 在线时间
- 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 编辑
2 P2 O1 b/ H- e+ B2 _9 O. n
6 Z _' K/ s9 I7 n- m% f2 J5 k; k这是一篇文献里给出的script,但我放到S-plus里跑不出来,有错误,有可能是我不会用S-plus,所以请教大神下面的script是否有错,没有错的话怎么解决错误?能不能改动换到R里运行?8 e) h2 U4 e2 C9 D. U9 d
; }1 z! e% n0 d
library(Design) 2 \7 ]! I) [+ C1 o! F% i, s
reri <- function(datsam)
% p5 G% d5 C% N1 S* T; z{ ' b }: R$ p# _* f! U# M
fitlr <- glm(outcome ~ determinant A * determinant B, family=binomial, data=datsam)
1 ]4 x' d& j' I1 Ireri <- exp(fitlr$coef[2]+fitlr$coef[3]+fitlr$coef[4]) - exp(fitlr$coef[2]) - exp(fitlr$coef[3]) + 1
9 a6 d0 l+ T! p} - P" k; p3 V/ |( \
summary.bootstrap(bootstrap(dataset, reri(dataset), B=10000), probs=c(0.025,0.5, 0.975))
9 P2 Z1 E0 N" z7 j" V$ Q# a8 L" T6 c( ^, u7 O, F# R1 q
————————————我是分割线——————————————————————————————
5 i" z4 _- R9 x
# X& _$ }+ [4 `4 d# F0 @
( q9 P8 d1 q7 H9 Q$ F# T+ O* B: ]) R9 b$ e ^# O# ~+ {& f0 [
library(Design)) b8 E4 u. @7 a4 U5 Z: O& X
> reri <- function(datsam)# H# e7 ^. Z1 S, q( ]. J4 F, M" I; b
+ {fitlr <- glm(DM.0813 ~ smoking * rs163184, famliy=binomial, data=datsam)! Y/ B) [6 J- E% K* Y" b( J. z4 `
+ reri <- exp(fitlr$coef[2]+fitlr$coef[3]+fitlr$coef[4]) - exp(fitlr$coef[2]) -exp(fitlr$coef[3]) + 1}
- B% Z {5 _. U; \ Last expression in function is an assignment/ Q5 r, g6 \& }0 N3 x
(You probably wanted to return the left-hand side)
- x0 D# \2 E# n3 y( { in: reri <- exp(fitlr$coef[2] + fitlr$coef[3] + fitlr$coef[4]) - exp(fitlr$coef[2]) - exp(fitlr$coef[3]) + 1
4 V* [9 l/ v6 N1 p4 J" b> summary.bootstrap(bootstrap(smoking, reri(smoking), B=10000), probs=c(0.025,0.5, 0.975)) : u \' h" {4 U8 g9 G
Warning messages:
. u* A. b5 c( T3 q b6 m Numerical expression has 7028 elements: only the first used in: model.frame(formula, data, na.action, dots)9 M- W# K3 s7 ? i
Problem: Object "DM.0813" not found 8 E, V+ q! {: E* O
Use traceback() to see the call stack
) ~7 n8 t- V. _1 y
4 L' F+ Z- T; `8 x# ?. ?( i! E4 o1 M
不明白这个错误什么意思,小白求指导……; W/ }: F {7 U0 \2 w/ x4 P
0 W7 @! [) R$ l* K" T. E3 S) P
6 J, b, [) h3 H# o
|
zan
|