Merge pull request #17 from master/fix_lda

Bob Carpenter · Bob Carpenter · commit a48bd9c8c247 · 2015-01-08T12:57:32.000-05:00
fix sim-lda script

Thanks --- that one was my fault.  I keep forgetting to clear the environment before testing the scripts.
diff --git a/misc/cluster/lda/sim-lda.R b/misc/cluster/lda/sim-lda.R
@@ -13,6 +13,9 @@ avg_doc_length <- 10;
 doc_length <- rpois(M,avg_doc_length);
 N <- sum(doc_length);
 
+alpha <- rep(1/K,K);
+beta <- rep(1/V,V);
+
 theta <- rdirichlet(M,alpha);
 
 w <- rep(NA,N);
@@ -27,7 +30,4 @@ for (m in 1:M) {
   }
 }
 
-alpha <- rep(1/K,K);
-beta <- rep(1/V,V);
-
 dump(c("K","V","M","N","z","w","doc","alpha","beta"),"lda.data.R");

Original file line number	Diff line number	Diff line change
`@@ -13,6 +13,9 @@ avg_doc_length <- 10;`
`13`	`13`	`doc_length <- rpois(M,avg_doc_length);`
`14`	`14`	`N <- sum(doc_length);`
`15`	`15`
	`16`	`+alpha <- rep(1/K,K);`
	`17`	`+beta <- rep(1/V,V);`
	`18`	`+`
`16`	`19`	`theta <- rdirichlet(M,alpha);`
`17`	`20`
`18`	`21`	`w <- rep(NA,N);`
`@@ -27,7 +30,4 @@ for (m in 1:M) {`
`27`	`30`	`}`
`28`	`31`	`}`
`29`	`32`
`30`		`-alpha <- rep(1/K,K);`
`31`		`-beta <- rep(1/V,V);`
`32`		`-`
`33`	`33`	`dump(c("K","V","M","N","z","w","doc","alpha","beta"),"lda.data.R");`