R-Scripte:


library(RODBC)

### Example1: observed vs. predicted data

ch10<-odbcConnectAccess("C:/.../.../CCB.mdb")

values<-sqlQuery(ch10,"SELECT rs_table.FL_ID, rs_table.YEAR_, rs_table.MEAS_VALUE, rs_table.SIM_VALUE FROM rs_table INNER JOIN field_description ON rs_table.FL_ID = field_description.FL_ID WHERE (((rs_table.M_IX)=7) AND ((rs_table.YEAR_NUM)>0) AND ((rs_table.SIM_VALUE)>-99)) ORDER BY rs_table.FL_ID, rs_table.YEAR_")
obs<-values$MEAS_VALUE
pre<-values$SIM_VALUE

abb <- layout(matrix(1:1,1,1))
par(oma=c(0,0,0,0)+0.1, mar=c(4,4.5,4,1)+0.1)
layout.show(abb)

plot(pre, obs, ylim=c(0,6),xlim=c(0,6), cex.axis=1.2, cex.lab=1.2, cex.main=1.2 , xlab="predicted values [M%]", ylab="observed values [M%]" )
grid(ny=NULL,nx=NA,col = "gray", lty = "dotted")
grid(ny=NA,nx=NULL,col = "gray", lty = "dotted")
abline(0,1)
cor.test(pre, obs, method="pearson")
lm.obspre <- lm(obs ~ pre)
b <- round(coefficients(lm.obspre)[1],3)
a <- round(coefficients(lm.obspre)[2],3)
lm.obspre.sum <- summary(lm.obspre)
r.q <- round(lm.obspre.sum$r.squared,3)
text(4.5,2.5, expression(R^2), adj=c(1,1.5))
text(4.51,2.5, label=paste("=", r.q), adj=c(-0.2,2.2))


### Example 2: ANOVA
anova_tab2corg_1<-sqlQuery(ch10,"select 'BA'+str(Kl_BOART) as Kl_BOART1, MW_delta_Corg from rs_table_delta_4anova_fehler")
attach(anova_tab2corg_1)
anova_tab2corg_1$Kl_BOART1<-as.factor(anova_tab2corg_1$Kl_BOART1)
a<-lm(MW_delta_Corg~Kl_BOART1, data=anova_tab2corg_1)
anova(a)
summary(a)

library(stats)
aov1<-aov(MW_delta_Corg~Kl_BOART1, data=anova_tab2corg_1)
TukeyHSD(aov1, "Kl_BOART1",ordered=TRUE, conf.level = 0.95)# TUKEY-Test
plot(TukeyHSD(aov1, "Kl_BOART1"))


### Example3: boxplot soil type
png(filename = "C:/.../.../RMSErel.png", width = 480, height = 480, units = "px", pointsize = 12, bg = "white", res = NA, restoreConsole = TRUE)
boxplot(rmsprz~Kl_BOART1, data=anova_tab2corg_1, xlab="category", ylab=expression(paste(RMSErel," [%]")), ylim=c(0,35), main="RMSE rel ~ soil texture", cex.axis=1.4, cex.lab=1.4)
grid(ny=NULL,nx=NA,col = "gray", lty = "dotted", lwd=3)
dev.off()


### Example4: Shapiro-Wilk test of ME
ME<-sqlQuery(ch10,"select med from ccb_stat_corg")
attach(ME)
shapiro.test(ME$med)