similar to: (Censboot, Z-score, Cox) How to use Z-score as the statistic within censboot?

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

Dear R users, Below I have written 4 functions CIT1, CIT2a and CIT2b and CIT3 which recode a variable CLD_ISCH into 3 new variables(T1 T2 T3), I wish to use T1, T2 and T3 based on the values of tf1 and tf2. (NOTE:- T2a is used to create T2 in a long winded manner due to my lack of programming experience) I then attach T1 T2 and T3 to a dataset KidneyT that contains other variables i wish to use

Dear all, I run the example of "censboot" contained in "boot" package. But, I can't find the confidence interval of the resulted "censboot" object. Any idea ? > aml.fun <- function(data) { + surv <- survfit(Surv(time, cens)~group, data=data) + out <- NULL + st <- 1 + for (s in 1:length(surv$strata)) { +

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

Hi, I tried to do a nested Anova with the attached Data. My response variable is "survivors" and I would like to know the effect of (insect-egg clutch) "size", "position" (of clutch on twig) and "clone" (/plant genotype) on the survival of eggs (due to predation). Each plant was provided with three different sizes of clutches (45,15,5) and had

Dear R users, I am trying to convert a dataset from one format to several rectangular datasets. A consultant helped design the data entry program for our survey using Delphi/Pascal and for each household the information is stored in a file called "EA-HM-HH.TXT" where EA is the enumeration area number, HM is the homestead number and HH is the household number. Within this file the

Hi I was looking through the help on how to get standardized scores and came across an anomaly. The help.search("zscore") reports there is such a function in R.base. However, when I try to get help on zscore it reports no such function and says that R.base in built for win32. I'm running the Mac OS X version of R with the smorgasboard of packages that are included in that

Hi, I am writing a basic function to extract the z scores for some linear regression coefficients: zscore<-function( y, x) { lm<-lm( y ~ x ) z <- coef(lm)/sqrt(diag(vcov(lm))) return(z) } I would like to pass a dataframe to the function as a argument so the function call changes from zscore(df$y1,df$x1) to zscore(y1,x1,data=df) but I am not sure how to reference the data

Hi, R 1.0.1, on Linux(mandrake/RedHat). I am trying to write my first function using R and I need help figuring out the error I have indicated in the Subject line above. The purpose of the function is not to "wow" anyone, but just to serve as a stepping stone so that I might write more complicated functions in the future. First I created my function named "CIfunc" and saved

Dear mailing list, I tried to run the example for the conditional bootstap written in the help file of censboot. I got the following result: STRATIFIED CONDITIONAL BOOTSTRAP FOR CENSORED DATA Call: censboot(data = aml, statistic = aml.fun, R = 499, F.surv = aml.s1, G.surv = aml.s2, strata = aml$group, sim = "cond") Bootstrap Statistics : original bias std. error t1*

Forwarded Message: > To: Steve Arthur <sarthur at crick.protogene.com> > From: Thomas Lumley <thomas at biostat.washington.edu> > Subject: Re: [R] couldn't find function "..." > Date: Wed, 5 Jul 2000 16:16:25 -0700 (PDT) > ----- > On Wed, 5 Jul 2000, Steve Arthur wrote: > > > Hi, > > > > R 1.0.1, on Linux(mandrake/RedHat).

Hi All: I was looking for, but could not locate in the packages, or in the R archive searches if there exists an R implementation of Cuzick's test of trend. The test is described as follows: An extension of the Wilcoxon rank-sum test is developed to handle the situation in which a variable is measured for individuals in three or more (ordered) groups and a non-parametric test for trend

Hello, I have been trying to use the censored bootstrap function 'censboot' from the library boot, with the function 'survfit' from the library survival5. My code is: surv.fun <-function(data) { sur<-survfit(Surv(data[,1],data[,2]), data=data) } c.boot<-censboot(m,surv.fun,R=99) m is my data set I get the error message: Error in

Hi, The following code, from Angelo Canty article on line "Resampling Methods in R: the boot Package, 2002", works fine for Angelo Canty using R 2.6.0 on Windows XP. It also works for me using R 1.2.1 and S-PLUS 2000 on Windows XP after installing the S-PLUS bootstrap library, with slight differences in my outputs. > library(boot) > library(survival) >

library(survival) library(boot) data=NULL lambda=NULL result=NULL pat=rep(1:102,each=1) trt=rep(c(1,0),51) status=rep(1,102) site=rep(1:51, each=2) nr.datasets=100 seed=2006 beta=log(1/2) for (i in 1:51) { lambda[i]=1+((3-1)/50)*(i-1)} lambda1=rep(lambda, each=2) dummy=rep(c(exp(beta),1),51) elf=lambda1*dummy r=70 #the number of bootstrap replicates

Hi, The following code, from Angelo Canty article on line "Resampling Methods in R: the boot Package, 2002", works fine for Angelo Canty using R 2.6.0 on Windows XP. It also works for me using R 1.2.1 and S-PLUS 2000 on Windows XP after installing the S-PLUS bootstrap library, with slight differences in my outputs. > library(boot) > library(survival) >

Hi, The following code, from Angelo Canty article on line "Resampling Methods in R: the boot Package, 2002", works fine for Angelo Canty using R 2.6.0 on Windows XP. It also works for me using R 1.2.1 and S-PLUS 2000 on Windows XP after installing the S-PLUS bootstrap library, with slight differences in my outputs. > library(boot) > library(survival)

Hi, The following code, from Angelo Canty article on line "Resampling Methods in R: the boot Package, 2002", works fine for Angelo Canty using R 2.6.0 on Windows XP. It also works for me using R 1.2.1 and S-PLUS 2000 on Windows XP after installing the S-PLUS bootstrap library, with slight differences in my outputs. > library(boot) >

-----BEGIN PGP SIGNED MESSAGE----- Hash: SHA1 Hi, I just installed R (2.11.1) 64bit under Windows 7 64 bit When trying to readin files I get the following error message: > > source(file="C:\\Users\\me\\Documents\\My PhD\\Modelling\\R\\Scripts\\Functions\\functions 20jan10.r") Error in file(file, "rt") : cannot open the connection In addition: Warning message: In

Dear useRs, The seqinR package is a library of utilities to retrieve and analyse biological sequences. A new version of seqinR, seqinR 1.0-7, has been released on CRAN. Here is a summary of changes: o A new *experimental* function extractseqs() to download sequences thru zlib compressed sockets from an ACNUC server is released. Preliminary tests suggest that working with about 100,000