Displaying 1 result from an estimated 1 matches for "x_red".
Did you mean:
a_red
2008 Jan 25
1
increasing speed for permutations of glm
...ake weeks to run, unless I try something else.
Could you give me any suggestions on how to improve the speed of the
code below, or any general suggestions on how I may accomplish this task.
Thanks for your time,
Juliet
To run this code, first enter the model matrices (matrices given at bottom):
# X_red <- as.matrix(read.table("clipboard",header=F),nrow=18,byrow=T)
# X_full <- as.matrix(read.table("clipboard",header=F),nrow=18,byrow=T)
library(combinat)
# make some example data; the actual data is 700x800
myData <- matrix(sample(c(1:3),500,replace=TRUE),nrow=100,ncol...