[R] Help with code
Petr PIKAL
petr.pikal at precheza.cz
Thu Dec 22 08:23:32 CET 2011
Hi
>
> I do not know how to use dput, i am attaching the txt file for the data
dput(any.object)
puts a structure of this object to console. You can copy it to your email
and anybody can copy it back to R.
Or you can transfer the structure to file see
?dput, ?dget
> http://r.789695.n4.nabble.com/file/n4222616/foo.txt foo.txt
>
> c1<-read.dlim('foo.txt')
> c2<-c1
>
> any_comp<-NULL
>
> for( i in 1:dim(c1)[1])
> {
> num_comp<-0
> for (j in 1:dim(c1)[2])
> if (c1[i,j]==2) num_comp=num_comp+1 #"Y"=2
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
this line is never true. There is no 2 in a c1.
Maybe you shall really tell us what is your intention with some simple
example.
Here is a subset of your file c1
my.c1 <- structure(list(HTN = structure(c(2L, 1L, 1L), .Label = c("",
"Y"), class = "factor"), HTN_FDR = structure(c(2L, 2L, 2L), .Label = c("",
"Y"), class = "factor"), Dyslipidemia = structure(c(2L, 1L, 2L
), .Label = c("", "Y"), class = "factor"), CAD = c(NA, NA, NA
), CAD_FDR = structure(c(2L, 2L, 1L), .Label = c("", "Y"), class =
"factor"),
Prior_MI = c(NA, NA, NA)), .Names = c("HTN", "HTN_FDR",
"Dyslipidemia",
"CAD", "CAD_FDR", "Prior_MI"), row.names = c("1", "5", "9"), class =
"data.frame")
Forget your complicated code and tell us what is you desired output?
Regards
Petr
> for (j in 1:dim(c1)[2])
> if(num_comp>0)
> {
> if (data$t1d_ptype[i] == "T1D" && c1[i ,j] == 2)
c2[i,j]<-"T1D_w"
> if (data$t1d_ptype[i] == "T1D" && c1[i, j] == 1)
c2[i,j]<-"T1D_oc"
> if(substr(data$t1d_ptype[i],1,4) == "Ctrl" && c1[i,j] == 2)
> c2[i,j]<-"Ctrl_w"
> if (substr(data$t1d_ptype[i],1,4) == "Ctrl" && c1[i,j] == 1)
> c2[i,j]<-"Ctrl_oc"
> }
> else
> {
> if(data$t1d_ptype[i] == "T1D") c2[i,j]<-"T1D_noc"
> if(substr(data$t1d_ptype[i],1,4) == "Ctrl")
c2[i,j]<-"Ctrl_noc"
> }
> }
>
> this gives me whether the particular t1d_ptype has a specific
complication
> as well as there is another complication as well.
>
> I will appreciate help very much.
> thanks
> sharad
>
> --
> View this message in context: http://r.789695.n4.nabble.com/Help-with-
> code-tp4218989p4222616.html
> Sent from the R help mailing list archive at Nabble.com.
>
> ______________________________________________
> R-help at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide
http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.
More information about the R-help
mailing list