R程序求助

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

M<-500
K<-100
n<-2
N<-n*K
alfa<-0.25
delta<-0.2  
beta<-0
stoptime<-2
condest<-matrix(0,M,5)
condsc<-matrix(0,M,2)
for(i in 1:M){  
print(i)

x<-Simulationposstgl(K,alfa,beta)  
clid<-as.vector(t(matrix(rep((1:K),2),K,2)))
time<-as.vector(t(x[,,1]))
status<-as.vector(t(x[,,2]))
z<-as.vector(t(x[,,3]))
s.tider<-sort(time[status[]==1])
k<-length(s.tider)
condest[i,5]<-coxph(Surv(time,status)~z)$coef
condest[i,4]<-coxph(Surv(time,status)~z+strata(clid))$coef  
condest[i,3]<-condest[i,5]/condest[i,4]
beta.I<-condest[i,5]
Lam<-cumsum(breslow.I(beta.I))
delta<-(s.tider[exp(-Lam)<=(1-0.05)])[1]
print(delta)
res<-c(0,0)
V<-1:4
U.c.res<-U.c(beta/alfa,alfa)
condsc[i,]<-as.vector(U.c.res[[4]])
condest[i,(1:2)]<-as.vector(NewR(beta/alfa,alfa))
print(condest[i,])
print(condsc[i,])
}
这是一段R程序，但是不知道为什么实现不了，还请高手帮忙解决，不胜感激

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：R程序 求助 程序

> M<-500
> K<-100
> n<-2
> N<-n*K
> alfa<-0.25
> delta<-0.2  
> beta<-0
> stoptime<-2
> condest<-matrix(0,M,5)
> condsc<-matrix(0,M,2)
> for(i in 1:M){  
+ print(i)
+
+ x<-Simulationposstgl(K,alfa,beta)  
+ clid<-as.vector(t(matrix(rep((1:K),2),K,2)))
+ time<-as.vector(t(x[,,1]))
+ status<-as.vector(t(x[,,2]))
+ z<-as.vector(t(x[,,3]))
+ s.tider<-sort(time[status[]==1])
+ k<-length(s.tider)
+ condest[i,5]<-coxph(Surv(time,status)~z)$coef
+ condest[i,4]<-coxph(Surv(time,status)~z+strata(clid))$coef  
+ condest[i,3]<-condest[i,5]/condest[i,4]
+ beta.I<-condest[i,5]
+ Lam<-cumsum(breslow.I(beta.I))
+ delta<-(s.tider[exp(-Lam)<=(1-0.05)])[1]
+ print(delta)
+ res<-c(0,0)
+ V<-1:4
+ U.c.res<-U.c(beta/alfa,alfa)
+ condsc[i,]<-as.vector(U.c.res[[4]])
+ condest[i,(1:2)]<-as.vector(NewR(beta/alfa,alfa))
+ print(condest[i,])
+ print(condsc[i,])
+ }
[1] 1
错误: 没有"Simulationposstgl"这个函数
应该是少了个什么包吧？

> M<-500
> K<-100
> n<-2
> N<-n*K
> alfa<-0.25
> delta<-0.2  
> beta<-0
> stoptime<-2
> condest<-matrix(0,M,5)
> condsc<-matrix(0,M,2)
> for(i in 1:M)  
+
+ x<-Simulationposstgl(K,alfa,beta)  
错误: 没有"Simulationposstgl"这个函数
> clid<-as.vector(t(matrix(rep((1:K),2),K,2)))
> time<-as.vector(t(x[,,1]))
错误于x[, , 1] : 量度数目不对
> status<-as.vector(t(x[,,2]))
错误于x[, , 2] : 量度数目不对
> z<-as.vector(t(x[,,3]))
错误于x[, , 3] : 量度数目不对
> s.tider<-sort(time[status[]==1])
错误于sort(time[status[] == 1]) : 找不到对象'status'
> k<-length(s.tider)
错误: 找不到对象's.tider'
> condest[i,5]<-coxph(Surv(time,status)~z)$coef
错误: 没有"coxph"这个函数
> condest[i,4]<-coxph(Surv(time,status)~z+strata(clid))$coef  
错误: 没有"coxph"这个函数
> condest[i,3]<-condest[i,5]/condest[i,4]
> beta.I<-condest[i,5]
> Lam<-cumsum(breslow.I(beta.I))
错误: 没有"breslow.I"这个函数
> delta<-(s.tider[exp(-Lam)<=(1-0.05)])[1]
错误: 找不到对象's.tider'
> print(delta)
[1] 0.2
> res<-c(0,0)
> V<-1:4
> U.c.res<-U.c(beta/alfa,alfa)
错误: 没有"U.c"这个函数
> condsc[i,]<-as.vector(U.c.res[[4]])
错误于as.vector(U.c.res[[4]]) : 找不到对象'U.c.res'
> condest[i,(1:2)]<-as.vector(NewR(beta/alfa,alfa))
错误于as.vector(NewR(beta/alfa, alfa)) : 没有"NewR"这个函数
> print(condest[i,])
[1]   0   0 NaN   0   0
> print(condsc[i,])
[1] 0 0

少了个什么包吧？希望楼主提供文章等更多的信息。

4# zhangtao


我也是想知道用什么软件包的，是关于positive stable frailty model的一个模拟程序，但是不知道用的是哪一个软件包

google搜不到这个函数

6# qoiqpwqr


恩，是呢，我也是没找见才会请教大家的，希望高手帮忙解决哦

找一下这个程序是谁写的，到他主页看看
如果R里面找不到Simulationposstgl，肯定是自编的

8# ywh19860616
谢谢诶，应该是自编的

楼主从哪儿搞到这一段程序的？