similar to: Relative Risk/Hazard Ratio plots for continuous variables

Dear R-Help, I am using the 'mfp' package. It produces three plots (as I am using the Cox model) simultaneously which can be viewed together using the following code: fit <- mfp(Surv(rem.Remtime,rem.Rcens)~fp(age)+strata(rpa),family=cox,data=nearma,select=0.05,verbose=TRUE) par(mfrow=c(2,2)) plot(fit) They can be viewed separately but the return key must be pressed after each graph

All, I am struggling to create a survival plot using LTRC data for each year of a 10 year period. I have a set of individuals (birds) where 'entry' is the day of the year (1-365) they are released (let out of pens) into the wild (2 year data snip below). 'Entry' (e.g., day of year the first bird is released for each year) is highly variable, ranging from 48 to >250. When I

Hi, I am trying to convert the following NLMIXED code to NLME, but am running into problems concerning 'Singularity in backsolve'. As I am new to R/S-Plus, I thought I may be missing something in the NLME code. NLMIXED *********** proc nlmixed data=kidney.kidney; parms delta=0.03 gamma=1.1 b1=-0.003 b2=-1.2 b3=0.09 b4=0.35 b5=-1.43 varu=0.5; eta=b1*age+b2*sex+b3*gn+b4*an+b5*pkn+u;

Dear list: working with date/times I have come across a problem that ISOdate and ISOdatetime are too slow on large vectors of data. I was surprised just until I looked at the implementation and the man page: "ISOdatetime and ISOdate are convenience wrappers for strptime". In other terms, they convert data to character representation first in order to create a POSIXlt object that is then

Hi All, This sounds a relatively simple query, and I hope it is! I am looking at a continuous variable, age. I am looking at time to 12-month remission and can calculate the HR and 95% confidence interval are follows: coxfita = coxph(Surv(rem.Remtime,rem.Rcens)~nearma$all.age,data=nearma) exp(coxfita$coefficients) exp(confint(coxfita)) However, because I am looking at age as a continuous

Hi, I'm working with timeseries data. The values are every 5 seconds and each series can last up to 4-5 days. To generate the x-axis labels, I'm doing the following: ========================= # Variable for displaying hours on the x-axis rtime <<- as.POSIXct(round(range(timedata), "hours")) # Variable for displaying days on the x-axis stime <<-

Hi there, Here is a record of an R session: > library(survival4) Autoloading required package: splines > coxph(Surv(c(1,2,3,4), c(1,1,1,1)) ~ c(1,0,1,0)) Error in .C("coxfit2", iter = as.integer(iter.max), as.integer(n), as.integer(nvar), : C/Fortran function not in load table Why is "coxfit2" not in the load table? How do I fix it? System is Linux. Goran

Dear All, Sorry to bother you again. I have a model: coxfita=coxph(Surv(rem.Remtime/365,rem.Rcens)~all.sex,data=nearma) and I'm trying to do a plot of Schoenfeld residuals using the code: plot(cox.zph(coxfita)) abline(h=0,lty=3) The error message I get is: Error in plot.window(...) : need finite 'ylim' values In addition: Warning messages: 1: In sqrt(x$var[i, i] * seval) : NaNs

Dear All, I am contemplating centering the covariates in my Cox model to reduce multicollinearity between the predictors and the interaction term and to render a more meaningful interpretation of the regression coefficient. Suppose I have two indicator variables, x1 and x2 which represent age categories (x1 is patients less than 16 while x2 is for patients older than 65). If I use the following

Dear all, I am facing problem with inserting the scheduled day of Observation in the dataset. In the dataset I have only  relative time (table 1) and not scheduled day of observation (day 1, 4, 8, 15, 22, 29, 36, 43). I would appreciate if any one could suggest me how to proceed. Eg: Table 1 The real dataset looks like this with Time, DV ... etc RTime   DV 1        101 4        95 8       

I need to have an AD user that need to *write* in an users LDAP field. The user case is a MFP (a set of MFP, indeed) that have RFID auth, and so need to 'register' the RFID cards. Seems to me that i have to use dsacl/samba-tool acl ds, but i don't found a way to set the property for every user. EG, assign write permission to user 'mfp' to field 'pager' for every

Hello, I am testing a new samba share (Version 4.7.1) in order to replace an old samba share (Version 3.5.6). Especially, I am facing issues with Canon MFP "scan to share". With model Canon IR1730 (SMB1 only capable), "scan to share" seems to not work with the new 4.7.1 samba share. The canon error code returns by the Canon iR1730 is #0753 with this 4.7.1 samba share. It

-- begin included message ----- I am looking at a continuous variable, age. I am looking at time to 12-month remission and can calculate the HR and 95% confidence interval are follows: coxfita = coxph(Surv(rem.Remtime,rem.Rcens)~nearma$all.age,data=nearma) exp(coxfita$coefficients) exp(confint(coxfita)) However, because I am looking at age as a continuous variable I cannot draw a Kaplan-Meier

Dear R users, The package 'mfp' that fits fractional polynomial terms to predictors. Example: data(GBSG) f <- mfp(Surv(rfst, cens) ~ fp(age, df = 4, select = 0.05) + fp(prm, df = 4, select = 0.05), family = cox, data = GBSG) print(f) To describe the association between the original predictor, eg. age and risk for different values of age I can plot it the polynomials

Mandi! Rowland Penny via samba In chel di` si favelave... > I think that is what Andrew is trying to tell you, the printer needs to > support SASL over TLS/SSL or it will never work. I don't think there is > anything you can do, but I am surprised that the print doesn't already > support it, after all, it isn't something new ;-) Mi confusion grow. ;-) As stated in my

Dear Sir, ??????????????? Please know me how to function hp 1213nf mfp printer (printer + scanner + fax + copy) in centOS 6.2. The computer shows the printer but cannot oprate. the message displayed "Missing printer driver"? etc. ??????????????? This will be very helpful to me if anyone can inform me the technique to operate this printer in centOS 6.2. ????????? Thanking you,

I am doing some coxPH model fitting and would like to have some idea about how good the fits are. Someone suggested to use Frank Harrell's C-index measure. As I understand it, a C-index > 0.5 indicates a useful model. I am probably making an error here because I am getting values less than 0.5 on real datasets. Can someone tell me where I am going wrong please ? Here is an example using

Dear All, I am attempting to use forward and/or backward selection to determine the best model for the variables I have. Unfortunately, because I am dealing with patients and every patient is receiving treatment I need to force the variable for treatment into the model. Is there a way to do this using R? (Additionally, the model is stratified by randomisation period). I know that SAS can be

I'm looking at memory efficiency for some of the survival code. The following fragment appears in coxph.fit coxfit <- .C("coxfit2", iter=as.integer(maxiter), as.integer(n), as.integer(nvar), stime, sstat, x= x[sorted,] , ... Does this make a second copy of x to pass to the routine (my

Alle, I have an HP L7780 MFP, which is supported under hplip 1.7.2 but not the current install of 1.6.7. Does anyone know of a way just to add that particular driver to 1.6.7, or will it require removing the standard 1.6.7 and doing a custom install of 1.7.2. If the latter, what would you consider the dangers/other implications of going "rogue"? Best Regards, Camron -- Camron W.