stefaneng · August 29, 2015 14:23
diff --git a/sig_variables_glm.R b/sig_variables_glm.R
 # From http://stackoverflow.com/a/16154053/1351718

 # x is a 10 x 3 matrix
 x = matrix(rnorm(3*10), ncol=3)
 y = rnorm(10)
 res = glm(y~x)
 # ignore the intercept pval
 summary(res)$coeff[-1,4] < 0.05


 ## This answer http://stackoverflow.com/a/16154609/1351718 show how this can be automated to create a new glm
	# From http://stackoverflow.com/a/16154053/1351718

	# x is a 10 x 3 matrix
	x = matrix(rnorm(3*10), ncol=3)
	y = rnorm(10)
	res = glm(y~x)
	# ignore the intercept pval
	summary(res)$coeff[-1,4] < 0.05


	## This answer http://stackoverflow.com/a/16154609/1351718 show how this can be automated to create a new glm