QQ登录

只需要一步,快速开始

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

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

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

1

主题

1

听众

1

积分

升级  20%

该用户从未签到

跳转到指定楼层
1#
发表于 2018-3-16 12:27 |只看该作者 |倒序浏览
|招呼Ta 关注Ta
本帖最后由 419595891 于 2018-3-16 13:51 编辑 7 S9 @$ \: e/ E; w/ J
0 x0 C% B! {. R! h6 {  o1 ^  x
这是一篇文献里给出的script,但我放到S-plus里跑不出来,有错误,有可能是我不会用S-plus,所以请教大神下面的script是否有错,没有错的话怎么解决错误?能不能改动换到R里运行?
3 k3 i+ U, d2 I( v9 g" {
  T; S  {- j% G# v. G1 z7 ilibrary(Design)                                                        
9 \' S! }: _: o- f  K# ?+ treri <- function(datsam)                                                        - |! X  _* O; W9 W: O# X# E
{                                                        , s' ?* D' k1 E0 l) \! p
fitlr <- glm(outcome ~ determinant A * determinant B, family=binomial, data=datsam)                                                        
6 M& A- `# ?& y& m$ \0 ]reri <- exp(fitlr$coef[2]+fitlr$coef[3]+fitlr$coef[4]) - exp(fitlr$coef[2]) - exp(fitlr$coef[3]) + 1                                                        
$ N  F' W6 W0 t}                                                        7 m. O. i+ v: K! T
summary.bootstrap(bootstrap(dataset, reri(dataset), B=10000), probs=c(0.025,0.5, 0.975))                                              % o" p. i+ V- o+ h& [

9 \5 k" }4 m4 H————————————我是分割线——————————————————————————————      
* Q, a1 z2 o. A- H
; W( f, B, ]$ L* I" P& W1 r$ v: S* f2 }) M
8 X! E/ ]. ~# M3 [, \5 K1 g7 M
library(Design)5 f0 Q7 F& h- _" F9 I
> reri <- function(datsam)
1 H" s- X1 |* {4 m+ {fitlr <- glm(DM.0813 ~ smoking * rs163184, famliy=binomial, data=datsam): K  A0 O' ?8 G' V% {
+ reri <- exp(fitlr$coef[2]+fitlr$coef[3]+fitlr$coef[4]) - exp(fitlr$coef[2]) -exp(fitlr$coef[3]) + 1}
9 ]- g0 `& e7 ]    Last expression in function is an assignment
+ ?8 f3 x& n9 r1 D9 b* c' _                (You probably wanted to return the left-hand side)! s1 Y; }5 e. B9 O, v2 M  o0 W/ z
         in: reri <- exp(fitlr$coef[2] + fitlr$coef[3] + fitlr$coef[4]) - exp(fitlr$coef[2]) - exp(fitlr$coef[3]) + 1

4 j0 i$ B# q0 B3 t( P> summary.bootstrap(bootstrap(smoking, reri(smoking), B=10000), probs=c(0.025,0.5, 0.975))                                               
1 S' @3 [9 _2 tWarning messages:' V+ x0 g+ }% p9 U2 W! ~6 P
  Numerical expression has 7028 elements: only the first used in: model.frame(formula, data, na.action, dots)
0 p% }' Z- I% S# B# H* [0 RProblem: Object "DM.0813" not found
+ n/ W( p$ `  w& p5 tUse traceback() to see the call stack1 y7 p9 K6 L0 v

. O4 Y. S, e9 q/ |% ?& l: A, D# E: A% k% @- j9 ?
不明白这个错误什么意思,小白求指导……; x- \: x8 w- C" c

: s7 r( N' z& d7 y% G3 [7 C. e  J, N5 x, X& p. Q! s( u
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-10-23 06:17 , Processed in 0.721437 second(s), 61 queries .

回顶部