Hi,
 
I am using globaltest to test the association of 3 different pData
columns with the entire gene set from 11 Rat 230_2 Affy chips.
 
My pData looks as follows:
 > pData(x.gcrma)
                 Sample.Name Tissue    State
Rat230_2_432.CEL        N7ep     ep  healthy
Rat230_2_433.CEL        N7ct     ct  healthy
Rat230_2_434.CEL       8W2ep     ep diseased
Rat230_2_436.CEL        N5ct     ct  healthy
Rat230_2_437.CEL        N5ep     ep  healthy
Rat230_2_438.CEL        N6ct     ct  healthy
Rat230_2_439.CEL        N6ep     ep  healthy
Rat230_2_440.CEL       8W5ct     ct diseased
Rat230_2_441.CEL       8W5ep     ep diseased
Rat230_2_442.CEL       8W8ct     ct diseased
Rat230_2_443.CEL       8W8ep     ep diseased

I get the following errors:
 
> gt.all <- globaltest(x.gcrma, "State", model='logistic')
Error in "[<-"(`*tmp*`, 1, , value = c(TRUE, TRUE, TRUE, TRUE, TRUE,
TRUE,  : 
        number of items to replace is not a multiple of replacement
length
> gt.all <- globaltest(x.gcrma, "Tissue")
Error in "[<-"(`*tmp*`, 1, , value = c(TRUE, TRUE, TRUE, TRUE, TRUE,
TRUE,  : 
        number of items to replace is not a multiple of replacement
length

Any ideas?
 
Thanks
 
Tarun Nayar
Genome Sciences Centre, Vancouver, BC

	[[alternative HTML version deleted]]

