R - survdiff'in çıktısını almaya çalıştığım Sweave

Question

benim lateks pdf dosyasında çıkıyorum. R, Sweave ve latx kullanarak survdiff'in pdf formatında çıktısını almaya çalışıyorum.

Ben aşağıdaki alakalı tüm kod taklit ediyorum, ama bu çıktısını gerek parçasıdır: Burada

 N Observed Expected (O-E)^2/E (O-E)^2/V 
male=0 10  3  8.57  3.624  4.67 
male=1 65  48 42.43  0.732  4.67 

Chisq= 4.7 on 1 degrees of freedom, p= 0.0306 

istediğim neredeyse beni nerede alır kod şudur:

s.male <- survfit(Surv(futime, status)~ strata(male), shou3, conf.type="log-log") 
plot(s.male, lty=c(1,3), xlab="Time", ylab="Survival Probability") 
legend(30, 1.1, c("Female","Male") 
     , lty=c(1,3)) 
sma <- print.data.frame(survdiff(Surv(futime, status) ~ male 
     , data=shou3, rho=0)) 
sma2 <- print.data.frame(sma, digits = 2) 

> sma2 
Call: 
survdiff(formula = Surv(futime, status) ~ male, data = shou3, 
    rho = 0) 

     N Observed Expected (O-E)^2/E (O-E)^2/V 
male=0 10  3  8.57  3.624  4.67 
male=1 65  48 42.43  0.732  4.67 

Chisq= 4.7 on 1 degrees of freedom, p= 0.0306 

İşte her şey bir şey yapmıyor deniyorum.

<<results=tex, echo=FALSE>>= 
latex(print(sma), file='', landscape=TRUE, round=1) 
#latex(sma2, file = " ") 
#latex(print.data.frame(sma2, cdec = 2, file='')) 
#malesd <- survdiff(as.formula(sma$call$formula)) 
#latex(print(sma2, latex=TRUE, coefs=FALSE, file = " ")) 
#latex(sma2, file = " ") 
#  , cdec=c(0,0,2,2,2), na.blank=TRUE, booktabs=TRUE, table.env=FALSE, center="none", file="", title="") 
@ 

I (OE) elde etmek gibi olamaz

^2/D (OE)^2/V nesnenin üzerinden elemanları, ve dışarı Chisq hattı alamayan, istediğim (Cbind/rbind denedim). R Studio kullanıyorum, bir .rnw dosyası. Bu iyi biçimlendirilmemişse özür dilerim. Ben ABD'de sorabilir meslektaşları için bir zaman farkı ile vurgulayarak bir son tarih & karşı, bu yüzden HERHANGİ yardım çok takdir, teşekkür ederim! -Jen

Answer 1

İki çözümüm var. Ne de çok zarif.

Bazı örnek veriler.

n <- 1000 
set.seed(731) 
age <- 50 + 12*rnorm(n) 
label(age) <- "Age" 
sex <- factor(sample(c('Male','Female'), n, 
      rep=TRUE, prob=c(.6, .4))) 
cens <- 15*runif(n) 
h <- .02*exp(.04*(age-50)+.8*(sex=='Female')) 
dt <- -log(runif(n))/h 
label(dt) <- 'Follow-up Time' 
status <- ifelse(dt <= cens, 1, 0) 
futime <- pmin(dt, cens) 
# Check it out: 
survfitObj <- survfit(Surv(futime, status)~ strata(sex), conf.type="log-log") 
plot(survfitObj, lty=c(1,3), xlab="Time", ylab="Survival Probability") 
survdiffObj <- survdiff(Surv(futime, status) ~ sex, rho=0) 

---

Çözüm 1:

<<results=verbatim, echo=FALSE>>= 
survdiffObj 
@ 

---

Çözüm 2: Yerine çıkışına lateks sonuçları çalışmakla, sadece ham çıkışını kullanmak "elle"

tablo yapın

Bu, yine de biçimlendirilmelidir.