Question: eqtl using ggtools
0
gravatar for archie
7 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 • 178 views
ADD COMMENTlink modified 7 months ago • written 7 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: 1814 users visited in the last hour