search for: data11

Hello, I'm hoping someone with a wide experience with R may be able to see what the program is trying to tell me. I've got an array: y1=rnorm(41,0.2) y2=rnorm(41,0.2) y3=rbind(y1,y2) data11<-array(0,c(41,2,2)) data11[,1,]=y3 data11[,2,]=y3 rownames(data11)<-rownames(data11, do.NULL = FALSE, prefix = "Obs.") colnames=c("V","R") varnames=c("one","two") dimnames(data11)<-list(rownames(data11), varnames, colnames) data11a<-a...

Hello, I'm working with an array; I'm trying to make it so that an array of dim(42,2,2) has names whose length corresponds to that of the array, and am hoping someone with experience with this can see what I'm not doing correctly: data11 = array(0,c(41,2,2)) y = lsoda(x0,times,fhn$fn.ode,pars)#This is make.fhn() from colloc infer package# y = y[,2:3] data11<-array(0,c(41,2,2)) m8<-cbind(.29,1:41) m9<-as.array(m8,dim=c(41,2)) data11[,1,] = y+m9 data11[,2,] = y+m9 rownames(data2, do.NULL = FALSE, prefix=&qu...

Suppose I have following dataset : > head(data1) Date Return 1 03/31/00 0.14230650 2 04/28/00 -0.03276228 3 05/31/00 -0.06527890 4 06/30/00 -0.04999873 5 07/31/00 -0.01447902 6 08/31/00 0.22265729 Now I convert it to zoo object : > data11 = zoo(data1[,2], as.Date(data1[,1], format="%m/%d/%y")) > head(data11) 2000-03-31 2000-04-28 2000-05-31 2000-06-30 2000-07-31 2000-08-31 0.14230650 -0.03276228 -0.06527890 -0.04999873 -0.01447902 0.22265729 Clearly those are monthly data. Therefore I want to convert i...

...lustr-03:/mnt/data06 clustr-01:/mnt/data07 clustr-02:/mnt/data07 clustr-03:/mnt/data07 clustr-01:/mnt/data08 clustr-02:/mnt/data08 clustr-03:/mnt/data08 clustr-01:/mnt/data09 clustr-02:/mnt/data09 clustr-03:/mnt/data09 clustr-01:/mnt/data10 clustr-02:/mnt/data10 clustr-03:/mnt/data10 clustr-01:/mnt/data11 clustr-02:/mnt/data11 clustr-03:/mnt/data11 clustr-01:/mnt/data12 clustr-02:/mnt/data12 clustr-03:/mnt/data12 clustr-01:/mnt/data13 clustr-02:/mnt/data13 clustr-03:/mnt/data13 clustr-01:/mnt/data14 clustr-02:/mnt/data14 clustr-03:/mnt/data14 clustr-01:/mnt/data15 clustr-02:/mnt/data15 clustr-03:/mn...

...file(file, ifelse(append, "a", "w")) : cannot open the connection In addition: Warning message: In file(file, ifelse(append, "a", "w")) : cannot open file 'c:/data1.csv': Permission denied where dataa is a zoo object > head(dataa) data11 data22 Nov 1980 988.25 194841 Dec 1980 942.38 205732 Jan 1981 935.90 226501 Feb 1981 968.79 227402 Mar 1981 932.77 233490 Apr 1981 906.18 233447 Can please tell me why this error is coming? I am using Windows vista [[alternative HTML version deleted]]

