On Sat, Mar 5, 2011 at 4:08 PM, Stat Consult <statconsul...@gmail.com>wrote:
> Dear ALL > > I cannot run this line > > stat.obs <- apply (GS, 2, function(z) Hott2(t(DATA[which(z==1),]), cl)) > > Error in Hott2 (t(DATA[which(z == 1), ]), cl) : subscript out of bounds > > I will be glade if you guide me. > > > ******************************************************************************* > > *GS is a matrix 1857*200 > > *DATA is a matrix 1857*79 > > *cl <- as.factor(y) > > > *y<-*c(1,0,1,0,1,0,0,0,0,0,0,0,1,1,0,0,0,1,0,1,0,0,0,0,0,0,0,1,1,1,1,1,1,0,0,0,0,1,1,0,0,1,1,0,0,1,1,0,0,0,1,0,1,0*,0,1,1,1,0,1,1,1,1,0,1,0,0,0,1,0,1,1,1,1,0,1,0,1,1) > > *Hott2 <- function(x, y, var.equal=TRUE) > > > ******************************************************************************** > > Best Regards, > [[alternative HTML version deleted]] ______________________________________________ R-help@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.