<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 "seoul_3month" contains hourly pm10 concentration for around 105 sites collected for 3 months. <a href="https://drive.google.com/file/d/0ByY3OAw62EShWVlJNFR4ZWU1WW8/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 3 months (total 2208 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?  [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 2208 variogram :</p><p> </p><p>library(sp)<br>library(gstat)<br>library(rgdal)<br>library(automap)</p><p><br>seoul_3month<-read.csv("seoul_3month.csv")<br>seoul_3month[seoul_3month==0] <-NA<br>seoul_3month <- na.omit(seoul_3month)</p><p><br>seoul_3month_hour<-split(seoul_3month,seoul_3month$time) <br>seq(seoul_3month_hour)</p><p> </p><p>vars<-lapply(seq(seoul_3month_hour), function(i)<br>{<br>  dat<-seoul_3month_hour[[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)<br>  plot<- plot(variogram,plotit=FALSE, asp=1)</p><p><br>  return(plot)<br>})<br>vars[[1]]</p><p>.</p><p>.<br>vars[[2208]] </p><p> </p><p>If you have any further query please let me know. Thanks in advance.</p><p> </p><p>Orpheus<br> </p></div><br><br/><br/><table border=0 cellpadding=0 cellspacing=0 ><tr><td style='vertical-align: to
p;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%2BEHzqmIW2tGDc4Qy5wwQejSb0bmTwY32tiKlJWFNa3qDm5RlDzhT9YkWq%2Fg%0Apnpso%2FTSq5z9cGa%2FphYmGU6O2PK8p84Akvk3px%2FzHHzJ3hjnxPW%2BCfkFkpWgaB%2Fhp%2FH4Q%2Bk0yg%3D%3D%0A')"></div>