<style type='text/css'>.TerraceMsg { font-size: 12px; font-family:Dotum, Arial, Verdana, Sans-Serif;}.Bold { font-weight: bold; }</style><div class='TerraceMsg'><p>Hello, I have a csv file named "seoul1to7" contains hourly pm10 concentration for around 105 sites collected for 7 days. <a href="https://drive.google.com/file/d/0ByY3OAw62EShdEdqcUVHVTNabnM/view?usp=sharing" target="_self">Please download from here</a></p><p> </p><p> By using this dataset, I have plotted variogram for every hour of 7 days (total 161 variogram) using automap package in R. From every variogram, I can see that the range, sill and nugget value is written on the variogram plot. But I want to extract all the psill,range and nugget value in a table for every individual variogram (hour). Is it possible in R? </p><p> </p><p> [Actually, in future I want to make some graph using those range, sill and nugget to compare the changes]. I wrote the following code to plot 161 variogram :</p><p> </p><p>library(sp)<br>library(gstat)<br>library(rgdal)<br>library(automap)<br>library(latticeExtra)<br>library(ggplot2)</p><p>seoul1to7<-read.csv("seoul1to7.csv")<br>seoul1to7[seoul1to7==0] <-NA<br>seoul1to7 <- na.omit(seoul1to7)<br>seoul1to7_splt<-split(seoul1to7,seoul1to7$time) #a list contain 161 dataframe<br>seq(seoul1to7_splt)<br>names(seoul1to7_splt)</p><p>hours<-as.POSIXct(names(seoul1to7_splt), format="%Y%m%d%H")<br>hours</p><p>vars<-lapply(seq(seoul1to7_splt), function(i)<br>{<br>  dat<-seoul1to7_splt[[i]]  <br>  coordinates(dat)<-~LON+LAT<br>  proj4string(dat) <- "+proj=longlat +datum=WGS84" <br>  dat <- spTransform(dat, CRS("+proj=utm +north +zone=52 +datum=WGS84"))<br>  variogram<-autofitVariogram(log(PM10)~1,dat, model="Sph")<br>  plt<- plot(variogram,plotit=F,asp=1)<br>  plt<- update(plt, main=paste("Variogram for", as.character(hours)[i]))<br>  <br>  return(plt)<br>})<br>vars[[1]]</p><p>.</p><p>.</p><p><br>va
rs[[161]]</p><p> </p><p>If you have any further query please let me know. Thanks in advance.</p><p> </p><p>Orpheus </p></div><br><br/><br/><table border=0 cellpadding=0 cellspacing=0 ><tr><td style='vertical-align: top;font-size:9pt;'><div jquery1411706373835="566"><br></div></td></tr></table><br>
<div id='TMSMDN' style="background:url('http://mail.gist.ac.kr:80/mail/receiveMDN.do?mdnData=oGhEpHGdw7UD%2BvJWe%2BEHzgqOXcw4myi2waJaG4Pju9aaVlyyA2XSFYDXoGEXi1VizG0Zfz9uUSUP%0AaHo7ZvtPfIHsGfzaLxo7H81NhQPURjYYfBt9NM0iRgc%2Bf07TGk0%2FGSq2Gt4%2BuJy5dzEzBEbOVQ%3D%3D%0A')"></div>