Comments (3)
I think this is the expected behavior. eta=1 is uniform over matrices, not correlations. So even with K=3, you'll see concentration of mass around rho=0, just because it's not possible to build valid matrices with independent correlations.
library(rethinking)
R <- rlkjcorr(1e4,K=3,eta=1)
plot( R[,1,2] , R[,1,3] , col=col.alpha("black",0.2) , pch=16 )
The Stan functions behave the same way:
mc <- "
parameters{
corr_matrix[3] Rho;
}
model{
Rho ~ lkj_corr(1);
}"
m <- stan(model_code=mc,chains=1,warmup=100,iter=1e4)
R <- extract.samples(m)$Rho
plot( R[,1,2] , R[,1,3] , col=col.alpha("black",0.2) , pch=16 )
from rethinking.
Indeed, even in 3 dimensions the space of possible correlation
matrices looks like a stuffed Dorito rather than a full cube.
On Thursday, January 21, 2016, Richard McElreath [email protected]
wrote:
—
Reply to this email directly or view it on GitHub
#17 (comment).
from rethinking.
Ah, you're right. I mis-remembered or misunderstood the Stan documentation. Sorry for the false alarm.
Hope you're enjoying Germany!
from rethinking.
Related Issues (20)
- Random Effect in Zero-Inflated Model
- Rethink package is NOT getting installed HOT 8
- Sampling issue for the Cat Adoptions Survival model.
- Random walk in ulam HOT 2
- problem in precis command
- Can multiple outcomes be regressed on the same linear model?
- ulam fits missing @stanfit slot? HOT 2
- WAIC not computing from ulam() model, even with log_lik = TRUE HOT 4
- Problem installing cmdstanr HOT 1
- ulam fits giving odd errors HOT 4
- errata 2nd Edition page 130
- log_lik fails with matrices in models fit by ulam HOT 6
- Where is sim_globe, and the rest of the code shown in the videos? HOT 2
- Panda_nuts Nut-Cracking Frequency? HOT 6
- The garden of forking data - 2 blue from 1? HOT 3
- quap allows incorrect(?) definition with categorical variable HOT 1
- Consistent syntax error on map2stan/ulam issues HOT 2
- Some typos here and there
- rethinking cannot be installed: Rcpp fails to compile HOT 1
- sim() not working for models with simplex
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from rethinking.