R/scatter.R

In mgraph: Graphing map attributes and non-map variables in R

inscatter <- function(source,layer='',attributes,type='',label='',col='')
{

#read data
if(layer=="gisobject"||is.null(layer)||layer=="gisdata"||layer==""){
object <- source

}
else{

file=c("[.]csv$","[.]txt$","[.]tab$","[.]dat$")
i=1
nonspatial=0
while (i <= length(file)) {
log=grep(file[i],source)
if(length(log)>0)
{
file.type=file[i]
thisfile=source
nonspatial=nonspatial+1
}
else{
log=grep(file[i],layer)
if(length(log)>0){
file.type=file[i]
thisfile=paste(source,layer,sep = "/")
nonspatial=nonspatial+1
}
}
i=i+1
}
if(nonspatial<1){
if(layer==""||layer=="nofile"||layer=="notfile"||layer=="nonfile"){
object=source
}
else{
#check input data
if(source==""){
print("Full Folder path is needed:source")
return
}
if(layer==""){
print("Map name is needed: layer")
return
}

object <- readOGR(source, layer)
}
}
else{

log=grep('[.]csv$',file.type)
if(length(log)>0)
{
file.type=file[i]
object=read.csv(thisfile,header=TRUE)

}
else{
object=read.table(thisfile,header=TRUE)
}

}
}
#scatter region
#split
attribute=strsplit(attributes,",")
attribute1=attribute[[1]][1]
attribute2=attribute[[1]][2]

slice1=object[[attribute1]]
slice2=object[[attribute2]]
slices=object[[attribute2]]

main=label
xlab=toupper(attributes)
#default label
if(label==""){
main=paste("Scatter plot of", attributes)
}

#label only names without %

#define colours
colour=col
if(col==""){
colour="black"
}
ylab=attribute2
#construct scatter
if(type=="map"||type=="spatial"){
if(attribute2!=""){
#par(mfrow=c(1,2)) 
lists=c(attribute1,attribute2)
plot(spplot(object,lists,scales = list(draw = TRUE)))
plot(spplot(object,attribute1,scales = list(draw = TRUE)))
}
plot(spplot(object,attribute2,scales = list(draw = TRUE)))
}
else{
plot(slice1,slice2, xlab=attribute1,ylab=attribute2,col=colour,  
main=main,pch=19)
}
#lines(slice1,slice2)
data=object[[attributes]]
list(data=data,table=slices,source=source,layer=layer,
attributes=attributes,type=type,label=label,x=slice1,y=slice2,colour=col)
}

#generic function
scattermap<-function(source,layer='',attributes,type='',label="",col='') UseMethod ("scattermap")

#default function
scattermap.default<-function(source,layer='',attributes,type='',label="",col='')
{

factor<-inscatter(source,layer,attributes,type,label,col)

factor$call<-match.call()

class(factor)<-"scattermap"
factor
}

#summary function
summary.scattermap<-function(object,...)
{
x<-object
cat("Call:\n")
print(x$call)
cat("\n correlation between ")
print(x$attributes)
cor(x$x,x$y)
cat("\nLinear relationship\n")
summary(lm(x$x~x$y))
}

#print function
print.scattermap<-function(x,...)
{
cat("Call:\n")
print(x$call)
cat("\nCall: Data\n")
print(x$data)

cat("\n Table  \n")
print(x$table)

}

#plot function
plot.scattermap<-function(x,...)
{
scattermap(source=x$source,layer=x$layer,attributes=x$attributes,type=x$type,col=x$colour,label=x$label)

}

#graph=scattermap("E:/geokings/advancegis/R-process", type="map","KwabibiremAtiwa",attributes='Age,FArea',col='black',label="")