Question: pData function not found. substitute? "An end to end workflow for differential gene expression using Affymetrix microarrays"
1
gravatar for RNAseqer
8 days ago by
RNAseqer 100
RNAseqer 100 wrote:

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?

pdata sva pvca bioconductor R • 81 views
ADD COMMENTlink modified 8 days ago by zx87547.3k • written 8 days ago by RNAseqer 100
3
gravatar for zx8754
8 days ago by
zx87547.3k
London
zx87547.3k wrote:

From the linked manual:

The pData function of the Biobase package...

Do you have Biobase package installed/loaded?

ADD COMMENTlink written 8 days ago by zx87547.3k

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!

ADD REPLYlink written 8 days ago by RNAseqer 100
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1654 users visited in the last hour