Hi,
I apologize if this is the wrong place to ask this question or if it's been asked already but I am having an issue with ballgown. I'm trying to perform the analysis on data that is associated with the Pertea et al, 2016 nature protocol paper. My issue is with step 9. Here Is my code
library(ballgown) library(RSkittleBrewer) library(genefilter) library(dplyr) library(devtools) pheno_data <-read.csv("geuvadis_phenodata.csv") bg_chrX = ballgown(dataDir = "ballgown", samplePattern = "ERR", pData = pheno_data)
The error I get is
Fri Jan 25 12:20:27 2019: Reading linking tables
Fri Jan 25 12:20:27 2019: Reading intron data files
Fri Jan 25 12:20:28 2019: Merging intron data
Fri Jan 25 12:20:28 2019: Reading exon data files
Fri Jan 25 12:20:29 2019: Merging exon data
Fri Jan 25 12:20:29 2019: Reading transcript data files
Fri Jan 25 12:20:29 2019: Merging transcript data
Error in ballgown(dataDir = "ballgown", samplePattern = "ERR", pData = pheno_data) :
is.null(pData) | class(pData) == "data.frame" are not all TRUE`
However if I run all(pheno_data$ids ==list.files("ballgown"))
I get TRUE. So I'm a bit confused. Any help is appreciated.
A simple
as.dataframe(phenodata
) fixes the problem. Thanks!!