Justin Silverman

Wonderful! (1) you are correct (2) I think you are write calling this `\upsilon` would be better. I had distinguished gamma and upsilon only because they are different priors and...

Great Idea! I will do this. Justin > On Nov 5, 2020, at 11:47 AM, Michael McLaren wrote: > > > (Though perhaps if there are already published papers /...

Hi Daniel, That sounds like a great idea. Maltipoo has been a work in progress for a long time. Honestly I think it works well (may not be as fast...

U is specified by stacking each variance component. If each variance component is N x N (the number of samples) and you have p variance components... then U is Np...

I would probably say some correlation ρ between biological replicates 1 on the diagonal and 0 other wise. There are many options though. Justin > On Sep 25, 2020, at...

Very sorry for the delay. Have you checked the dimensions of the other parts? e.g., are you sure that nrow(X) == ncol(Theta) == ncol(U)? Justin

Oh shit! I said the wrong thing, U has dimension Q x Q - number of covariates (its just for many of these models Q == N). My apologies! Same...

So in general mixed effects models (and associated variance components do not need a distinct random effect for each sample). You can, for example, have multiple samples sharing a random...