search for: t_n

I have the following xts objetct "temp" > str(temp) An ?xts? object from 2010-12-26 to 2011-03-05 containing: Data: num [1:70, 1] 2.95 0.852 -0.139 1.347 2.485 ... - attr(*, "dimnames")=List of 2 ..$ : NULL ..$ : chr "t_n" Indexed by objects of class: [POSIXct,POSIXt] TZ: GMT xts Attributes: NULL > temp t_n 2010-12-26 2.9500000 2010-12-27 0.8520833 2010-12-28 -0.1390625 ........... I would like to associate another column with the day of week in the form of 1=Mon, 2=Tue, ..., 7=...

Hello, i got an warning message in the following code: f<-1:100 t<-1:100 b<-100 ll2 <- function(b,f,t) { t<-cumsum(t) tn<-t[length(t)] i<-seq(along=f) s1<-(tn*exp(-b*tn)*sum(f[i]))/(1-exp(-b*tn)) s2<-sum((f[i]*(t[i]*exp(-b*t[i])-t[i-1]*exp(b*t[i-1])))/(exp(-b*t[i-1])-exp(-b*t[i]))) s1-s2 } ll2(b,f,t) i think, the problem here is, that t[0] doesn't

...which is not yet supposed to work, it seems that glm does not work either. Without Survival glm seems to work. I am not sure whether this gives a clue concerning the problems with Survival. One interesting application is to Weibull data which is only observed at single inspection points T_1, ..., T_n, not all the same and contained in insp.times. The only thing that is known about the actual lifetimes X_1, ..., X_n is whether X_i < T_i or not, giving me the 0,1 vector status. One can then get the mle's of the Weibull scale (alpha) and shape (beta) parameters by using dat <- data.fra...