similar to: Loop and loop output [Cox model, for, function, loglik]

Dear R Help list, My question is regarding extracting the standard error or Z-score from a cph or coxph call. My Cox model is: - modz=cph(Surv(TSURV,STATUS)~RAGE+DAGE+REG_WTIME_M+CLD_ISCH+POLY_VS, data=kidneyT,method="breslow", x=T, y=T) I've used names(modz) but can't see anything that will let me extract the Z scores for each coefficient or the standard errors in the same

Hi all, I'm wanting to be able to find and store the z-score of coxph below: - modz=coxph(Surv(TSURV,STATUS)~RAGE+DAGE+REG_WTIME_M+CLD_ISCH+POLY_VS, data=kidneyT,method="breslow") I know summary(modz) will give me this, but how do i extract the standard error or z-score values in a similar way to obtaining the coefficients by coef(modz) ? I think it must be something to do with

On 4/16/2009 10:29 AM, ?????? n/a wrote: > I want write article by russian language using Sweave. For cyrillic text > LaTeX use T2A encoding > > \usepackage[T2A]{fontenc} > > But in Sweave.sty we find: > > \RequirePackage[T1]{fontenc} > > It is source of critical problem. > > For example Rnw file > > $ cat estimation.Rnw > >

Hello R Freaks, I calculate the difference in days between two events with the following litte R expresseion: T1a <- strptime(T1,"%m/%d/%y %H:%M:%S"); T2a <- strptime(T2,"%m/%d/%y %H:%M:%S"); T1b <- as.Date(T1a); T2b <- as.Date(T2a); days <- T2b-T1b; time <- T2a - T1a; In the project I would like to calculate only working day. I the a possibility to count

I want write article by russian language using Sweave. For cyrillic text LaTeX use T2A encoding \usepackage[T2A]{fontenc} But in Sweave.sty we find: \RequirePackage[T1]{fontenc} It is source of critical problem. For example Rnw file $ cat estimation.Rnw \documentclass[A4paper]{article} \usepackage[T2A]{fontenc} \usepackage[utf8]{inputenc} \usepackage[russian,english]{babel}

Hi, I've looked through the posts but couldn't find a solution to this. I'd be really grateful if someone could help, I'd like to convert a data file of mutual information that is formatted as a matrix:             TF1    TF2    TF3    TF200... Gene1    0.0    0.2    0.2 Gene2    1.4    0.0    2.8 Gene3    0.3    0.6    1.7 Gene6000.... To a list: Gene1    TF1    0.0 Gene1   

-----BEGIN PGP SIGNED MESSAGE----- Hash: SHA1 > Message: 5 > Date: 15 Aug 2003 08:14:58 -0500 > From: "Azelton Sean (RBNA/CIT1)" <sean.azelton@us.bosch.com> > Subject: [Samba] Repost: group membership limitations and Linux kernel > To: samba@lists.samba.org > Message-ID: <1060953297.1972.1.camel@sbdgecko.sbd.us.bosch.com> > Content-Type: text/plain >

Dear R-helpers, (1) I can't figure out how to tell Dotplot that I want the colors of the CI bars to be the same as the colors of the dots. For example: t2a.ci3 <- data.frame(est = c(7, 20, 75), lower = c(-9, 0.5, 42), upper = c(22, 39, 109)) mypal <- c('skyblue3', 'mistyrose3') condNames <- c('noise', 'horn', 'marimba') Dotplot(1:2 ~

Hi, if I got it right then the survival-time we expect for a subject is the integral over the specific survival-function of the subject from 0 to t_max. If I have a trained cox-model and want to make a prediction of the survival-time for a new subject I could use survfit(coxmodel, newdata=newSubject) to estimate a new survival-function which I have to integrate thereafter. Actually I thought

Hi Chris, Your were totally right with your suggestion. I have implemented the code that : a) does not merge multiple TokenFactor nodes in the DAGCombiner::visitTokenFactor(), if the resulting TF node would contain more than 64 operands. b) produces a bunch of TokenFactor nodes with at most 64 operands, instead of one huge TokenFactor in the SelectionDAGLowering::getRoot(). If we have n

