R
vicky428
这个作者很懒,什么都没留下…
展开
-
savedata.R-20170704
write.table(des1, file = "/Users/vicky/Documents/code/R/des1.txt")write.table(des2, file = "/Users/vicky/Documents/code/R/des2.txt")write.table(des3, file = "/Users/vicky/Documents/code/R/des3.txt...原创 2019-08-31 15:09:30 · 221 阅读 · 0 评论 -
logit2.R-20170913
library(caret) library(MASS)library(glmnet)data=as.matrix(oadata2)colnames(data)=c('y','x1','x2','x17',paste("x",3:16,sep=""))c=10 #验证c次set.seed(5000)r3=matrix(0,c,1)a=matrix(0,c,2)merror=m...原创 2019-09-02 12:46:42 · 191 阅读 · 0 评论 -
nnls.R-20170908
cdata=condata2cdata=as.matrix(cdata)colnames(cdata)=c('y','x1','x2','x7^2','x7','x17',paste("x3",0:1,sep="="),paste("x4",0:6,sep="="), paste("x5","=",0:1,sep=""),paste("x6","="...原创 2019-09-02 12:45:01 · 541 阅读 · 0 评论 -
selectvar.R-20170907
# datap=cbind(condata[,1],pc)# colnames(datap)=c('y',paste("x",1:9,sep=""))# l=lm(y~.,data=data.frame(datap))# summary(l)condata2=condata3condata2=as.matrix(condata2)colnames(condata2)=c('y','...原创 2019-09-02 12:44:29 · 162 阅读 · 0 评论 -
lm.R-20170905
#colnames(condata)=c('y','x1','x2','x17',paste("x",3:16,sep=""))condata=as.matrix(condata)colnames(condata)=c('y','x1','x2','x17',paste("x3",0:1,sep="="),paste("x4",0:6,sep="="), ...原创 2019-09-02 12:43:52 · 122 阅读 · 0 评论 -
dividedata2.R-20170905
library(caret) library("kernlab")data=as.matrix(oadata)colnames(data)=c('y','x1','x2','x17',paste("x",3:16,sep=""))c=1000 #验证c次set.seed(1000)r1=matrix(0,c,1)r2=matrix(0,c,1)r4=matrix(0,c,1)...原创 2019-09-02 12:43:17 · 130 阅读 · 0 评论 -
fit.R-20170905
yinterval=as.matrix(yinterval)f=as.matrix(yinterval[,3:18])xint=seq(5000,155000,10000)par(mfrow=c(2,2))plot(xint,f[4,])m=matrix(0,nrow(f),1)v=matrix(0,nrow(f),1)xmu=matrix(0,nrow(f),1)xsigma...原创 2019-09-02 12:42:29 · 132 阅读 · 0 评论 -
bootstrap.R-20170901
# #bootstrap验证# #返回回归系数向量的自编函数# bs <- function(formula, data, indices){# d <- data[indices,]# fit <- lm(formula,data=d)# return(coef(fit))# }# #自助抽样# set.seed(1234)# a = lm(y~0+...原创 2019-09-02 12:41:47 · 217 阅读 · 0 评论 -
logit.R-20170901
library(caret) data=as.matrix(oadata2)colnames(data)=c('y','x1','x2','x17',paste("x",3:16,sep=""))c=5000 #验证c次set.seed(5000)r1=matrix(0,c,1)r2=matrix(0,c,1)p1=matrix(0,12,c)for (k in 1:c) {...原创 2019-09-02 12:41:08 · 196 阅读 · 0 评论 -
dividedata.R-20170830
#data=cbind(ypart,xcovpart2)#sample(2, nrow(data), replace=T, prob=c(0.7,0.3))library(caret) data=as.matrix(oadata)colnames(data)=c('y','x1','x2','x17',paste("x",3:16,sep=""))c=5000 #验证c次set....原创 2019-09-02 12:40:27 · 159 阅读 · 0 评论 -
covdesign2.R-20170829
data=as.matrix(oadata)colnames(data)=c('y','x1','x2','x17',paste("x",3:16,sep=""))# library(DoE.base)# oa.design(ID=L9.3.4)# oa.design(nfactors=6, nlevels=3)#leave-one-out验证library(caret) #r...原创 2019-09-02 12:39:18 · 138 阅读 · 0 评论 -
covdesign-20170826
x1<-xcovpart[,1] x2<-xcovpart[,2] x3<-xcovpart[,3] x4<-xcovpart[,4] x5<-xcovpart[,5] x6<-xcovpart[,6] x7<-xcovpart[,7] x8<-xcovpart[,8] x9<-xcovpart[,9] x10<-xc...原创 2019-09-02 12:38:03 · 103 阅读 · 0 评论 -
svm.R-20170824
library(e1071)m1<-svm(xcovpart2,ypart)summary(m1)ypart.pre3<-predict(m1,xcovpart2)resi3=abs(ypart-ypart.pre3)/ypartmean(resi3) #误差yall.pre3=predict(m1,newdata=data.frame(xcovall2))m2<...原创 2019-08-31 15:10:28 · 92 阅读 · 0 评论 -
assignment1.R-20170914
#数据预处理data=data.frame(XY_train)colnames(data)=c('y',paste("x",1:34,sep=""))#离散因素数据处理data[,'y']=factor(data[,'y'])data[,'x33']=factor(data[,'x33'])data[,'x34']=factor(data[,'x34'])#转换因子水平为数值le...原创 2019-09-03 12:43:21 · 187 阅读 · 0 评论