Entering edit mode
SCAN.UPC was used to normalize Hu-133A arrays, and the run worked fine
using the syntax:
normalized=SCAN(celFilePath, "output_file.txt")
Next, I would like to implement ComBat in SVA using the syntax:
library("sva")
exprdata<- as.matrix(read.table("output_file.txt", header=T))
batch<- as.matrix(read.table("batch.txt", header=T))
ComBat(exprdata,batch)
(FYI, "exprdata" exists in memory as a 22283x919 double matrix, and
"batch" is 919 x 1 character matrix with the GEO dataset name, e.g.,
"GSE1456")
However, an exception is being thrown in ComBat which reads as
follows:
> ComBat(data,batch)
Reading Sample Information File
Error in file(file, "rt") : invalid 'description' argument
Any suggestions?
LP
Houston Methodist. Leading Medicine.
Ranked by U.S.News & World Report as one of America's "Best Hospitals"
in 13 specialties. Named to FORTUNE? Magazine's "100 Best Companies to
Work For?" list eight years in a row. Designated as a Magnet hospital
for excellence in nursing. Visit us at houstonmethodist.org. Follow us
at twitter.com/MethodistHosp and www.facebook.com/HoustonMethodist
***CONFIDENTIALITY NOTICE*** This e-mail is the property of Houston
Methodist Hospital and/or its relevant affiliates and may contain
restricted and privileged material for the sole use of the intended
recipient(s). Any review, use, distribution or disclosure by others is
strictly prohibited. If you are not the intended recipient (or
authorized to receive for the recipient), please contact the sender
and delete all copies of the message. Thank you.
Hi,
I see that you got it working.. Can you tell me how you overcame the 1st error - "Error in file(file, "rt") : invalid 'description' argument"
I am having the same problem. Thanks.