我想用R中的一个预测器来估计指数危险模型.出于某种原因,当我使用带有偏移log t的glm poisson并且当我只使用生存包中的幸存函数时,我得到具有相反符号的系数.我确信这个解释非常明显,但我无法弄清楚.
例
t <- c(89,74,23,74,53,3,177,44,28,43,25,24,31,111,57,20,19,137,45,48,9,17,4,59,7,26,180,56,36,51,6,71,23,6,13,28,16,180,16,25,6,25,4,5,32,94,106,1,69,63,31)
d <- c(0,1,1,0,1,1,0,1,1,0,1,1,1,1,0,0,1,0,1,1,1,0,1,0,1,1,0,0,1,1,1,1,1,1,1,1,1,0,1,1,1,1,1,1,1,0,1,1,1,1,1)
p <- c(1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,0,0,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,0,0,0,1,1,1)
df <- data.frame(d,t,p)
# exponential hazards model using poisson with offest log(t)
summary(glm(d ~ offset(log(t)) + p, data = df, family = "poisson"))
生产:
Coefficients:
Estimate Std. Error z value Pr(>|z|)
(Intercept) -5.3868 0.7070 -7.619 2.56e-14 ***
p 1.3932 0.7264 1.918 0.0551 .
相比
# exponential hazards model using survreg exponential
require(survival)
summary(survreg(Surv(t,d) ~ p, data = df, dist = "exponential"))
生产:
Value Std. Error z p …我想使用 knit 和（理想情况下）stargazer 自动生成带有变量组标题的描述性表格。由于我需要加权描述，因此我不使用 stargazer 的内置汇总函数，而是生成包含统计数据的数据框并使用summary=FALSE来打印数据帧。
问题 1：将变量和标题作为行、将汇总统计数据作为列的 df 不起作用，因为 stargazer 将NA标题行上的 s 转换为$$s，这会破坏编织过程。
问题 2：作为解决方法，我生成了一个数据框，其中变量和标题作为列，汇总统计数据作为行，并使用参数flip=TRUE在观星仪输出中翻转行和列。虽然这允许我为标题使用空字符向量，为变量使用数字向量，但 stargazer 不会以数学模式输出数字向量，而是（似乎）将它们视为字符。
例子：
# create example df
df <- data.frame(heading=c(" "," "," "),var1=c(1,2,3),var2=c(4,5,6))
df$heading <- as.character(df$heading)
# output using stargazer
stargazer(df, summary = FALSE, flip = TRUE)
% Table created by stargazer v.5.2 by Marek Hlavac, Harvard University. E-mail: hlavac at fas.harvard.edu
% Date and time: Fri, Aug 12, 2016 - 10:39:01
\begin{table}[!htbp] \centering
\caption{}
\label{} …