Dear List, Here is an example of survival data in counting process format (detailed record of each day)> data[data$Id == 11,]# extracted one person's record Id Event Fup Start Stop sex Drug1 601 11 0 6 0 1 0 0 602 11 0 6 1 2 0 0 603 11 0 6 2 3 0 0 604 11 0 6 3 4 0 0 605 11 0 6 4 5 0 1 606 11 1 6 5 6 0 1 which is compressed in the following format (unchanged records of drug exposure merged):> compressed.data[compressed.data$Id ==11,]# compressed same person's record Id Event Fup Start Stop sex Drug1 21 11 0 6 0 4 0 0 22 11 1 6 4 6 0 1 My question is: since the provided information is the same, should I expect numerically exactly same results from the following coxph outputs? If no, then which format is recommended?> data <- read.csv("http://stat.ubc.ca/~e.karim/dd.csv") > compressed.data <- read.csv("http://stat.ubc.ca/~e.karim/cd.csv") > head(data) > head(compressed.data) > coef(coxph(Surv(Start, Stop, Event) ~ sex + Drug1 + cluster(Id), robust = T, data))sex Drug1 0.8696213 3.1755854> coef(coxph(Surv(Start, Stop, Event) ~ sex + Drug1 + cluster(Id), robust = T, compressed.data))sex Drug1 0.8674742 2.7147013 PS: discrete time analogue to Cox's (using cloglog link) also gives similar results corresponding to the dataset chosen. Any suggestions/references/direction to R-package will be highly appreciated. Thanks, Ehsan
Ehsan Karim
2020-Oct-08 00:25 UTC
[R] svyglm: p-values and confidence intervals contradictory
Dear list:
I am getting p-values and confidence intervals contradictory in
svyglm() output from the survey package.
Here is a reproducible example:
https://ehsanx.github.io/SurveyDataAnalysis/#114_Regression_analysis
This problem can be easily fixed by setting appropriate df.resid in
the summary function though. But as a default option, contradictory
p-values and confidence intervals or Inf SEs is probably a bug, or I
may have missed something.
Regards,
Ehsan
https://ehsank.com/
----------------------------
packageVersion("survey")
?4.0?
R.version
_
platform x86_64-w64-mingw32
arch x86_64
os mingw32
system x86_64, mingw32
status
major 4
minor 0.2
year 2020
month 06
day 22
svn rev 78730
language R
version.string R version 4.0.2 (2020-06-22)
nickname Taking Off Again