? stato filtrato un testo allegato il cui set di caratteri non era indicato... Nome: non disponibile URL: <https://stat.ethz.ch/pipermail/r-help/attachments/20130501/294047df/attachment.pl>
rpy2 postgres qgis problem
3 messages · Enzo Cocca, R. Michael Weylandt
On Wed, May 1, 2013 at 3:05 PM, Enzo Cocca <enzo.ccc at gmail.com> wrote:
Hello every body,
I am using rpy2_2.0.8 with postgres and Qgis.
The code that I wrote is the following:
def on_calcola_pressed (self):
# bottone per calcoli statistici
import rpy2
import rpy2.robjects as robjects
import rpy2.robjects as ro
Redundant no?
r = robjects.r
robjects.r
Why?
r.library("RPostgreSQL")
r.library("gstat")
r.library("zoo")
rcmd_drv = r["dbDriver"]
rcmd_con = r["dbConnect"]
rcmd = r["dbListTables"]
rcmd_archezoology_table = r["dbReadTable"]
png = r.png
drv = rcmd_drv("PostgreSQL")
con = rcmd_con(drv, host="127.0.0.1", dbname="pyarchinit",
port="5432", password="enzo", user="postgres")
archezoology_table =
rcmd_archezoology_table(con,"archeozoology_table")
r.png(file="/home/enzo/Scrivania/prova.png",width=512,height=512)
plot = r.plot
zoo = r.zoo
('VGM_PARAM_A3 <- gstat(id="bos_bison",
formula=combusto~1,locations=~coord_x+coord_y, data=archezoology_table,
nmax = 10) VGM_PARAM_A3 <- gstat(VGM_PARAM_A3, "calcinati", strie~1,
locations=~coord_x+coord_y, archezoology_table, nmax = 10)VGM_PARAM_A3 <-
gstat(VGM_PARAM_A3, "camoscio", cervo~1,
locations=~coord_x+coord_y,archezoology_table, nmax = 10) VGM_PARAM_A3 <-
gstat(VGM_PARAM_A3, model=vgm(1, "Sph", 5, 0), fill.all=TRUE ESV_A3 <-
variogram(VGM_PARAM_A3, cutoff=9) VARMODEL_A3 = fit.lmc(ESV_A3,
VGM_PARAM_A3)')
plot(zoo('ESV_A3, model = VARMODEL_A3,xlab=,ylab=,pch=20,
cex=0.7, col="red",main="Linear Model of Coregionalization for A3
variables"'))
As this appeared in my mailer, it's not valid python or R; perhaps you could try resending not in HTML? It's rather dificult to debug as is. Even better would be to see if you can get it working solely within R -- post that code -- and then work on translating into rpy2. Cheers, MW
The error is the following:
res = super(RFunction, self).__call__(*new_args, **new_kwargs)
RRuntimeError: Error in plot.window(...) : invalid 'ylim' value
Do you have a suggestion or examples about how to call a function with the
gistat module on a postgress tab using rpy2?
Thank you,
Enzo
--
Enzo Cocca (PhD Candidate)
Research Fellow
Universit? di Napoli "L'Orientale"
mail: enzo.ccc at gmail.com
cell: +393495087014
[[alternative HTML version deleted]]
______________________________________________ R-help at r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
? stato filtrato un testo allegato il cui set di caratteri non era indicato... Nome: non disponibile URL: <https://stat.ethz.ch/pipermail/r-help/attachments/20130502/2df074e0/attachment.pl>