Twin Workshop 2010
Posted on
mkeller
Joined: 08/04/2009
Forums
Hi all,
Please post comments & responses here regarding the Twin Workshop 2010 (schedule, ideas on what to teach, ideas on how to teach, random thoughts, etc.).
I'll start it off with an
Given that the workshop participants will now be working in R, I wonder if there should be a session on simulating data from the ground up (e.g., not using a particular program/script, but rather how to simulate twin/family data using base R). Simulating data in R (which is quite easy) and then running it in OpenMx is good not only for checking that models are working as expected, but also for helping to reinforce conceptual understanding of these models. Thoughts?
Log in or register to post comments
In reply to I'll start it off with an by mkeller
Simulation is a good
the power & sample size session on Thursday afternoon. However, it perhaps needs more attention than that because it is such a fundamental technique for establishing the statistical validity of anything. Mainly because we're not as smart as Fisher, so showing things from first principles is less of an option :). Putting a smattering of mvrnorm() into the Tuesday morning R session might be a good place to start.
Log in or register to post comments
In reply to Simulation is a good by neale
Hi Neale, I am a
I am a computational biologist very new to Twin analyses. I wanted to get some suggestions from you regarding the Twin analysis. In my data set I have 24 MZ twins with a mix of some same sex male twins and some female and their age ranges from 21- 90 yrs old. Similarly, I have 10 DZ pairs with similar combination of male and female twin pairs and age range. How do I adjust the age and gender effect which might have effect on the variance estimate.
Any input would be helpful
Thanks,
Sanchita
Log in or register to post comments
The main web page for the
The schedule is at http://ibg.colorado.edu/dokuwiki/doku.php?id=workshop:2010:schedule
Most workshop faculty will have logins for the wiki, but if not (or if you don't remember your password) go to https://ibg.colorado.edu/dokuwiki/doku.php?id=workshop:2010:schedule&do=login§ok=85d53fd732bdf4492fbd3fd6c1dd3d5b and register or request a new password.
Log in or register to post comments
In reply to The main web page for the by jlessem
I wonder if this particular
<5 mins later>
Ok... hardly any sooner noted than done!
Log in or register to post comments
I wanted to check the
Are we suggesting people save the R-workspace or are you going to introduce some of the "save to text" functions?
I'm particularly thinking of the situation in which R crashes or your unix job hits a wall time before the script completes.
With Mx1.0 output wrote at checkpoints so scripts that crashed always left a 'paper trail' but I haven't been able to get this running smoothly under openMx
Log in or register to post comments
In reply to I wanted to check the by smedland
Checkpoints (and restoring
Log in or register to post comments
In reply to Checkpoints (and restoring by tbrick
checkpoint is implemented in mxRun now
mxRun(..., checkpoint = TRUE)
There are numerous other options in mxRun (suppressWarnings, useSocket etc), see ?mxRun
You can set things like checkpoints and parallel execution on multiple cores as global options like this:
# set up number of cores
mxOption(NULL, "Number of Threads", omxDetectCores() - 1)
mxOption(NULL, "Always Checkpoint", "Yes")
mxOption(NULL, "Checkpoint Count" , 1)
mxOption(NULL, "Checkpoint Units" , "evaluations")
Log in or register to post comments