Account Options

  1. Sign in
The old Google Groups will be going away soon, but your browser is incompatible with the new version.
Google Groups Home
« Groups Home
mapthis.clean in genetic map construction tutorial
There are currently too many topics in this group that display first. To make this topic appear first, remove this option from another topic.
There was an error processing your request. Please try again.
flag
  2 messages - Collapse all  -  Translate all to Translated (View all originals)
The group you are posting to is a Usenet group. Messages posted to this group will make your email address visible to anyone on the Internet.
Your reply message has not been sent.
Your post was successful
 
From:
To:
Cc:
Followup To:
Add Cc | Add Followup-to | Edit Subject
Subject:
Validation:
For verification purposes please type the characters you see in the picture below or the numbers you hear by clicking the accessibility icon. Listen and type the numbers you hear
 
James E Specht  
View profile  
 More options Feb 21 2012, 6:54 pm
From: James E Specht <jespe...@unlnotes.unl.edu>
Date: Tue, 21 Feb 2012 17:54:33 -0600
Local: Tues, Feb 21 2012 6:54 pm
Subject: mapthis.clean in genetic map construction tutorial

Karl:  My students have been working their way thru your nicely documented
tutorial for Genetic map construction with R/qtl, using your supplied
command code for this tutorial.  However, one student who followed up on
the mapthis.clean set of cmds (see below) with his own extra follow-up
plot.geno command indicated to me that the eight genotypes did not seem to
have been cleaned up (i.e., converted to NA or missing) in the
mapthis.clean object.  I did note to the student that the solid or open
symbols red-boxed in the 1st plot did disappear in the 2nd plot (in some
cases replaced by an x symbol), but red boxes were still present in the
2nd plot at the same positions.  Can you offer us some insight or
clarification on this?  Thanks!  - Jim

> ### chunk number 122: plotgeno eval=FALSE
> ###################################################
> ## #line 1920 "geneticmaps.Rnw"
> ## plot.geno(mapthis, chr=1, ind=toperr$id[toperr$chr==1],
> ##           cutoff=6, include.xo=FALSE)

> ###################################################
> ### chunk number 123: plotgenoplot
> ###################################################
> #line 1927 "geneticmaps.Rnw"
> par(mar=c(4.1,4.1,0.6,0.6), las=1, cex.axis=0.9)
> plot.geno(mapthis, chr=1, ind=toperr$id[toperr$chr==1], main="",

cex=0.8,
+           include.xo=FALSE, cutoff=6)

> ###################################################
> ### chunk number 124: dropgenotypes
> ###################################################
> #line 1953 "geneticmaps.Rnw"
> mapthis.clean <- mapthis
> for(i in 1:nrow(toperr)) {

+   chr <- toperr$chr[i]
+   id <- toperr$id[i]
+   mar <- toperr$marker[i]
+   mapthis.clean$geno[[chr]]$data[mapthis$pheno$id==id, mar] <- NA
+ }

> plot.geno(mapthis.clean, chr=1,

ind=c("id200","id36","id236","id217","id235","id261","id87","id72"),cutoff= 6,include.xo=FALSE)

 
You must Sign in before you can post messages.
To post a message you must first join this group.
Please update your nickname on the subscription settings page before posting.
You do not have the permission required to post.
Karl Broman  
View profile  
 More options Feb 22 2012, 9:34 am
From: Karl Broman <kbro...@gmail.com>
Date: Wed, 22 Feb 2012 08:34:31 -0600
Local: Wed, Feb 22 2012 9:34 am
Subject: Re: [Rqtl-disc] mapthis.clean in genetic map construction tutorial
You need to first re-run calc.errorlod.  plot.geno(mapthis.clean, ...) is looking at the previously calculated error LOD scores.

> mapthis.clean <- calc.errorlod(mapthis.clean, error.prob=0.005)
> top.errorlod(mapthis.clean, cutoff=6)

        No errorlods above cutoff.

karl

On Feb 21, 2012, at 5:54 PM, James E Specht wrote:


 
You must Sign in before you can post messages.
To post a message you must first join this group.
Please update your nickname on the subscription settings page before posting.
You do not have the permission required to post.
End of messages
« Back to Discussions « Newer topic     Older topic »