Question: eqtl using ggtools
0
gravatar for archie
19 months ago by
archie90
India
archie90 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 • 403 views
ADD COMMENTlink modified 19 months ago • written 19 months ago by archie90
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: 1769 users visited in the last hour