Displaying 3 results from an estimated 3 matches for "polytim".
Did you mean:
polytime
2012 Apr 25
2
Accessing a list
...<- c (c(1:12), c(1:12), c(1:12), c(1:12), c(1:12), c(1:12), c(1:12),
c(1:6) )
length(ts)
data.use = cbind(ts, data.use)
#a.out2 <- amelia(data.use, m = 1000, idvars = "year")
n.times = 100
a.out.time <- amelia(data.use, m = n.times, ts="ts", idvars="year",
polytime=2)
constant.col = dischargex1.col = pressurex2.col = windx3.col =
rep(0,n.times)
for (i in 1: n.times)
{
x = c("imp",i)
df = paste(x, collapse = "")
data1 = a.out.time[[1]]$df
attach(data1)
y = as.numeric(y)
dischargex1 = as.numeric(dischargex1)
pressurex2 = as.numeric(pressu...
2013 Feb 15
0
Ho w Do I Get Cox Model Convergence After Multiple Imputation
...rcentage of
missing values
myimp <- amelia(impdat, m=M, p2s=0,
noms=c("postrantb","vlsupp","whostage","resistance"),cs=c("id"),ts=c("X_t0"),bounds=matrix(c(3,0,70,
4,0,100, 11,0,400, 12,0,50
),ncol=3,nrow=4,byrow=T),logs="nvp",polytime=3,splinetime=3,empri=1000,incheck=TRUE,tolerance=0.001)
The I do cox model as follows
# Categorization of relevant variables
for(i in 1:10){
myimp$imputations[[i]] <-
cbind(myimp$imputations[[i]],cut(myimp$imputations[[i]]$t0rvl,breaks=c(-1,50,400,1000000000)))
colnames(myimp$imputations[[i]])[...
2013 Jan 28
6
Thank you your help.
Hi,
temp3<- read.table(text="
ID CTIME WEIGHT
HM001 1223 24.0
HM001 1224 25.2
HM001 1225 23.1
HM001 1226 NA
HM001 1227 32.1
HM001 1228 32.4
HM001 1229 1323.2
HM001 1230 27.4
HM001 1231 22.4236 #changed here to test the previous solution
",sep="",header=TRUE,stringsAsFactors=FALSE)
?tempnew<- na.omit(temp3)
?grep("\\d{4}",temp3$WEIGHT)
#[1] 7 9 #not correct