全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学 R语言论坛
2770 1
2016-03-10
我用coxph做了一个回归,算出来beta为(-0.383,-0.302),然后我写了一个cox的似然函数,将beta的值带入始终算不出来coxph$loglik里面beta=final的那个loglik值。我用beta=init带入我的程序可以得出beta=init的那个loglik值。求大神指点是否程序有问题。
复制代码


二维码

扫码加我 拉你入群

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

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

全部回复
2016-3-10 14:01:06
另外,我还参照coxph源代码里的程序写了一个算loglik的程序
复制代码
复制代码
上面两段程序跑出来的结果是一样,但是第一段程序zbeta用的是beta*x算的,而第二段程序zbeta用的是beta*(x-mean(x))算的,为什么这两段程序跑出来的结果会是一样的呢?
二维码

扫码加我 拉你入群

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

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

相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群