[R-sig-Geo] Small problem in extraction of a table of range, sill and nugget in loop

Uzzal uzzal at gist.ac.kr
Fri Sep 11 14:15:15 CEST 2015


Hello, I have a csv file named "seoul1to7" contains hourly pm10 concentration for around 105 sites collected for 7 days. Please download from here  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?   [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 : library(sp)
library(gstat)
library(rgdal)
library(automap)
library(latticeExtra)
library(ggplot2)seoul1to7hoursvars..
vars[[161]] my problem starts from here: Now I want to extract the table of sill, range and Nugget.So I wrote this code table But as output here I got this type of table below


 
  
  
  
  
  
  
  
  
  
  
 
 
  
  
  
  
  
  
  
  
  
  
 

  sill
  
  range
  
  nugget
  
  Sill.1
  
  Range.1
  
  Nugget.1
  
  Sill.2
  
  Range.2
  
  Nugget.2
  
  ???
  
   
  
   
  
   
  
   
  
   
  
   
  
   
  
   
  
   
  
   
  



But I want this type of table below.


 
  
  
  
  
 
 
  
  
  
  
 
 
  
  
  
  
 
 
  
  
  
  
 
 
  
  
  
  
 
 
  
  
  
  
 
 
  
  
  
  
 

  hours
  
  Psill
  
  Range
  
  Nugget
  
  2012030101
  
   
  
   
  
   
  
  2012030102
  
   
  
   
  
   
  
  .
  
   
  
   
  
   
  
  .
  
   
  
   
  
   
  
  .
  
   
  
   
  
   
  
  2012030731
  
   
  
   
  
   
  

What correction I should do in my loop? If you have any further query please let me know. Thanks in advance. Orpheus  




-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://stat.ethz.ch/pipermail/r-sig-geo/attachments/20150911/3d36ab4c/attachment.html>


More information about the R-sig-Geo mailing list