Displaying 8 results from an estimated 8 matches for "survminer".
2018 Feb 14
2
Fleming-Harrington weighted log rank test
Hi all,?
The survdiff() from survival package has an argument "rho" that implements Fleming-Harrington weighted long rank test.?
But according to several sources including "survminer" package (https://cran.r-project.org/web/packages/survminer/vignettes/Specifiying_weights_in_log-rank_comparisons.html), Fleming-Harrington weighted log-rank test should have 2 parameters "p" and "q" to control the weighting for earlier vs later times in the follow-up.
For...
2018 Feb 15
0
Fleming-Harrington weighted log rank test
...2018, at 4:02 PM, array chip via R-help <r-help at r-project.org> wrote:
>
> Hi all,
>
> The survdiff() from survival package has an argument "rho" that implements Fleming-Harrington weighted long rank test.
>
> But according to several sources including "survminer" package (https://cran.r-project.org/web/packages/survminer/vignettes/Specifiying_weights_in_log-rank_comparisons.html), Fleming-Harrington weighted log-rank test should have 2 parameters "p" and "q" to control the weighting for earlier vs later times in the follow-up.
>...
2018 Feb 15
1
Fleming-Harrington weighted log rank test
...chip via R-help <r-help at r-project.org> wrote:
>>
>> Hi all,
>>
>> The survdiff() from survival package has an argument "rho" that implements Fleming-Harrington weighted long rank test.
>>
>> But according to several sources including "survminer" package (https://cran.r-project.org/web/packages/survminer/vignettes/Specifiying_weights_in_log-rank_comparisons.html), Fleming-Harrington weighted log-rank test should have 2 parameters "p" and "q" to control the weighting for earlier vs later times in the follow-up.
>...
2018 Jan 27
3
error en función ggadjustedcurves, paquete survminer
Tengo un modelo de regresión de Cox y quiero obtener el plot ajustado por una covariable (sexo) con la función ?ggadjustedcurves?, pero me da el siguiente error:
> cox2 <- coxph(os ~ imc_25 + sexo.1, data = datos)
> cox2
Call:
coxph(formula = os ~ imc_25 + sexo.1, data = datos)
coef exp(coef) se(coef) z p
imc_25 -0.621 0.537 0.299 -2.08 0.038
sexo.1M 0.714
2017 Oct 07
2
Adjusted survival curves
For adjusted survival curves I took the sample code from here:
https://rpubs.com/daspringate/survival
and adapted for my date, but got error.
I would like to understand what is my mistake. Thanks!
#ADAPTATION FOR MY DATA
library(survival)
library(survminer)
df<-read.csv("F:/R/data/base.csv", header = TRUE, sep = ";")
head(df)
ID start stop censor sex age stage treatment
1 1 0 66 0 2 1 3 1
2 2 0 18 0 1 2 4 2
3 3 0 43 1 2 3 3 1
4 4 0 47 1 2 3 NA 2
5 5 0 26 0 1 4 3 NA
S <- Surv(
time = df$start,
time2 = df$stop,
event = df$c...
2017 Oct 07
2
Adjusted survival curves
For adjusted survival curves I took the sample code from here:
https://rpubs.com/daspringate/survival
and adapted for my date, but ... have a QUESTION.
library(survival)
library(survminer)
df<-read.csv("base.csv", header = TRUE, sep = ";")
head(df)
ID start stop censor sex age stage treatment
1 1 0 66 0 2 1 3 1
2 2 0 18 0 1 2 4 2
3 3 0 43 1 2 3 3 1
4 4 0 47 1 2 3 NA 2
5 5 0 26 0 1 4 3 NA
# THANKS, DAVID WINSEMIUS for remark!!! ("R regression functions al...
2017 Oct 09
0
Adjusted survival curves
...t work, because the covariates are not factors.
2. Jeff Newmiller: "Change the columns into factors before you give them to the coxph function, e.g. df$treatment <- as.factor(df$treatment)"
And I study David Winsemius's instructions.THANKS!!!
Code works:
library(survival)
library(survminer)
df<-read.csv("F:/R/data/edgr-orig.csv", header = TRUE, sep = ";")
head(df)
# "age" means the age groups
ID start stop censor sex age stage treatment
1 0 66 0 2 1 3 1
2 0 18 0 1 2 4 2
3 0 43 1 2 3 3 1
4...
2017 Oct 09
0
Adjusted survival curves
Adjusted survival curves. (Sample code here: https://rpubs.com/daspringate/survival )
Deep gratitude?to Moderator/Admin!
At?David Winsemius prompt, more elegant working code:Thanks, Ted :)
library(survival)
library(survminer)
df<-read.csv("F:/R/data/edgr-orig.csv", header = TRUE, sep = ";")
df2 <- df
df2[,c('treatment', 'age', 'sex', 'stage')] <- lapply(df2[,c('treatment', 'age', 'sex', 'stage')], factor)
model <- coxph (S...