数学建模社区-数学中国
标题:
请教大神这个S-plus script是否正确？可以改动转换到R运行吗？
[打印本页]
作者:
419595891
时间:
2018-3-16 12:27
标题:
请教大神这个S-plus script是否正确？可以改动转换到R运行吗？
本帖最后由 419595891 于 2018-3-16 13:51 编辑
: l- P, p' q3 U
2 s) M. Y' z# U7 _9 @6 z2 A# E
这是一篇文献里给出的script，但我放到S-plus里跑不出来，有错误，有可能是我不会用S-plus，所以请教大神下面的script是否有错，没有错的话怎么解决错误？能不能改动换到R里运行？
! c0 X( N0 h5 w2 g( J8 n9 r% q
8 X1 s7 R3 P! f+ P
library(Design)
8 u1 Y9 i% V& x% t! ~2 q: K5 N* B3 r' q
reri <- function(datsam)
7 f: _( ^4 ]/ {+ [: y6 T+ x; a" d
{
7 f2 u7 F4 n# U* J5 m" n2 U% K
fitlr <- glm(outcome ~ determinant A * determinant B, family=binomial, data=datsam)
. v# d& k* A* L: B
reri <- exp(fitlr$coef[2]+fitlr$coef[3]+fitlr$coef[4]) - exp(fitlr$coef[2]) - exp(fitlr$coef[3]) + 1
( k5 v; }3 v' R0 ]3 _" {
}
$ V+ e: a1 n7 M. V. r5 V! s# E
summary.bootstrap(bootstrap(dataset, reri(dataset), B=10000), probs=c(0.025,0.5, 0.975))
1 Q7 D9 G' y u1 c1 t: H% i7 W
$ I+ v* q& X3 C' K+ U' V/ @
————————————我是分割线——————————————————————————————
, \& M% c# w" t+ d3 M z7 P
5 l' k9 e( D! E+ E: o& [* y
: T9 B3 ]) x5 f
" A) |. l0 i3 V t# }
library(Design)
1 X6 X% U* h# ?- X7 n
> reri <- function(datsam)
% b6 A9 f7 q4 X" v$ N. z( Y) I
+ {fitlr <- glm(DM.0813 ~ smoking * rs163184, famliy=binomial, data=datsam)
: P4 {) p; u( ~7 o; @
+ reri <- exp(fitlr$coef[2]+fitlr$coef[3]+fitlr$coef[4]) - exp(fitlr$coef[2]) -exp(fitlr$coef[3]) + 1}
' ]+ A8 }: C8 t+ i! H
Last expression in function is an assignment
$ B/ n$ t) W; r0 j# \, \' ~
(You probably wanted to return the left-hand side)
. M/ c5 q: Z$ z% P5 t4 X5 q' r" I
in: reri <- exp(fitlr$coef[2] + fitlr$coef[3] + fitlr$coef[4]) - exp(fitlr$coef[2]) - exp(fitlr$coef[3]) + 1
) V: o4 R' {' O1 ?7 U9 u
> summary.bootstrap(bootstrap(smoking, reri(smoking), B=10000), probs=c(0.025,0.5, 0.975))
% N# z' K6 Y$ R- b' Q4 }
Warning messages:
5 W/ R4 e/ s% M. J) t0 k7 Y
Numerical expression has 7028 elements: only the first used in: model.frame(formula, data, na.action, dots)
3 h4 x6 \; T" j8 E7 {
Problem: Object "DM.0813" not found
1 ^! B7 b7 i3 K
Use traceback() to see the call stack
4 t8 M3 O G% C. p
. `+ b1 ]& G( {
$ d7 O# @7 t, F( @
不明白这个错误什么意思，小白求指导……
" }+ B$ J" M' @6 I0 n- x
) [% A* V. | d% Z
. u E* B+ o7 p' z& w! R5 X8 u
欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)
Powered by Discuz! X2.5