QQ登录

只需要一步,快速开始

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

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

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

1

主题

1

听众

1

积分

升级  20%

该用户从未签到

跳转到指定楼层
1#
发表于 2018-3-16 12:27 |只看该作者 |倒序浏览
|招呼Ta 关注Ta
本帖最后由 419595891 于 2018-3-16 13:51 编辑 + @) Z) B7 R8 g% z- y; R8 [, n  V

1 O+ B1 o* K$ h- J" u% W这是一篇文献里给出的script,但我放到S-plus里跑不出来,有错误,有可能是我不会用S-plus,所以请教大神下面的script是否有错,没有错的话怎么解决错误?能不能改动换到R里运行?
6 R+ [/ S+ }/ a, O. h' k' r+ J
library(Design)                                                        ; [7 ?- L, [' y7 T' E* [* z+ J" N( P
reri <- function(datsam)                                                        
; v/ b  j3 G, z7 I& W( J& y{                                                        
$ {( h5 X. ^  I2 N: c7 ffitlr <- glm(outcome ~ determinant A * determinant B, family=binomial, data=datsam)                                                        
2 |$ y# U% q0 d: |reri <- exp(fitlr$coef[2]+fitlr$coef[3]+fitlr$coef[4]) - exp(fitlr$coef[2]) - exp(fitlr$coef[3]) + 1                                                        
0 V% ~3 H4 T% s4 c}                                                        $ _5 m$ w2 d7 \. G9 d
summary.bootstrap(bootstrap(dataset, reri(dataset), B=10000), probs=c(0.025,0.5, 0.975))                                             
) X0 b) T# t0 U
; k. x% b9 [# a% R& O4 I. k————————————我是分割线——————————————————————————————       0 ~  S) d- |8 |/ y# D0 h& `

; }6 ]9 H, d$ A7 m- f$ j0 [0 n9 N4 z5 G, @2 @# j* u

8 d1 U$ v$ T; A library(Design)& S8 n( w* n0 d  ~: j
> reri <- function(datsam)
4 T, U$ B* B, x( W/ L8 N+ {fitlr <- glm(DM.0813 ~ smoking * rs163184, famliy=binomial, data=datsam)
! \$ R: e7 ?, _; b+ reri <- exp(fitlr$coef[2]+fitlr$coef[3]+fitlr$coef[4]) - exp(fitlr$coef[2]) -exp(fitlr$coef[3]) + 1}9 |* ?5 l" x( T; z' z* @
    Last expression in function is an assignment
; G; J5 W; G. V* |$ o' O                (You probably wanted to return the left-hand side)
- R) c9 j. }0 p" Z         in: reri <- exp(fitlr$coef[2] + fitlr$coef[3] + fitlr$coef[4]) - exp(fitlr$coef[2]) - exp(fitlr$coef[3]) + 1

( ]$ w" E0 Q, G/ D8 [! {> summary.bootstrap(bootstrap(smoking, reri(smoking), B=10000), probs=c(0.025,0.5, 0.975))                                               
0 R- z$ q  p! ^6 u1 n2 kWarning messages:- s) ^2 o& i( ]9 J2 ^
  Numerical expression has 7028 elements: only the first used in: model.frame(formula, data, na.action, dots): c; h. C8 ^+ R( g
Problem: Object "DM.0813" not found ' s0 O. @! W8 {' G- b; }! Q
Use traceback() to see the call stack
. I( l1 c! w/ |' X% x- }) ~
: X" U/ L0 M2 x) X

" Y7 k3 D' ~2 P. p, D$ K% |, ^' p* m$ O不明白这个错误什么意思,小白求指导……: P9 s1 j& J7 U" e1 d, t# p

$ @9 R4 y+ q) h# Y6 b  q2 k! o" N( l8 u8 `( U; T, d1 j
zan
转播转播1 分享淘帖0 分享分享0 收藏收藏0 支持支持0 反对反对0 微信微信
2#
无效楼层,该帖已经被删除
3#
无效楼层,该帖已经被删除
您需要登录后才可以回帖 登录 | 注册地址

qq
收缩
  • 电话咨询

  • 04714969085
fastpost

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

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

蒙公网安备 15010502000194号

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

GMT+8, 2026-4-10 15:05 , Processed in 0.438118 second(s), 65 queries .

回顶部