similar to: Isolating a single plot from plots produced simultaneously

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

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

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

-- 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

Hi all, I'm very confused! I've been using the same code for many weeks without any bother for various covariates. I'm now looking at another covaraite and whenever I run the code you can see below I get an error message: "Error in rep(0, nrow(data)) : invalid 'times' argument" This code works: # remove 'missing' cases from data # snearma <-

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 using Windows and R 2.9.2 for my analyses. I have a large dataset and I am particularly interested in looking at time to an event for a continuous variable. I would like to produce a plot of log(relative risk) or relative risk (also known as hazard ratio) against the continuous variable. I have spent a long time looking for advice on how to do this but my search has proved

Dear All, I have analysed time to event data for continuous variables by considering the multivariable fractional polynomial (MFP) model and comparing this to the untransformed and log transformed model to determine which transformation, if any, is best. This was possible as the Cox model was the underlying model. However, I am now at the situation where the assumption that the competing risks

Dear R-help, I am using R 2.14.1 on Windows 7 with the 'gfcure' package (cure rate model). I have included the treatment variable in the cure part of the model as shown below: Ø ref_treat <- gfcure(Surv(rem.Remtime,rem.Rcens)~1,~1+strata(drpa)+factor(treat(delcure)),data=delcure,dist="loglogistic") >From that I can obtain the coefficients, standard errors etc as per

Hi, at my local university network, I have to make sure no student uploads more than x GB/day. So far, I give them unlimited bandwidth until they have more than y GB (y<x) upload. After that, I use the u32 filter to associate the affected IP address with a HTB class limited to the rate (remainingtraffic/remainingtime). Since the accounting is done with ipt_ACCOUNT right now and the netfilter

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

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,

Dear R-users, I am trying to modify the mfp() function in the mfp package to model Weibull survival using fractional polynomials approach. However, I keep getting into trouble when mfp.fit and other "hidden" functions can't be found. I did find some of them in Splus but it's getting nowhere. I wonder if any of you can give me some tips on how to modify it or any experience

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 hpa, please review and on ack merge belows patchset. thanks maks git pull git://git.debian.org/~maks/klibc.git maks for the changes: maximilian attems (8): [klibc] Add initial mntent.h with setmntent() and endmntent() [klibc] add getmntent() [klibc] mount: add help arg for usage() [klibc] mount: list

I have just published a pre-alfa of a patch that has the goal to make OpenSSH aware with PMI. Reference site : http://nutmay.sourceforge.net -- Vincenzo Sciarra

[ gsi.de pretends this is massmailing so removed from CC's ] hello hpa, please review and on ack merge belows patchset. changes from v1, thanks to review by kzak: * add strtok_r() * use it in getmntent_r() * getment_r() parses mnt_freq and mnt_passno too isn't that hard. thanks maks git pull git://git.debian.org/~maks/klibc.git maks for the changes: maximilian attems (9):

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

V; Sent via BlackBerry from T-Mobile -----Original Message----- From: r-help-request at r-project.org Date: Tue, 30 Sep 2008 12:00:06 To: <r-help at r-project.org> Subject: R-help Digest, Vol 67, Issue 31 Send R-help mailing list submissions to r-help at r-project.org To subscribe or unsubscribe via the World Wide Web, visit https://stat.ethz.ch/mailman/listinfo/r-help or, via email,

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