Matches in Nanopublications for { <https://w3id.org/ro-id/enrichment/8aa0a807-9dd5-4904-aa8a-198bf65a479c> ?p ?o ?g. }
Showing items 1 to 4 of
4
with 100 items per page.
- 8aa0a807-9dd5-4904-aa8a-198bf65a479c type Sentence assertion.
- 8aa0a807-9dd5-4904-aa8a-198bf65a479c normScore "29.71768202080238" assertion.
- 8aa0a807-9dd5-4904-aa8a-198bf65a479c score "100.0" assertion.
- 8aa0a807-9dd5-4904-aa8a-198bf65a479c name "setwd("/Users/francescodeleo/Desktop/rastafile") require(quantreg) require(splines) require(raster) require(maps) require(Rserve) invasive<-read.csv(file.choose("invasive"),header=TRUE,dec=",",sep=";",quote="\"") names(invasive) str(invasive) sp<-unique(invasive$species) for(k in 1:length(sp)){ png(paste(sp[k],'.png'),height=1000,width=750) par(cex=3,mfrow=c(3,3),lwd=2,mar=c(4,5,4,6)) pel_data<-invasive[invasive$species==sp[k],] pel_data$year<-as.numeric(pel_data$year) coordinates(pel_data)<-c('x.long','y.lat') m<-c(2009,2010,2011,2012,2013,2014,2015) for(i in 1:length(m)){ if(TRUE%in%(pel_data$year<=m[i]&pel_data$year<=m[i])){ pel_dataz<-pel_data[pel_data$year<=m[i]&pel_data$year<=m[i],] rasta<-raster(ext=extent(pel_data)*1.1,res=.33) rasta_fari<-rasterize(pel_dataz,rasta,field='species',fun="count",zlim=c(0,2)) } else rasta_fari<-raster(ext=extent(pel_data)*1.1,res=.33,vals=0) plot(log10(rasta_fari),main=paste('Year',m[i]), col=rev(heat.colors(20)),cex.main=2,xaxt='n',yaxt='n',legend.mar=10,legend.shrink =.95,legend.width=1.5, legend.args=list(text='log10(Number of blooms)',cex=.75,side=4,font=2,line=3.25)) map(add=T,col='grey',fill=T) } dev.off()" assertion.