search for: qst

Displaying 5 results from an estimated 5 matches for "qst".

Did you mean: dst
2002 Aug 28
0
Extracting variance component estimates from lme
...e... The short form of my main question is: how do I extract variance components from an lme object? The longer form (plus optional supplementary question!): I'm looking at some quantitative genetics, and want to estimate two variance components so that I can then calculate a statistic called Qst from them. So I have this: reg1 <- lme(y ~ temp*food, random =~1|POP/SIRE/DAM, na.action=na.exclude)) summary(reg1) which shows me the random effects as standard deviations. I then want the SIRE and POP variance components (as Vsire and Vpop respectively) to put into the calculation of Q...
2008 Mar 18
1
Updating selected variables in a data frame
...replace=TRUE), q2=sample(1:5, size=10, replace=TRUE)) # Generate a corresponding recodes dataset recodes <- reshape(main, idvar="id", direction="long", varying=list(2:3), v.names="value", timevar="qst") recodes <- recodes[sample(seq(nrow(recodes)), 10, replace=FALSE),] recodes$value <- recodes$value + sign(runif(1)-0.5) I could easily use a for() loop to update the main dataset. The recodes dataset can get quite large, however, so I would rather avoid this. The q...
2007 Jan 03
1
mcmcsamp and variance ratios
Hi folks, I have assumed that ratios of variance components (Fst and Qst in population genetics) could be estimated using the output of mcmcsamp (the series on mcmc sample estimates of variance components). What I have started to do is to use the matrix output that included the log(variances), exponentiate, calculate the relevant ratio, and apply either quantil...
2009 Oct 07
1
output validations errors form causing type error?
...must be a U.S. state selected from the drop down list for customers in United States"], "card_verifiation"=>[["is required"]]}, @base=#<Order id: nil, created_at: nil, updated_at: nil, customer_id: nil, gst: #<BigDecimal:105d342c0,''0.0'',9(18)>, qst: #<BigDecimal: 105d33ed8,''0.0'',9(18)>, shipping_ca: #<BigDecimal:105d33ac8,''0.0'',9(18) >, shipping_us: #<BigDecimal:105bfb840,''0.7E1'',9(18)>, currency: "USD", status: nil, comments: nil, customer_ip: "127.0.0.1&q...
2012 Jan 14
2
How can I doing Quality adjusted survival analysis in R?
Hi R users I need to estimate, with kaplan Meier methodology, a Quality adjusted survival analysis. It is possible doing this at R? Thanks in advance. Best Regards Pedro Mota Veiga -- View this message in context: http://r.789695.n4.nabble.com/How-can-I-doing-Quality-adjusted-survival-analysis-in-R-tp4295868p4295868.html Sent from the R help mailing list archive at Nabble.com.