Stack Overflow
Stack Overflow

Questions tagged [asreml]

10 questions
4
votes
1 answer

Mean variance of a difference of BLUEs or BLUPs in `lme4`

Given below is the code for analysis of a resolvable alpha design (alpha lattice design) using the R package asreml. # load the data library(agridat) data(john.alpha) dat <- john.alpha # load asreml library(asreml) # model1 - random…
Crops
  • 5,024
  • 5
  • 38
  • 65
2
votes
0 answers

How did using a pedigree influence the degrees of freedom available to test fixed effects in asreml-R?

I am using a half-sibling pedigree (individuals, Sires, Dams and groups (populations) assigned). When I run the model I received the message “Terms with zero df listed …” and when I ask for the wald test of fixed effect one of my effects is missing.…
Rebecca Moss
  • 21
  • 1
2
votes
1 answer

R with asreml package: Setting a keyword for multivariate model, however it just stays keyword is an "object not found"

I am using a dataset from an online practice tutorial and the code can be found at the bottom of Page 7 (https://tomhouslay.files.wordpress.com/2017/02/indivvar_mv_tutorial_asreml.pdf) In the tutorial, they state that they are introducing "trait" as…
qsxsqsxs
  • 161
  • 1
  • 8
2
votes
1 answer

Standard error for random effect (BLUP) of a mixed effect model in r

Can we calculate standard errors (SE) for BLUPs from a mixed model in R? I am using a custom package called ASReml and it can calculate SE from the predicted value. However I am not sure where to find the SE for the BLUP values. So far I have…
FrosyFeet456
  • 349
  • 2
  • 12
1
vote
0 answers

How to run ASReml via shell in R

I am trying to run a .bat file via R, which contains input files for ASReml. The .bat file is working correctly when executed in the Windows cmd prompt, but I can't get it working via R. The code that I use was working before (some months ago) but…
Marjolein Verweij
  • 59
  • 5
1
vote
0 answers

What does my "Error: factor level [4] is duplicated" in AS-Reml R mean?

For my PhD I am currently running a bivariate variance component estimation in AS Reml R with the following command: Test <- asreml(cbind(Reflex,lakt_1) ~ trait + Betriebsnummer + at(trait,1):gestLength + at(trait,1):halfyear, random=…
Helen Schneider
  • 11
  • 1
1
vote
1 answer

Predictions of missing response values in ASreml/ASreml-R

Does anyone know if there is a simple way of generating predictions from an ASreml-R model for cases where the response variable is missing but all values for fixed effects and levels of random effect factors are known. The generic "fitted" function…
PS1
  • 11
  • 2
1
vote
3 answers

Hiding information when fitting model in asreml-R

I'm using asreml-R package for my analysis. Whenever I stored a fitted model object asreml-R gives some extra information which I want to hide. See the code and the information below: library(asreml) dat <- data.frame(y=rnorm(20),x=seq(1,20)) ex.asr…
MYaseen208
  • 22,666
  • 37
  • 165
  • 309
0
votes
2 answers

How to convert a text (with the same colunms) into a table?

[1] genos IDV_V 24 0.506472 14.0206 1.17 0 P [2] Lcrop IDV_V 6 0.768434E-06 0.212724E-04 0.00 0 B [3] Lgenos IDV_V 24 0.768434E-06 0.212724E-04 0.00 0 B [4] …
Aoxiang Xin
  • 87
  • 4
0
votes
0 answers

Error in installing R package asreml on mac

I'm not too familiar with R and are using it for Uni, and we need this asreml package for one of our classes. When I try installing it this error comes up: Error: package or namespace load failed for ‘asreml’: .onLoad failed in loadNamespace() for…
Cami
  • 1