Hi,

We have 4 cases and 5 controls. We would like to add RIN number as one of the covariates to the data analysis. Here is my code which gives errors:

> countsTable=read.table("4Cases_5Controls_counts.txt",header=TRUE)
> cdsFull<- newCountDataSet( countsTable, nisc1Design )
> cdsFull<- estimateSizeFactors( cdsFull )
> cdsFull<- estimateDispersions(cdsFull,"pooled")
> Error in rowSums(sapply(tapply((1:ncol(counts))[replicated_sample], factor(conditions[replicated_sample]),  :
>    'x' must be an array of at least two dimensions
>

Any suggestions will be greatly appreciated.

Thank you very much.

Regards,
Nirmala


------------------------------------------------------------------------------------
Nirmala Akula, MS, PhD
Contractor,
Human Genetics Branch
NIMH/NIH
Blg 35, Rm 1A/205
Bethesda, MD - 20892
Phone# 301-451-4258


	[[alternative HTML version deleted]]

