Hi,
I am having a number of problems using rpy with locfit. One example of a
problem is the following code:
import re
import rpy
from scipy import rand
from math import sin
r = rpy.r
rpy.set_default_mode (rpy.NO_CONVERSION)
r.library("locfit")
x = []
y = []
for i in range(0, 1000):
x.append(i)
y.append(sin(i/300) + rand())
r.assign("x", x)
r.assign("y", y)
alpha = 63.0
r.assign("alpha", alpha)
# Calculation 1
residuals = r("residuals(locfit(y~x, alpha = alpha, deg = 1, maxk=1000))")
rpy.set_default_mode(rpy.BASIC_CONVERSION)
print "residuals calculation 1"
print rpy.r.summary(residuals)
# Calculation 2
rpy.set_default_mode(rpy.NO_CONVERSION)
test1 = r("locfit(y~x, alpha = alpha, deg = 1, maxk=1000)")
residuals = r.residuals(test1)
rpy.set_default_mode (rpy.BASIC_CONVERSION)
print "residuals calculation 2"
print rpy.r.summary(residuals)
# Calculation 3
rpy.set_default_mode(rpy.NO_CONVERSION)
fit = r.locfit(r("y~x"), alpha = alpha, deg = 1, maxk=1000)
residuals = r.residuals(fit)
rpy.set_default_mode(rpy.BASIC_CONVERSION)
print rpy.r.summary(residuals)
The program crashes in calculation 3 when I try to pass fit to r.residuals.
I include the print out from Eclipse.
I would appreciate any help.
Thanks,
Jurek.
RHOME= C:\ProgramFiles\R\R-2.2.1pat
RVERSION= 2.2.1
RVER= 2021
RUSER= C:\Documents and Settings\jurek
Loading the R DLL C:\ProgramFiles\R\R- 2.2.1pat\bin\R.dll .. Done.
Loading Rpy version 2021 .. Done.
Creating the R object 'r' .. Done
Loading required package: akima
Loading required package: lattice
locfit 1.5-3 2006-04-06
residuals calculation 1
{'Min.': -1.226, '1st Qu.': -0.29999999999999999, '3rd Qu.':
0.35510000000000003, 'Median': 0.019230000000000001, 'Max.':
0.94179999999999997, 'Mean': 5.8970000000000005e-007}
residuals calculation 2
{'Min.': -1.226 , '1st Qu.': -0.29999999999999999, '3rd Qu.':
0.35510000000000003, 'Median': 0.019230000000000001, 'Max.':
0.94179999999999997, 'Mean': 5.8970000000000005e-007}
Traceback (most recent call last):
File "C:\cygwin\home\jurek\RPYProblemExample.py", line 37, in ?
residuals = r.residuals(fit)
rpy.RException: Error in m[[1]] == "locfit.raw" : comparison (1) is possible
only for atomic and list types
-------------------------------------------------------------------------
Take Surveys. Earn Cash. Influence the Future of IT
Join SourceForge.net's Techsay panel and you'll get the chance to share your
opinions on IT & business topics through brief surveys - and earn cash
http://www.techsay.com/default.php?page=join.php&p=sourceforge&CID=DEVDEV
_______________________________________________
rpy-list mailing list
rpy-list@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/rpy-list
