similar to: Pec function in R

Hello everyone, (This is my second question posted today on the R list). I am carrying out a competing risks analysis using the cuminc function...this takes the form: cuminc(ftime,fstatus,group) In my study, fstatus has 3 different causes of failure (1,2,3) there are also censored cases (0). "group" has two levels (0 and 1). I therefore have 6 different cumulative incidence curves:

Dear all, I would like to obtain the Brier score prediction error at different times t for an extended Cox model. Previously I have used the 'pec' function (pec{pec}) to obtain prediction error curves for standard Cox PH models but now I have data in the counting process format (I have a covariate with a time-varying effect) and it seems that the pec function does not support the counting

Hi, I am attempting to evaluate the prediction error of a coxph model that was built after feature selection with glmnet. In the preprocessing stage I used na.omit (dataset) to remove NAs. I reconstructed all my factor variables into binary variables with dummies (using model.matrix) I then used glmnet lasso to fit a cox model and select the best performing features. Then I fit a coxph model

Dear Raja, A Fine and Gray model can be fitted using the standard coxph function with weights that correct for right censoring and left truncation. Hence I guess any function that allows to perform stepwise regression with coxph should work. See e.g. my article in Biometrics https://doi.org/10.1111/j.1541-0420.2010.01420.x, or the vignette "Multi-state models and competing risks" in the

Dear Geskus, I want to develop a prediction model. I followed your paper and analysed thro' weighted coxph approach. I can develop nomogram based on the final model also. But I do not know how to do internal validation of the model and subsequently obtain calibration plot. Is it possible to use Wolbers et al Epid 2009 approach 9 (R code for internal validation and calibration) . It is

Im trying to write a controller spec, but when I run it I get an error like the one below. I found Bug6707 (http://rubyforge.org/pipermail/rspec-devel/2007-March/002147.html ), but I get the error no matter how I run the spec ( whether I use rake or not) Any ideas? I am using 0.8.2 The Error: C:/rails/InstantRails/ruby/lib/ruby/gems/1.8/gems/rspec-0.8.2/lib/spec/expectati ons/sugar.rb:14:in

Hi, I sent this in a couple of days ago and have only gotten one suggestion. Can someone read the below problem and try to help me out? I have searched everywhere for a solution to this and have tried numerous 'passwd chat' strings to no avail. Help is much appreciated! James Herschel JSquared Network Solutions (905)847-0799 james@jsquared.ca -----Original Message----- From: James

Dear List, I have a theoretical question related to epidemiological data analysis: If the treatment status (tx = 0,1) changes over time for the patients in a non-randomized cohort, is there a way to estimate the treatment effect? (i.e., after joining the study, some patients may have to wait for a period of time before receiving the treatment, i.e., the situation of patient with id == 2 for the

...something is still wrong. On the shared folders the commands "select" and "examine" work only the first time: -------%<-------------------------------------------------------- 3 select "shared/user1 at pec.inail.gov.it" * FLAGS (\Answered \Flagged \Deleted \Seen \Draft Junk) * OK [PERMANENTFLAGS (\Answered \Flagged \Deleted \Seen \Draft Junk \*)] Flags

Hi, when I generated a survfit() object, I can get number of patients at risk at various time points by using summary(): fit<-survfit(Surv(time,status)~class,data=mtdata) summary(fit) class=1 time n.risk n.event survival std.err lower 95% CI upper 95% CI 9.9 78 1 0.987 0.0127 0.963 1 41.5 77 1 0.974 0.0179 0.940 1 54.0 76

Hi, i have a problem with dovecot 2.0.13. I have one dovecot in Front-end that has only: passdb { driver = static args = proxy=y host=10.0.0.6 nopassword=y } In Back-end i have one dovecot that does auth and exports imap/pop3 ports. In dovecot's log of Front-end i see: /Apr 02 14:33:34 imap-login: Info: proxy(//user at example.com/ <mailto:amministrazionevendite at

>>>>> Pearce, Robert <Pearce.Robert at epa.gov> >>>>> on Tue, 5 Dec 2017 21:31:09 +0000 writes: > I am resubmitting this bug report but with additional information. I am running this with windows 10: w64-mingw32 with R Under development (unstable) (2017-12-04 r73829). There is no such message (as you cite below) in the R-devel archives:

I've developed an indexer which embeds a keyframe index track in Ogg files. It embeds the index in its own track, so that players that don't understand or don't want to use the index can just ignore it. Ogg needs this to make seeking over networks faster and more efficient. Currently we must do a bisection search when seeking, which usually takes aound 6 HTTP requests, give or

Hi all, I just uploaded a new version of the bpcp package. It calculates confidence intervals for a survival distribution for right-censored data using the newly developed beta product confidence procedure. Previously developed methods can have substantial error rate inflation for the lower limit, especially at the right end of the curves when there are small numbers of events. The bpcp method

On Tue, Jun 1, 2010 at 12:56 PM, Monty Montgomery <monty at xiph.org> wrote: > On Mon, May 31, 2010 at 10:51 PM, Timothy B. Terriberry > <tterribe at email.unc.edu> wrote: >> Chris Pearce wrote: >>> ? Hi Guys & Gals, >>> >>> I need you guys to decide whether we want to include extra granulepos >>> fields to Skeleton 4. Given the

Hello, I have a question concerning fitting a cox model with a random intercept, also known as a frailty model. I am using both the coxme package, and the frailty statement in coxph. Often 'shared' frailty models are implemented in practice, to group people who are from a cluster to account for homogeneity in outcomes for people from the same cluster. I am more interested in the classic

Robert, To expand a bit on Goran's point, is it possible that you have an old library with packages installed that your R-devel is hitting. The R headers changes relatively recently, so any packages with compiled code that were built long enough ago (I don't have the exact date off the top of my head) need to be rebuilt before they can be safely used. If this is happening in a state

Hi, when ploting Kaplan-Meier estimate curves as below, the censoring symbols (crosses) to not change thickness along the lines plot(survfit(surv ~ I(x>=cut.off) ),lty=c(1,2), lwd=2) is there any strightforward way to make it happen? thanks robert -- View this message in context: http://r.789695.n4.nabble.com/censoring-symbols-on-survfit-plot-tp3311283p3311283.html Sent from the R help

On 31 May 2010 20:51, Chris Pearce <chris at pearce.org.nz> wrote: > Ok, thanks Silvia. I'll keep working on OggIndex/Skeleton 4.0 without > the new granulepos fields, and if they're ready in time I'll include > them, otherwise they can wait until Skeleton 4.x. I was waiting for Monty so summarize his ideas too, but from the irc discussion, it sounds like the extra

I'm hoping someone can suggest a solution for this interesting problem I've run across.... I work within a multi-Unix environment. We've got HP-UX, Solaris, Linux, and Tru64 (!) here. Because of this, some of our users have multiple logins, so that they can have multiple user environments for the product development scripts that they run. One of our Sun boxes serves as our Samba