# for reproducibility: change seed to change noise
set.seed(3)
plot(w.spec$freq,w.spec$spec,t="l",ylim = c(0,1.2),
xlab="frequency (cycles/sample interval)",ylab="spectral density")
plot(w.spec$freq,w.spec$spec,t="l",ylim = c(0,1.2),
xlab="frequency (cycles/sample interval)",ylab="spectral density")
.libpaths()
.libPaths()
.libPaths()
# for reproducibility: change seed to change noise
set.seed(3)
# create noise time series
w <- rnorm(1000)
plot(as.ts(w),ylab="White noise")
# verify w is zero mean
mean(w)
var(w)
hist(w,main="",col="white")
# verify w is uncorrelated
acf(w,main="",lag=20,ylim=c(-0.2,1))
# verify w is really white
w.spec <- spectrum(w,log="no",span=200,plot=FALSE)
plot(w.spec$freq,w.spec$spec,t="l",ylim = c(0,1.2),
xlab="frequency (cycles/sample interval)",ylab="spectral density")
hist(w,main="",col="white")
# for reproducibility: change seed to change noise
set.seed(3)
# create noise time series
w <- rnorm(1000)
plot(as.ts(w),ylab="White noise")
# verify w is zero mean
mean(w)
var(w)
hist(w,main="",col="white")
# verify w is uncorrelated
acf(w,main="",lag=20,ylim=c(-0.2,1))
# verify w is really white
w.spec <- spectrum(w,log="no",span=200,plot=FALSE)
plot(w.spec$freq,w.spec$spec,t="l",ylim = c(0,1.2),
xlab="frequency (cycles/sample interval)",ylab="spectral density")
install.packages("agricolaeplotr")
library(agricolaeplotr)
library("agricolae")
library(agricolae)
library(agricolaeplotr)
library(ggplot2)
library(ggplot2)
library(agricolaeplotr)
library(agricolae)
library(ggplot2)
help("agricolae-package")
help("design.crd")
esempio <- rnorm(1000,0,1)
plot(esempio)
hist(esempio)
esempio <-runif(1000,0,1)
esempio1<-runif(1000,0,1)
plot(esempio1)
hist(esempio1)
esempio1<-runif(120)
plot(esempio1)
hist(esempio1)
esempio1
esempio1<-runif(30)
plot(esempio1)
hist(esempio1)
esempio1
esempio1<-runif(30)
esempio1
set.seed(1)
esempio1<-runif(30)
plot(esempio1)
hist(esempio1)
esempio1
set.seed(1)
esempio1<-runif(30)
plot(esempio1)
hist(esempio1)
esempio1
1/120
1/30
library(agricolaeplotr)
library(agricolae)
library(ggplot2)
esempio <- rnorm(1000,0,1)
plot(esempio)
hist(esempio)
set.seed(1)
esempio1<-runif(30)
plot(esempio1)
hist(esempio1)
esempio1
trt<-c(3,2) # factorial 3x2
outdesign <- design.ab(trt, r=3, serie=2,design = 'crd')
trt<-c(3,2) # factorial 3x2
outdesign <- design.ab(trt, r=3, serie=2,design = 'crd')
head(outdesign$book,10)
plot_design.factorial_crd(outdesign,ncols=7,nrows=3, width = 1, height = 1)
library(agricolaeplotr)
library(agricolae)
library(ggplot2)
# esempio <- rnorm(1000,0,1)
# plot(esempio)
# hist(esempio)
# set.seed(1)
# esempio1<-runif(30)
# plot(esempio1)
# hist(esempio1)
# esempio1
trt<-c(3,2) # factorial 3x2
outdesign <- design.ab(trt, r=3, serie=2,design = 'crd',
width = 1, height = 1, reverse_y = TRUE)
outdesign <- design.ab(trt, r=3, serie=2,design = 'crd', width = 1, height = 1, reverse_y = TRUE)
library(agricolaeplotr)
library(agricolae)
library(ggplot2)
# esempio <- rnorm(1000,0,1)
# plot(esempio)
# hist(esempio)
# set.seed(1)
# esempio1<-runif(30)
# plot(esempio1)
# hist(esempio1)
# esempio1
trt<-c(3,2) # factorial 3x2
outdesign <- design.ab(trt, r=3, serie=2,design = 'crd')
head(outdesign$book,10)
plot_design.factorial_crd(outdesign,ncols=7,nrows=3, width = 1, height = 1)
trt<-c(3,2) # factorial 3x2
outdesign <- design.ab(trt, r=3, serie=2,design = 'crd')
head(outdesign$book,10)
plot_design.factorial_crd(outdesign,ncols=7,nrows=3, width = 1, height = 1)
plot_design.factorial_crd(outdesign,ncols=7,nrows=3, width = 1, height = 1,reverse_y = TRUE)
plot_design.factorial_crd(outdesign,ncols=7,nrows=3, width = 6, height = 1,reverse_y = TRUE)
plot_design.factorial_crd(outdesign,ncols=7,nrows=3, width = 6, height = 1,reverse_y = FALSE)
plot_design.factorial_crd(outdesign,ncols=7,nrows=3, width = 6, height = 1,reverse_y = FALSE)
outdesign <- design.ab(trt, r=3, serie=2,design = 'crd')
head(outdesign$book,21)
plot_design.factorial_crd(outdesign,ncols=7,nrows=3, width = 6, height = 1,reverse_y = FALSE)
library(sp)
library(gstat)
library(quantmod)
# =========================================
setwd("C:/RPA/code/spatial_analysis")
setwd("~/Didattica/R_class_5/exercises/spatial_analysis")
library(sp)
library(gstat)
library(quantmod)
help(sp)
help(gstat)
help("quantmod")
help(quantmod)
# =========================================
setwd("~/Didattica/R_class_5/exercises/spatial_analysis")
# =========================================
# load dat
geodata1 = read.csv("data/transect_1.csv")
geodata2 = read.csv("data/transect_2.csv")
edit(geodata1)
plot(geodata1$s,geodata1$z,xlab=("Distance [m]"),ylab=("Silt content [-]"),ylim=c(0,0.5))
plot(geodata2$s,geodata2$z,xlab=("Distance [m]"),ylab=("Silt content [-]"),ylim=c(0,0.5))
mu1<- mean(geodata1$z) #mean
V1<- var(geodata1$z) #variance
S1<- sd(geodata1$z) #standard deviation
mu2<- mean(geodata2$z) #mean
V2<- var(geodata2$z) #variance
S2<- sd(geodata2$z) #standard deviation
output -> c(mu1,V1,S1,mu2,V2,S2)
output <- c(mu1,V1,S1,mu2,V2,S2)
output
geodata1$lag1<-Lag(geodata1$z, k = 1)
geodata1$lag1<-as.numeric(geodata1$lag1)
geodata1$lag2<-Lag(geodata1$z, k = 2)
geodata1$lag2<-as.numeric(geodata1$lag2)
geodata1$lag3<-Lag(geodata1$z, k = 3)
geodata1$lag3<-as.numeric(geodata1$lag3)
geodata1$lag4<-Lag(geodata1$z, k = 4)
geodata1$lag4<-as.numeric(geodata1$lag4)
geodata1$lag5<-Lag(geodata1$z, k = 5)
geodata1$lag5<-as.numeric(geodata1$lag5)
edit(geodata1)
geodata2$lag1<-Lag(geodata2$z, k = 1)
geodata2$lag1<-as.numeric(geodata2$lag1)
geodata2$lag2<-Lag(geodata2$z, k = 2)
geodata2$lag2<-as.numeric(geodata2$lag2)
geodata2$lag3<-Lag(geodata2$z, k = 3)
geodata2$lag3<-as.numeric(geodata2$lag3)
geodata2$lag4<-Lag(geodata2$z, k = 4)
geodata2$lag4<-as.numeric(geodata2$lag4)
geodata2$lag5<-Lag(geodata2$z, k = 5)
geodata2$lag5<-as.numeric(geodata2$lag5)
edit(geodata2)
##Plot scatter plots
str(geodata1)
plot(geodata1$z,geodata1$lag1,xlim=c(0.2,0.4),ylim=c(0.2,0.4),xlab=("x+h"),ylab=("x"))
abline(a=0,b=1,col="blue")
plot(geodata1$z,geodata1$lag2,xlim=c(0.2,0.4),ylim=c(0.2,0.4),xlab=("x+2h"),ylab=("x"))
abline(a=0,b=1,col="blue")
plot(geodata1$z,geodata1$lag3,xlim=c(0.2,0.4),ylim=c(0.2,0.4),xlab=("x+3h"),ylab=("x"))
abline(a=0,b=1,col="blue")
plot(geodata1$z,geodata1$lag4,xlim=c(0.2,0.4),ylim=c(0.2,0.4),xlab=("x+4h"),ylab=("x"))
abline(a=0,b=1,col="blue")
plot(geodata1$z,geodata1$lag5,xlim=c(0.2,0.4),ylim=c(0.2,0.4),xlab=("x+5h"),ylab=("x"))
abline(a=0,b=1,col="blue")
corlag1<-cor(geodata1$z,geodata1$lag1,use="complete.obs")
corlag2<-cor(geodata1$z,geodata1$lag2,use="complete.obs")
corlag3<-cor(geodata1$z,geodata1$lag3,use="complete.obs")
corlag4<-cor(geodata1$z,geodata1$lag4,use="complete.obs")
corlag5<-cor(geodata1$z,geodata1$lag5,use="complete.obs")
correlations1<-c(corlag1,corlag2,corlag3,corlag4,corlag5)
correlations1
str(geodata1)
plot(geodata1$z,geodata1$lag1,xlim=c(0.2,0.4),ylim=c(0.2,0.4),xlab=("x"),ylab=("x+h"))
abline(a=0,b=1,col="blue")
plot(geodata1$z,geodata1$lag2,xlim=c(0.2,0.4),ylim=c(0.2,0.4),xlab=("x"),ylab=("x+2h"))
abline(a=0,b=1,col="blue")
plot(geodata1$z,geodata1$lag3,xlim=c(0.2,0.4),ylim=c(0.2,0.4),xlab=("x"),ylab=("x+3h"))
abline(a=0,b=1,col="blue")
plot(geodata1$z,geodata1$lag4,xlim=c(0.2,0.4),ylim=c(0.2,0.4),xlab=("x"),ylab=("x+4h"))
abline(a=0,b=1,col="blue")
plot(geodata1$z,geodata1$lag5,xlim=c(0.2,0.4),ylim=c(0.2,0.4),xlab=("x"),ylab=("x+5h"))
abline(a=0,b=1,col="blue")
str(geodata1)
plot(geodata1$z,geodata1$lag1,xlim=c(0.2,0.4),ylim=c(0.2,0.4),xlab=("x"),ylab=("x+h"))
abline(a=0,b=1,col="blue")
plot(geodata1$z,geodata1$lag2,xlim=c(0.2,0.4),ylim=c(0.2,0.4),xlab=("x"),ylab=("x+2h"))
abline(a=0,b=1,col="blue")
plot(geodata1$z,geodata1$lag3,xlim=c(0.2,0.4),ylim=c(0.2,0.4),xlab=("x"),ylab=("x+3h"))
abline(a=0,b=1,col="blue")
plot(geodata1$z,geodata1$lag4,xlim=c(0.2,0.4),ylim=c(0.2,0.4),xlab=("x"),ylab=("x+4h"))
abline(a=0,b=1,col="blue")
plot(geodata1$z,geodata1$lag5,xlim=c(0.2,0.4),ylim=c(0.2,0.4),xlab=("x"),ylab=("x+5h"))
abline(a=0,b=1,col="blue")
##Plot scatter plots
plot(geodata2$z,geodata2$lag1,xlim=c(0.2,0.4),ylim=c(0.2,0.4),xlab=("x"),ylab=("x+h"))
abline(a=0,b=1,col="blue")
plot(geodata2$z,geodata2$lag2,xlim=c(0.2,0.4),ylim=c(0.2,0.4),xlab=("x"),ylab=("x+2h"))
abline(a=0,b=1,col="blue")
plot(geodata2$z,geodata2$lag3,xlim=c(0.2,0.4),ylim=c(0.2,0.4),xlab=("x"),ylab=("x+3h"))
abline(a=0,b=1,col="blue")
plot(geodata2$z,geodata2$lag4,xlim=c(0.2,0.4),ylim=c(0.2,0.4),xlab=("x"),ylab=("x+4h"))
abline(a=0,b=1,col="blue")
plot(geodata2$z,geodata2$lag5,xlim=c(0.2,0.4),ylim=c(0.2,0.4),xlab=("x"),ylab=("x+5h"))
abline(a=0,b=1,col="blue")
