Question: eqtl using ggtools
gravatar for archie
2.4 years ago by
archie100 wrote:

Hi everyone

I am trying to run GGtools for prediction of eqtl based on the genotypes and expression information.

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 • 561 views
ADD COMMENTlink modified 2.4 years ago • written 2.4 years ago by archie100
Please log in to add an answer.


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