# add data
data(EuStockMarkets)
d <- as.vector(EuStockMarkets[,1])
# Calculate Z score
summary(d)
summary(Z_org <- z_normalization(d) )
summary(Z_mod <- z_normalization(d, method="modified") )
par(mfrow=c(3,1))
plot(density(d), main="original timeseries")
plot(density(Z_org), main="original z norm")
plot(density(Z_mod), main="modified z norm")
# Check NA behavior, insert some NA's
d[c(100, 500, 1000, 1400)] <- NA
length(d)
length( z_normalization(d, na="keep") )
length( z_normalization(d, na="remove") )
Run the code above in your browser using DataLab