• Source 
    1. futures = read.csv("D:/Dropbox/Public/futures.csv",header=T,sep=",")
    2. n = names(futures)[3:139]
    3. fmla = as.formula(paste(names(futures)[2],"~",paste0(n,collapse="+")))
    4.  
    5. model1<-glm(formula=fmla, data=futures, family=binomial(link="logit"), na.action=na.exclude)
    6. summary(model1)
    7. summary(model1)$coefficients[which(summary(model1)$coefficients[,4]<0.05),]
    8. ##想把p<0.05的變數列出來
    9.  
    10. p = predict(model1,futures,type = "response")
    11. ##這裡我要如何設定model1使他只會吃p<0.05的變數