similar to: Package "boot": How to construct CI from censboot object?

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*

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

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

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

Dear all, > data <- data.frame(id=seq(1:10),x=runif(10)) > data id x 1 1 0.3604464 2 2 0.4813987 3 3 0.0160058 4 4 0.7165909 5 5 0.6092248 6 6 0.2413049 7 7 0.7981568 8 8 0.6093960 9 9 0.2887064 10 10 0.3485780 > selected.id <- sample(data$id,3,replace=F) > selected.id [1] 9 7 1 I want to select data with corresponding selected.id, namely >

Thanks to all who response my problem. Here are my summary : 1. from Dirk Eddelbuettel <edd at debian.org> We could try a mixture of normals -- ie flip a coin (use a uniform with some cutoff c where 0 < c < 1 ) to choose between N(0, sigma_1) and N(0, sigma_2). 2. from Michaell Taylor <michaell.taylor at reis.com> We could use the gld library to specify the lambdas of

Thanks to all who response my problem. Here are my summary : 1. from Dirk Eddelbuettel <edd at debian.org> We could try a mixture of normals -- ie flip a coin (use a uniform with some cutoff c where 0 < c < 1 ) to choose between N(0, sigma_1) and N(0, sigma_2). 2. from Michaell Taylor <michaell.taylor at reis.com> We could use the gld library to specify the lambdas of

The latest version of the survival package has two important additions. In prior code the call coxph(Surv(time, status) ~ age + strata(inst), data=lung) could fail if a version of either Surv() or strata() existed elsewhere on the search path; the wrong function could be picked up. Second, a model with survival::strata(inst) in the formula would not do what users expect. These

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

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

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

"CALL FOR PAPERS" THE ISLAMIC SOCIETY IN SOUTH EAST ASIA CONFERENCE ON STATISTICAL AND MATHEMATICAL SCIENCES BANDUNG, 25 - 26 APRIL 2003 Fakultas Matematika dan Ilmu Pengetahuan Alam Universitas Islam Bandung Jalan Tamansari No. 1, Bandung 40116 Telepon +62 22 4203368 Ext. 136 or 135 Facsimile +62 22 4263895 E-mail: - fmipa_unisba at yahoo.com - wanda at

Dear R user, I am trying to convert the contents of a date.frame to a matrix. Since there are negative values in the date.frame, when I use data.matrix(x, rownames.force = NA), the resulting matrix is not the same as the original one. Basically I think R treats the numbers in the date.frame as character and converts it to corresponding numerics. Any idea on this issue? Many Thanks, Hongyuan

Dear all, Given mypi mypi <- c(0.1,0.2,0.2,0.1,0.3,0.4,0.4,0.4,0.4,0.2) I want to create myfreq as follows mypi myfreq 0.1 2 0.2 3 0.2 3 0.1 2 0.3 1 0.4 4 0.4 4 0.4 4 0.4 4 0.2 3 where myfreq is frequency of its corresponding observation. How to do that? Thank you, Regards, A. Kudus Institute for Math Research Univ Putra Malaysia [[alternative HTML version deleted]]

Hi, Please post any comments or questions if you get around to testing them. Your comments are very welcome! Best regards, -Gonglei > > >>> On 22.10.13 at 06:08, "Gonglei (Arei)" <arei.gonglei@huawei.com> > wrote: > > > Hi, guys. The new patch has been modified based on the principles you > > > suggested, thank you so much. > > >

Dear All, since years I am struggling with Surv objects in data.frames. The following seems to have to do with it. See below the modified example from the help page of survSplit. The original works, as expected. If, however, a Surv object is added to the data.frame, each record gets doubled. Is there some solution other than avoiding Surv objects in data.frames? Thanks, Heinz

Hello: Using the built-in dataset aml as an example: data(aml) If I use instead dummy variables: aml$x1 = (aml$x=="maintained")aml$x2 = (aml$x=="unmaintained") and I want to plot the survival curve using x1, x2, and I just want the 2 levels, rather than 4 curves from: fit <- survfit(Surv(time, status) ~ x1+x2, data=aml) plot(fit) I guess because there are 2 levels

Dear all, Right now, I'm doing research about outlier in statistical data (univariate and multivariate) and I want to simulate its behavior. My problem is : How to generate random data from distribution with zero skewness and some kurtosis values in R ? A. Kudus ===================== Dept. of Statistics Bandung Islamic University I n d o n e s i a ==========================