Hi Dan, Thanks for commenting on this topic. See my comments in-line. ----- Ursprüngliche Mail ---- > Von: Dan Gohman <gohman at apple.com> > An: LLVM Developers Mailing List <llvmdev at cs.uiuc.edu> > Gesendet: Mittwoch, den 30. April 2008, 21:38:26 Uhr > Betreff: Re: [LLVMdev] FoldingSetNodeID operations inefficiency > > > On Apr 28, 2008, at 6:21 AM, Roman

I hope you'll all forgive me for displaying my severe lack of knowledge on this topic, and I can't really provide much in the way of reproducible code. A colleague of mine has asked if I know how to import an S object into R. The object is stored in a file named 'pre3.f' When I open the file as a text document I get what's printed below. Being one of those young bucks with

Dear all, I've been learning biplot (Gabriel, 1971) and some days ago I sent for this list a procedural function with invitation for a collaborative package. Jari Oksanen made some suggestions and I agree with all. So, I reworked the function under the object-oriented programming (OOP/S3). I think it is now a good frame for more resources. Below it is the function and a small script to

Hello All, I am working on graph object using IGRAPH package wanted to do Bonacich Power. This is my graph object. The file 'Graph.RData' (4.2 MB) is available for download at http://dropbox.unl.edu/uploads/20090424/cfe4fcb854bb17f2/Graph.RData Graph size Vertices: 20984 Edges: 326033 Directed: FALSE No graph attributes. Vertex attributes: name. No edge attributes. When I use

On 4/18/2018 2:25 AM, Jiri Pirko wrote: > Wed, Apr 11, 2018 at 09:13:52PM CEST, sridhar.samudrala at intel.com wrote: >> On 4/11/2018 8:51 AM, Jiri Pirko wrote: >>> Tue, Apr 10, 2018 at 08:59:48PM CEST, sridhar.samudrala at intel.com wrote: >>>> This provides a generic interface for paravirtual drivers to listen >>>> for netdev register/unregister/link

On Wed, Apr 18, 2018 at 09:13:15PM +0200, Jiri Pirko wrote: > Wed, Apr 18, 2018 at 08:43:15PM CEST, sridhar.samudrala at intel.com wrote: > >On 4/18/2018 2:25 AM, Jiri Pirko wrote: > >> Wed, Apr 11, 2018 at 09:13:52PM CEST, sridhar.samudrala at intel.com wrote: > >> > On 4/11/2018 8:51 AM, Jiri Pirko wrote: > >> > > Tue, Apr 10, 2018 at 08:59:48PM CEST,

Hi Christian, Chuck (and lists) It seems that the problem may be the strange behaviour of 'unstack' inside a function. See this thread in the R mailing list: http://tolstoy.newcastle.edu.au/R/help/04/03/1160.html Anyway, I got round the problem by using 'aggregate' instead of converting to a list and then tapply to sum values of metric. Probably more efficient as well. Thanks

Inspired by colouring the dots of box-whisker plots I am trying to also fill the boxes (rectangles) with different colours. This seems not to work as I expected. Looking at the help page of panel.bwplot it says: 'fill - color to fill the boxplot'. Obviously it is only intended to fill all boxes with only one colour? Nevertheless the following example shows, that 'fill' from

# read.table is NOT part of the data.table package #library(data.table) DFM <- read.table( text= 'obs start end 1 2/1/2015 1/1/2017 2 4/11/2010 1/1/2011 3 1/4/2006 5/3/2007 4 10/1/2007 1/1/2008 5 6/1/2011 1/1/2012 6 10/5/2004 12/1/2004 ',header = TRUE, stringsAsFactors = FALSE) # cleaner way to compute D DFM$start <- as.Date( DFM$start, format="%m/%d/%Y" ) DFM$end

Hello R users and developers, I have a general question about reshaping a wide data frame using the "reshape" command. I have a data frame consisting of 108 columns that I would like to convert to a long table and remove the metadata (embedded in the column names of the wide table) to new metadata columns in the long format. #example data frame. NOTE column names contain metadata::