[R] Help with code

[Petr PIKAL]

Hi
> 
> hello gurus,
> 
> i have a data frame like this
>    HTN HTN_FDR Dyslipidemia CAD t1d_ptype[1:25]
> 1    Y       Y            Y                 T1D
> 2                                           T1D
> 3                                      Ctrl_FDR
> 4                                           T1D
> 5            Y                             Ctrl
> 6                                          Ctrl
> 7                                      Ctrl_FDR
> 8                                           T1D
> 9            Y            Y                 T1D
> 10                                          T1D
> 11                                     Ctrl_FDR
> 12           Y            Y                 T1D
> 13   Y       Y            Y                 T1D
> 14                                          T1D
> 15                                         Ctrl
> 16                                         Ctrl
> 17                                     Ctrl_FDR
> 18                                          T1D
> 19                                          T1D
> 20           Y                              T1D
> 21                                     Ctrl_FDR
> 22                                     Ctrl_FDR
> 23                                         Ctrl
> 24                                         Ctrl
> 25                                          T1D
> 
> i am converting it to define the groups more uniformly using this code:
> 
> 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
>   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"
>        }
> }
> 
> it is giving me error
> In `[<-.factor`(`*tmp*`, iseq, value = structure(c(NA,  ... :
>   invalid factor level, NAs generated
> 
> Also it there a simple way to do this.

To do what?
The only error i get is
Error: object 'c1' not found
 
You probably named your data frame c1 but I do not see any numbers in c1
With your programming style you shall probably use C+ or similar. In R it 
is almost always better to operate on whole objects.

Without knowing what you want to achive and without any data I can only 
suggest to look at

?factor
?interaction

Regards
Petr


> Thanks
> Sharad
> 
> --
> View this message in context: http://r.789695.n4.nabble.com/Help-with-
> code-tp4218989p4218989.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.