Question: eqtl using ggtools
0
gravatar for archie
11 months ago by
archie70
India
archie70 wrote:

Hi everyone

I am trying to run GGtools for prediction of eqtl based on the genotypes and expression information. https://bioconductor.org/packages/release/bioc/vignettes/GGtools/inst/doc/GGtools.pdf

I am trying to run function gwSnpTests which needs the exxpression object and phontype information

step 1

es = new("ExpressionSet", exprs=exprs1,phenoData=phenoData,featureData = fdata)

Expressionset variable (es)

ExpressionSet (storageMode: lockedEnvironment) assayData: 55 features, 19 samples element names: exprs protocolData: none phenoData sampleNames: F545 F546 ... F563 (19 total) varLabels: ID S Ph varMetadata: labelDescription featureData featureNames: AC090958.3 CR847794.2 ... MIR4774 (55 total) fvarLabels: GeneID fvarMetadata: labelDescription experimentData: use 'experimentData(object)' Annotation:

step 2

ss = make_smlSet(es, list("1" =xx$genotypes) ,annotation(es))

smlset variable (ss)

SnpMatrix-based genotype set: number of samples: 19 number of chromosomes present: 1 annotation:
Expression data dims: 14000 x 19 Total number of SNP: 4500 Phenodata: An object of class 'AnnotatedDataFrame' sampleNames: F545 F546 ... F563 (19 total) varLabels: ID S Ph varMetadata: labelDescription

step 3.

I am using the fvarLabels: genesym, also present in my expressionset fdata.

t1 = gwSnpTests(genesym("sym")~S , ss)

character(0) Failed with error: ‘'package' must be of length 1’ Error in revmap(get(paste(gsub(".db", "", annpack), "SYMBOL", sep = ""))) : error in evaluating the argument 'x' in selecting a method for function 'revmap': Error in get(paste(gsub(".db", "", annpack), "SYMBOL", sep = "")) : object 'SYMBOL' not found

How to define and create the instances genesym to run function gwSnpTests

I am stuck at this step. Please help me to sort this issue.

I will appreciate all the suggestion.

ggtools • 263 views
ADD COMMENTlink modified 11 months ago • written 11 months ago by archie70
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: 1846 users visited in the last hour