...autiful one. However, I want to change 5-year survival probability to 5-year failure probability. I couldn?t get hazard rate from Hazard(f1) because I used cph for the model. Here is my code: library(rms) f1 <- cph(Surv(retime,dfs) ~ age+her2+t_stage+n_stage+er+grade+cytcyt+Cyt_PCDK2 , data=data11, surv=T, x=T, y=T, time.inc=5) surv<- Survival(f1) haz<- Hazard(f1) Here is the Error in UseMethod("Hazard") : no applicable method for 'Hazard' applied to an object of class "c('cph', 'Design', 'coxph')" surv10 <- function(lp) surv(...

...d.delim("t (7).txt", quote="", as.is=TRUE) data8 <-read.delim("t (8).txt", quote="", as.is=TRUE) data9 <-read.delim("t (9).txt", quote="", as.is=TRUE) data10 <-read.delim("t (10).txt", quote="", as.is=TRUE) data11 <-read.delim("t (11).txt", quote="", as.is=TRUE) data12 <-read.delim("t (12).txt", quote="", as.is=TRUE) data13 <-read.delim("t (13).txt", quote="", as.is=TRUE) data14 <-read.delim("t (14).txt", quote="",...

...flt" "data10.flt" "data100.flt" "data101.flt" [5] "data102.flt" "data103.flt" "data104.flt" "data105.flt" [9] "data106.flt" "data107.flt" "data108.flt" "data109.flt" [13] "data11.flt" "data110.flt" "data111.flt" "data112.flt" [17] "data113.flt" "data114.flt" "data115.flt" "data116.flt" . . to . . [357] "data91.flt" "data92.flt" "data93.flt" "data94.flt" [3...

...x12 = rep(x1,nummod); #Put explanatory into a single vector data1 = data.frame(resp1, x12) #Data frame for kmeans cluster1 = kmeans(data1, 2, nstart=25)$cluster #Obtain cluster labels data1 = data.frame(data1 ,cluster1)#Cluster labels in third column data10 = subset(data1, cluster1==1) data11 = subset(data1, cluster1==2) model10 = lm(resp1 ~ x12, data10)#It works using the subset data frame model1 = lm(resp1 ~ x12, cluster1 == 1, data1); #Gives the following error Error in eval(expr, envir, enclos) : invalid 'envir' argument

...l data set* used for building the nomogram. *My question is how I get the formula for calculating the survival probability for this nomogram. Then I can use this formula to do validation by using other data set. * f1 <- cph(Surv(retime,dfs) ~ age+her2+t_stage+n_stage+er+cytcyt+Cyt_PCDK2 , data=data11, surv=T, x=T, y=T, time.inc=5) surv<- Survival(f1) surv10 <- function(lp) surv(10,lp) surv5 <- function(lp) surv(5,lp) quant <- Quantile(f1) at.surv <- c(0.1, 0.3, 0.5, 0.7, 0.9) nom<- nomogram(f1, conf.int=F, fun=list(surv5, surv10), funlabel=c('5-Year Survival Probabi...

...h_gen] posix7: Access to /mnt/data07//.. (on dev 16771) is crossing device (16705) [2010-05-04 10:50:30] W [posix.c:246:posix_lstat_with_gen] posix9: Access to /mnt/data09//.. (on dev 16771) is crossing device (2161) [2010-05-04 10:50:30] W [posix.c:246:posix_lstat_with_gen] posix11: Access to /mnt/data11//.. (on dev 16771) is crossing device (2113) [2010-05-04 10:50:30] W [posix.c:246:posix_lstat_with_gen] posix11: Access to /mnt/data11//.. (on dev 16771) is crossing device (2113) [2010-05-04 10:50:30] W [posix.c:246:posix_lstat_with_gen] posix13: Access to /mnt/data13//.. (on dev 16771) is crossin...

...o make a beautiful one. However, I want to change 5-year survival probability to 5-year failure probability. I couldn?t get hazard rate from Hazard(f1) because I used cph for the model. Here is my code: f1 <- cph(Surv(retime,dfs) ~ age+her2+t_stage+n_stage+er+grade+cytcyt+Cyt_PCDK2 , data=data11, surv=T, x=T, y=T, time.inc=5) surv<- Survival(f1) haz<- Hazard(f1) Here is the Error in UseMethod("Hazard") : no applicable method for 'Hazard' applied to an object of class "c('cph', 'Design', 'coxph')" surv10 <- function(lp) surv...