QQ登录

只需要一步,快速开始

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

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

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

1

主题

1

听众

1

积分

升级  20%

该用户从未签到

跳转到指定楼层
1#
发表于 2018-3-16 12:27 |只看该作者 |倒序浏览
|招呼Ta 关注Ta
本帖最后由 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
转播转播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-8-25 22:24 , Processed in 0.434258 second(s), 65 queries .

回顶部