Nothing
#' @export
######## DO NOT CALCULATE CP******************************************
modelval<-function(building.set,response.building, prediction.set, response.prediction){
building.set<-as.matrix(building.set)
response.building<-as.matrix(response.building)
prediction.set<-as.matrix(prediction.set)
response.prediction<-as.matrix(response.prediction)
model1<-lm(response.building~building.set)
model2<-lm(response.prediction~prediction.set)
smodel1<-summary(model1)$coeff[,1:2]
smodel2<-summary(model2)$coeff[,1:2]
#leaps(model.buildi11g.set,y,method = 'Cp')
#cp1<-cpc(model1)
#cp2<-cpc(model2)
p1<-pressc(model1)
p2<-pressc(model2)
a1<-AICp(model1)
a2<-AICp(model2)
s1<-SBCp(model1)
s2<-SBCp(model2)
r1<-summary(model1)$r.squared
r2<-summary(model2)$r.squared
r2a1<-summary(model1)$adj.r.squared
r2a2<-summary(model2)$adj.r.squared
Model.Training<-c(smodel1[,1],smodel1[,2],deviance(model1)/df.residual(model1),r1,r2a1,a1,s1,p1)
Model.Validation<-c(smodel2[,1],smodel2[,2],deviance(model2)/df.residual(model2),r2,r2a2,a2,s2,p2)
v<-cbind(Model.Training,Model.Validation)
namE<-paste('b',0:ncol(building.set),sep = "")
namSE<-paste('SE b',0:ncol(building.set),sep = "")
rownames(v)<-c(namE,namSE,'MSE','R2','R2.adj','AICp','SBCp','PRESSp')
return(v)
}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.