r - How to plot multiple variogram plots with different title using loop/lapply function? -


i have csv file contains hourly pm10 concentration of 1 march 7 march. please, download here. have plotted variogram (total 161) in loop automap package. 

library(sp) library(gstat) library(rgdal) library(automap) library(latticeextra)  seoul1to7<-read.csv("seoul1to7.csv") seoul1to7[seoul1to7==0] <-na seoul1to7 <- na.omit(seoul1to7) seoul1to7_split<-split(seoul1to7,seoul1to7$time)  seq(seoul1to7_split)    vars<-lapply(seq(seoul1to7_split), function(i) {   dat<-seoul1to7_split[[i]]     coordinates(dat)<-~lon+lat   proj4string(dat) <- "+proj=longlat +datum=wgs84"    dat <- sptransform(dat, crs("+proj=utm +north +zone=52 +datum=wgs84"))   variogram<-autofitvariogram(log(pm10)~1,dat, model="sph")   plot<- plot(variogram,plotit=false, asp=1)    return(plot) }) vars[[1]] vars[[2]]

here, can individual plot vars[[1]],vars[[2]]...etc variogram has same title. now, want plot variogram image different title in loop. want variograms title like, "variogram 2012-03-01 1.00", "variogram 2012-03-01 2.00" .....etc. 

a<-as.posixct(names(seoul1to7_split), format="%y%m%d%h") hours<-substr(a,1,16) hours

i keep desired different title in hours variable "2012-03-01 01:00", "2012-03-01 02:00", "2012-03-01 03:00" ...etc.

how can plot variograms (total 161) different title using loop?

the trick use png , dev.off:

library(automap) data(meuse) coordinates(meuse) = ~x+y  fit_for_model_type = function(model_type) {     return(autofitvariogram(log(zinc) ~ dist, meuse, model = model_type)) }  model_options = c("sph", "exp", "gau", "ste") list_of_variogram_models = lapply(model_options, fit_for_model_type) names(list_of_variogram_models) = model_options lapply(names(list_of_variogram_models), function(x) {     png(filename = sprintf('%s.png', x))     print(plot(list_of_variogram_models[[x]]))     dev.off() })

-

Comments

Post a Comment

Popular posts from this blog

c# - Binding a comma separated list to a List<int> in asp.net web api -

the default binder in web api expecst http://url.com/webapi/report/?pageids=3243&pageids=2365 to bind public ihttpactionresult report(list<int> pageids){ // exciting webapi code} i wish bind http://url.com/webapi/report/?pageids=3243,2365 as running out of space in url get . i have created class public class commaseparatedmodelbinder : system.web.http.modelbinding.imodelbinder { public bool bindmodel(httpactioncontext actioncontext, modelbindingcontext bindingcontext) { //binding in here } } and registered in webapiconfig.cs var provider = new simplemodelbinderprovider( typeof(list<int>), new commaseparatedmodelbinder()); config.services.insert(typeof(modelbinderprovider), 0, provider); i have altered method signature use model binder so public ihttpactionresult report( [modelbinder] list<int> pageids){ // exciting webapi code} however break point in binder not being ...
Read more

how to prompt save As Box in Excel Interlop c# MVC 4 -

i want export excel sheet using excel interloop... i have datatable convert excel sheet... now there way save file using save dialogue box in mvc 4 here code public string datatabletoexcel(datatable dt, string htmlheading, string key) { try { microsoft.office.interop.excel.application excel; microsoft.office.interop.excel.workbook excelworkbook; microsoft.office.interop.excel.worksheet excelsheet; microsoft.office.interop.excel.range excelcellrange; excel = new microsoft.office.interop.excel.application(); excel.visible = false; excel.displayalerts = false; excelworkbook = excel.workbooks.add(type.missing); excelsheet = (microsoft.office.interop.excel.worksheet)excelworkbook.activesheet; excelsheet.name = "sheet1"; dataset ...
Read more

xslt 1.0 - How to access or retrieve mets content of an item from another item? -

my use case have item has link in item. example, item 123456789/152 has metadata field dc.relation.hasversion=123456789/717 . within item view of 123456789/152 , how can retrieve values of metadata of 123456789/717 ? example, want retrieve dc.language.iso value of 123456789/717 123456789/152 . i looked @ related items feature in dspace don't know how metadata displayed in related items list pulled items in list. i using dspace version 5.3 mirage 2 theme . edit based on schweerelos answer, below actual code. using custom metadata field dc.relation.languageversion . basically, want link other version instead of displaying value, display dc.language.iso of other version in item-view.xsl . i've incorporated answer of schweerelos question in code just displaying value of dc.relation.languageversion . sample values of dc.relation.languageversion 10665.1/9843 ; 10665.1/9844 ie not complete uri handle. thanks in advance! actual code <xsl:template name=...
Read more