QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 4459|回复: 0
打印 上一主题 下一主题

请教大神这个S-plus script是否正确?可以改动转换到R运行吗?

[复制链接]
字体大小: 正常 放大
419595891        

1

主题

1

听众

1

积分

升级  20%

该用户从未签到

跳转到指定楼层
1#
发表于 2018-3-16 12:27 |只看该作者 |倒序浏览
|招呼Ta 关注Ta
本帖最后由 419595891 于 2018-3-16 13:51 编辑
8 A* W0 m7 D# X  L" f# k1 ?; i. [  |; g' [! N
这是一篇文献里给出的script,但我放到S-plus里跑不出来,有错误,有可能是我不会用S-plus,所以请教大神下面的script是否有错,没有错的话怎么解决错误?能不能改动换到R里运行?
; S: O1 `6 B; H: P* F7 L$ I/ c$ m
7 r; a' m: a" J' u/ S8 c: vlibrary(Design)                                                        1 \2 w+ H- \. Q1 W0 T/ ]- ^9 {* u4 w) ]5 U
reri <- function(datsam)                                                        
5 `" ~4 [# F: @. j& \{                                                        ; D& M) Q/ ]1 G& {5 N! Q
fitlr <- glm(outcome ~ determinant A * determinant B, family=binomial, data=datsam)                                                        
1 {, E3 |: b/ Z# S* u6 L0 jreri <- exp(fitlr$coef[2]+fitlr$coef[3]+fitlr$coef[4]) - exp(fitlr$coef[2]) - exp(fitlr$coef[3]) + 1                                                        * t1 k5 v3 j8 E2 p8 m  y) Z
}                                                        
5 D# A, r$ Z+ H# x6 }summary.bootstrap(bootstrap(dataset, reri(dataset), B=10000), probs=c(0.025,0.5, 0.975))                                              ( O3 R. U  V5 ^6 U0 P
4 P0 Y, h4 g; ], M
————————————我是分割线——————————————————————————————       ! w4 f$ [7 P" a/ f/ \; }& b. b# y
; ?3 s! r, J4 N8 _

# {) T; R7 @' ?7 R
- D! I$ h) x! {) R library(Design)* z8 o9 x+ q' h* O; J$ J
> reri <- function(datsam)
3 p9 V3 R1 {3 |) y9 w9 T$ B+ m+ {fitlr <- glm(DM.0813 ~ smoking * rs163184, famliy=binomial, data=datsam)/ N% k: R, d& g; A& d2 Z
+ reri <- exp(fitlr$coef[2]+fitlr$coef[3]+fitlr$coef[4]) - exp(fitlr$coef[2]) -exp(fitlr$coef[3]) + 1}5 d4 W6 T/ s1 @$ P
    Last expression in function is an assignment4 e2 M1 |: y& N. W+ R3 L) M
                (You probably wanted to return the left-hand side)
/ g5 J/ j' M% n         in: reri <- exp(fitlr$coef[2] + fitlr$coef[3] + fitlr$coef[4]) - exp(fitlr$coef[2]) - exp(fitlr$coef[3]) + 1
; [3 l( o1 Z. v% v" X% S* L
> summary.bootstrap(bootstrap(smoking, reri(smoking), B=10000), probs=c(0.025,0.5, 0.975))                                               
! X6 J0 ^+ |+ WWarning messages:# I) x, W) v$ |
  Numerical expression has 7028 elements: only the first used in: model.frame(formula, data, na.action, dots)
( {' v, B6 d2 A& N: \, @; R7 D0 RProblem: Object "DM.0813" not found : }9 r. w: G9 }' c
Use traceback() to see the call stack
0 T4 E9 E, w4 H1 p

' N7 y( R; g+ C
  o  C. }4 H0 r7 ]不明白这个错误什么意思,小白求指导……6 K. G3 ?1 s" d
" V- n* N+ `! g" S4 k
1 i7 b! B' A# S
zan
转播转播1 分享淘帖0 分享分享0 收藏收藏0 支持支持0 反对反对0 微信微信
您需要登录后才可以回帖 登录 | 注册地址

qq
收缩
  • 电话咨询

  • 04714969085
fastpost

关于我们| 联系我们| 诚征英才| 对外合作| 产品服务| QQ

手机版|Archiver| |繁體中文 手机客户端  

蒙公网安备 15010502000194号

Powered by Discuz! X2.5   © 2001-2013 数学建模网-数学中国 ( 蒙ICP备14002410号-3 蒙BBS备-0002号 )     论坛法律顾问:王兆丰

GMT+8, 2026-7-23 02:10 , Processed in 0.410498 second(s), 56 queries .

回顶部