similar to: Cox regression

This is a simple R program that I have been trying to run. I keep running into the "singular matrix" error. I end up with no sensible results. Can anyone suggest any changes or a way around this? I am a total rookie when working with R. Thanks, Rasika > library(survival) Loading required package: splines > args(coxph) function (formula, data, weights, subset, na.action, init,

Dear Everybody: I'm doing my usual "how does that work in R" thing with some Stata projects. I find a gross gap between the Stata and R in Cox PH models, and I hope you can give me some pointers about what goes wrong. I'm getting signals from R/Survival that the model just can't be estimated, but Stata spits out numbers just fine. I wonder if I should specify initial

Hi, I want to fit multiple cox models using the coxph() function. To do this, I use a for-loop and save the relevant results in a separate matrix. In the example below, only two models are fitted (my actual matrix has many more columns), one gives a warning message, while the other does not. Right now, I see all the warning message(s) after the for-loop is completed but have no idea which model

hi, I applied coxph to my matrix of 300 samples and 215 variables and got the following error Error in fitter(X, Y, strats, offset, init, control, weights = weights, : NA/NaN/Inf in foreign function call (arg 6) In addition: Warning message: Ran out of iterations and did not converge in: fitter(X, Y, strats, offset, init, control, weights = weights, 26% of time data is censored and here

Due to missing data with some of my predictor variables I first do multiple imputation as follows: library(foreign) library(Amelia) library(norm) set.seed(666) M=10 impdat <-

Dear R-help folks: Here is my platform: > version platform sparc-sun-solaris2.7 arch sparc os solaris2.7 system sparc, solaris2.7 status major 1 minor 5.0 year 2002 month 04 day 29 language R I have a

Hi, I have a script that takes a subset of genes on a microarray and tries to fit a coxph model to the expression values for each gene. This seems to work fine but in some cases it produces warnings and/or errors. For example: Error in fitter(X, Y, strats, offset, init, control, weights = weights, : NA/NaN/Inf in foreign function call (arg 6) In addition: Warning message: In fitter(X, Y,

Good day R help list!!! I've been trying to do Bootstrap in R on Censored data. I encountered WARNING/ERROR messages which I could not find explanation. I've been searching on the literature for two days now and still can't find answers. I hope there's anyone out there who can help me with these two questions: 1. If the "Loglik converged before variable..." message

Dear R-helpers: I am trying to fit a multivariate Cox proportional hazards model, modelling survival outcome as a function of treatment and receptor status. The data look like below: # structure of the data str(sample.data) List of 4 $ survobj : Surv [1:129, 1:2] 0.8925+ 1.8836+ 2.1191+ 5.3744+ 1.6099+ 5.2567 0.2081+ 0.2108+ 0.2683+ 0.4873+ ... ..- attr(*, "dimnames")=List of 2

Hello, I'm trying to estimate a cox model with a frailty variable and time-dependent covariate (below there is the statement I use and the error message). It's seems to be impossible, because every time I add the time-dependent covariate the model doesn't converge. Instead, if I estimate the same model without the time-dependent covariate it's converge. I'd like knowing if

Dear R-helpers: I am trying to fit a multivariate Cox proportional hazards model, modelling survival outcome as a function of treatment and receptor status. The data look like below: # structure of the data str(sample.data) List of 4 $ survobj : Surv [1:129, 1:2] 0.8925+ 1.8836+ 2.1191+ 5.3744+ 1.6099+ 5.2567 0.2081+ 0.2108+ 0.2683+ 0.4873+ ... ..- attr(*, "dimnames")=List of 2

Hi everyone, I encountered a problem using the coxph function for the conditional logistic regression. I am trying to do some simulations and I really don’t understand a mistake which happened maybe only 1 time among more than 1,000 simulations. What appeared on the screen is the following: Error in fitter(X,Y,strats,offset,init,control,weights=weights,: NA/NaN/Inf in a foreign function (arg

I have been attempting to use cph models to get better calibration of my models for which I had originally used logistic regression. I tried running with 40 repetitions and got an error. I then tried 500 repetitions (thinking that the NaNs in the output below might be caused by that choice) and then let my computer crunch for several hours and got only the same error message and

Hi, I am getting errors of the following kind. I can't seem to point the source of the error. I would greatly appreciate any advice. Many thanks and good day, -Melinda Error message : ---------------- "Ran out of iterations and did not converge in: fitter(X, Y, strats, offset, init, control, weights = weights,..." Details : --------- E is a vector of survival times (or censored

R list: I am using a 'for' loop to run a number of different models (stratified by different variables) with coxph. The data becomes sparse when some strata are used causing the model to become unstable. The following error occurs and the analysis is terminated. >Error in fitter(X, Y, strats, offset, init, control, weights = weights, : (converted from warning) Loglik

Hello, I have a couple of questions with regards to fitting a coxph model to a data set in R: I have a very large dataset and wanted to get the baseline hazard using the basehaz() function in the package : 'survival'. If I use all the covariates then the output from basehaz(fit), where fit is a model fit using coxph(), gives 507 unique values for the time and the corresponding cumulative

You asked about survival curves with age scale versus follow-up scale. > fit1 <- coxph(Surv(time/365.25, status) ~ t5 + id + age, data=stanford2) > surv1<- survfit(fit1) > surv1 n events median 0.95LCL 0.95UCL 157.000 102.000 1.999 0.898 3.608 > summary(surv1, times=3) time n.risk n.event survival std.err lower 95% CI upper 95% CI 3 46 85 0.451

Dear R-help members, Apologies - I am posting on behalf of a colleague, who is a little puzzled as STATA and R seem to be yielding different survival estimates for the same dataset when treating a variable as ordinal. Ordered() is used to represent an ordinal variable) I understand that R's coxph (by default) uses the Efron approximation, whereas STATA uses (by default) the Breslow. but we

hello, I just installed R-3.4.0 from scratch: $ sudo apt install r-base but when I try > library(survival, lib.loc = "/usr/lib/R/library") > fit <- coxph(Surv(exit, event) ~ x, data = mort) I get Error in fitter(X, Y, strats, offset, init, control, weights = weights, : object 'Ccoxmart' not found I was told on R-help that I need to >

Dear Sir/Ma, I Adelabu.A.A, one of the R-users from Nigeria. When am running a coxph command the below error was generated, and have try some idea but not going through. kindly please assist: > cox1 <- coxph(Surv(tmonth,status) ~ sex + age + marital + sumassure, X) Warning message: In fitter(X, Y, strats, offset, init, control, weights = weights, : Ran out of iterations and did not