Displaying 4 results from an estimated 4 matches for "uncen".
Did you mean:
unces
2008 Dec 28
1
Random coefficients model with a covariate: coxme function
...2,5,6,5,7,3,4,5,2,4,1,5,6,7,6,5,2,5,6,5,7,3,4,5,
6,5,7,8,5,6,8,6,5,6,8,6,5,8,6,9,6,5,7,8,5,6,8,6,5,6,8,6,5,8,6,9),2)
treat<-rep(c(0,1),each=64)
covar<-
rep(c(1,3,2,5,6,4,7,8,9,6,4,3,4,3,2,1,4,3,5,4,6,3,4,5,6,4,6,3,6,7,4,5,3,4,5,6,4,6,6,4,7,6,5,7,5,4,3,2,8,6,4,9,7,9,5,4,6,8,6,4,6,2,5,8),2)
uncens<-
rep(c(0,1,0,0,0,1,0,0,1,1,1,1,0,1,1,0,0,1,0,0,0,1,0,0,1,1,1,1,0,1,1,0
,0,1,1,0,0,1,0,0,1,0,1,0,0,0,1,1,0,1,1,0,0,1,0,0,1,0,1,0,0,0,1,1),2)
centers<- rep(rep(1:32, 2),2)
data1<-list(y, treat, covar, uncens, centers)
names(data1)<- c("y", "treat", "covar"...
2013 Mar 21
0
"[[i]]$" <- "" indexing and lapply
...mes of the
columns rather than the site name.
## modelEstimation -- this is the code that I use to do the loop through
modelEstimation
for(i in Daily) {
modelEstimation(localDaily = Daily[[i]], localSample = Sample[[i]],
localINFO = INFO[[i]], windowY=10, windowQ=2,
windowS=0.5,minNumObs=10,minNumUncen=6, env=parent.frame())
}
# shows the correct naming convention before being ran through
modelEstimation
> dput(Daily)
structure(list(`21NC02WQ.C1000000` = structure(list(Date =
structure(c(11231,
11232, 11233, 11234, 11235, 11236), class = "Date"), Q =
c(0.886534950681092,
0.8579370...
2012 Aug 31
3
fitting lognormal censored data
...#X2 Covariate=Type of treatment(1=chemo,0=radio) #
L1<-length(cens) #number of censored#
for (j in 1:L1)
{
if ((cens[j]==1)&(curd[j]==0)) {(cen[j]=1)&(cur[j]=1)}
else {(cen[j]=cens[j])&(cur[j]=curd[j])}
}
Now, the following is my data:
####### My Data only with uncensored and right censored ####################
data=data.frame(Ti=dat1$time,Censored=cen)
#################### Estimation using Surv pakage ############################
library(survival)
survreg(Surv(Ti, Censored)~1, data=data, dist="lognormal")
########### Seperating the data for simply...
2012 Aug 29
2
Estimation parameters of lognormal censored data
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