I am currently tinkering with an analysis bioconductor ppackage found here:
"An end to end workflow for differential gene expression using Affymetrix microarrays" https://bioconductor.org/packages/devel/workflows/vignettes/maEndToEnd/inst/doc/MA-Workflow.html
Under section 4.5 you perform the task of running a general linearized model where sva surrogate variables are modeled as a function of your covariates. The following does this by pulling from an expression set object created in the first part of the pipeline:
glm.sv1 <-glm(pData(inpData_sv)[,"sv1"]~pData(inpData_sv)[,"Batch"]+pData(inpData_sv)[,"Sex"])
Unfortunately, I am getting an error telling me:
Error in pData(inpData_sv) : could not find function "pData"
I have tried getting the pData function by downloading ballgown and phylobase, the only two packages where I have seen pData(). Neither of these have seemed to provide me with the pData function required here. I'm pretty new to R, so could anyone recommend a substitute function for this R command?
I sure thought I had biobase. I Re-installed it just in case.
Loading the library though...that I forgot.
This seems to have fixed the problem. Thanks!