par(mfrow=c(2,2)) # 2x2 matrix of graphs
 points(parana)

#Variogram estimation
 parana.cloud <- variog(parana,option="cloud")
 plot(parana.cloud)
 parana.bin <- variog(parana,uvec = max(parana.cloud$u)*seq(0,1,l=11))
 plot(parana.bin)

#Check for significant spatial correlation
 env.mc <- variog.mc.env(parana,obj.var=parana.bin)
 plot(parana.bin,envelope = env.mc)