similar to: Using censummary on censored data

Hi Forum, I am new and learning this research tool. was trying to run some tests on a censored dataset using the NADA library. Everything is hunky dory except 'censummary' throws an error and I am trying to figure out a reason for this. Please help! below you will see the error statement... thnx --nk > library(NADA) Loading required package: survival Loading required package: splines

HI I am trying to analyse data which is left-censored (i.e. has values below the detection limit). I have been using the NADA package of R to derive summary statistics and do some regression. I am now trying to carry out regression on paired data where both my X and Y have left-censored data within them. I have tried various commands in R: rega = cenreg(Cen(conc, cens_ind) ~ Gp_ident))? with

Hi , I am trying to get some estimator based on lognormal distribution when we have left,interval, and right censored data. Since, there is now avalible pakage in R can help me in this, I had to write my own code using Newton Raphson method which requires first and second derivative of log likelihood but my problem after runing the code is the estimators were too high. with this email ,I provide

Hi, I want to plot non parametric estimates of the empirical hazard function for right censored data. I've tried many functions from different packages (muhaz, Design, survival, eha, event), but none of them gave me what I wanted. Am I missing something? Here's what I want. The data below is the same used by Kiefer (J. Economic Literature, 1988), which in turns use a subset of the data

Your function is giving NaN's during the optimization. The R-forge version of optimx() has functionality specifically intended to deal with this. NOTE: the CRAN version does not, and the R-forge version still has some glitches! However, I easily ran the code you supplied by changing optim to optimx in the penultimate line. Here's the final output. KKT condition testing Number of

Hello, I wish to?censor 10% of my sample units of 50 from a Weibull distribution. Below is the code for it. I will need to know whether what i have done is correct and if not, can i have any suggestion to improve it? Thank you ?p=2;b=120 n=50 r=45 t<-rweibull(r,shape=p,scale=b) meantrue<-gamma(1+(1/p))*b meantrue cen<- runif(n-r,min=0,max=meantrue) cen Chris Guure Researcher,

Hello, When i run the code below from Weibull distribution with 30% censoring by using optim i get an error form R, which states that Error in optim(start, fn = z, data = q, hessian = T) :? ? objective function in optim evaluates to length 25 not 1 can somebody?help me remove this error. Is my censoring approach correct. n=25;rr=1000 p=1.5;b=1.2 for (i in 1:rr){ q<-c(t,cen)

Thank you for Prof.Ripley's suggestion. I fixed the program by adding a lower scope, and the program ran, but I still got warning messages, and don't know what is going on, would this affect my results? ... Step: AIC= 12337.74 Surv(tlfup, cen) ~ MI[[j]]$trt + MI[[j]]$agem40 + MI[[j]]$agem40sq + mhtypeed1 + mhtypeed2 Df AIC <none> 12338 -

Hi, Found this in logfiles: Jul 25 14:12:38 dovecot: pop3(cen at so.red): Error: UIDL: File name lookup failed: Message was expunged Jul 25 14:12:38 dovecot: pop3(cen at so.red): Disconnected: POP3 UIDLs couldn't be listed top=187/356348, retr=0/0, del=0/1207, size=789346199 Jul 25 14:12:38 dovecot: pop3(cen at so.red): Panic: file ioloop.c: line 39 (io_add_file): assertion failed: (fd

Hi, there: It's my first time to post question in this forum, so thanks for your tolerance if my question is too naive. I am using a nonparametric smoothing procedure in sm package to generate smoothed survival curves for continuous covariate. I want to truncate the suvival curve and only display the part with covariate value between 0 and 7. The following is the code I wrote:

Hi, I'm analyzing a large number of large simulation datasets, and I've isolated one of the bottlenecks. Any help in speeding it up would be appreciated. `dat` is a dataframe of samples from a regular grid. The first two columns are the spatial coordinates of the samples, the remaining 20 columns are the abundances of species in each cell. I need to calculate the species richness in

Hello all, I am trying to use stepwise procedure to select covariates in Cox model and use bootstrap to repeat stepwise selection, then record how many times variables are chosen by step() in bootstrap replications. When I use step() (or stepAIC) to do model selection, I got errors. Here is the part of my code for (j in 1:mm){ #<--mm=10 for (b in 1:nrow(reg.bs)){ #<--bootstrap 10

Hi all - (this is posted to r-help and R-SIG-MAC) OSX 10.3.7, R 2.0.1, lme4/Matrix/latticeExtra latest, fresh install of R. MASS loaded (or not). I am getting an error message for the print() and summary() commands with all lme models I try and run in lme4 (GLMM's work fine). Using the example from the lme help, summary and print produce the following errors, despite the model being

Dear list, i am a student of psychology and have to do a multilevelanalysis on some data. About that i have one general and one specific question. This is what i have copied from the help-file on lme: data(bdf) fm <- lme(langPOST ~ IQ.ver.cen + avg.IQ.ver.cen, data = bdf, random = ~ IQ.ver.cen | schoolNR) summary(fm) after summary(fm) i get the following error:

I would like to simulate a stochastic process by specifying an intensity function, say lambda(t)=0.05 + 3 exp(鈭抰) I have the following code, but it has a number of problems. One issue I have is I do not know where the time variable is. I would appreciate your help with this code. library("survival") n <- 5000 lamc<-0.08 k1<-0.05 k2<-3 ld<-function(t){k1+k2}

Hi, Got a new one on 2.2.13: Sep 19 04:11:39 moxi dovecot: imap(cen at cor.ed): Error: Log synchronization error at seq=3,offset=144 for /Maildir/.Spam/dovecot.index: Extension record update for invalid uid=494 Sep 19 04:11:39 moxi dovecot: imap(cen at cor.ed): Warning: fscking index file /Maildir/.Spam/dovecot.index Sep 19 04:11:39 moxi dovecot: imap(cen at cor.ed): Panic: file

<div align="center"> <strong><font face="Trebuchet MS" color="#ff0000" size="6">Credit Card Shoppers Are Here To Stay!</font></strong></div> <div align="center"><strong><font face="Trebuchet MS" color="#000080" size="6">Make Sure They Are Buying from

Would someone be able to help with this question? I'm using the Gcmrec, Survrec, and Design packages to do a power analysis on simulated data. I'm receiving an error after using the Survr function that all data must have a censoring time even after using the gcmrec function: newdata<-addCenTime(olddata). My program is below. I'd greatly appreciate any help!

Hi, I am trying to get the maximum likelihood estimator for lognormal distribution with censored data;when we have left, interval and right censord. I built my code in R, by writing the deriving of log likelihood function and using newton raphson method but my estimators were too high " overestimation", where the values exceed the 1000 in some runing of my code. is there any one can

Hi there, I have a question on how to extract the linear term in the penalized spline in coxph. Here is a sample code: n=100 set.seed(1) x=runif(100) f1 = cos(2*pi*x) hazard = exp(f1) T = 0 for (i in 1:100) { T[i] = rexp(1,hazard[i]) } C = runif(n)*4 cen = T<=C y = T*(cen) + C*(1-cen) data.tr=cbind(y,cen,x) fit=coxph(Surv(data.tr[,1],