var.init for generic effect

[Victor Papin]

Hello,

I have a simple model ran with breedR:

model <- remlf90(fixed = pheno ~ 1, generic = list(genetic= list(inc.mat, cov.mat)), data = mydata)

I want to specify initial variances by adding:

model <- remlf90(fixed = pheno ~ 1, generic = list(genetic= list(inc.mat, cov.mat)), list(generic_genetic = 2, residuals = 10), data = mydata) 

but I got this error message :

Error in check_var.ini(var.ini, random, responsem) : Some initial variances missing in var.ini. Please specify either all or none.

I don't seem to be naming my effect correctly,  especially the term "generic_genetic". I tried with "generic", "genetic", "random" but it never works. Can you tell me how this effect should be named properly in the var.init term ?

Best,
Victor

[Facundo Muñoz]

Dear Victor,

Include the specification for the initial variance of the generic component as a third element of the list with the incidence and covariance matrices, like so:

model <- remlf90(

  fixed = pheno ~ 1,

  generic = list(genetic = list(inc.mat, cov.mat, var.ini = 2)),

  var.ini = list(residuals = 10),

  data = mydata

)

Hope it helps.

ƒacu.-

--
Report issues at https://github.com/famuvie/breedR/issues
---
You received this message because you are subscribed to the Google Groups "breedR" group.
To unsubscribe from this group and stop receiving emails from it, send an email to breedr+un...@googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/breedr/72888f31-8150-4643-a73b-60bd6a6728d6n%40googlegroups.com.