Displaying 20 results from an estimated 10000 matches similar to: "loess() crashes R on my system"
2005 Jun 15
2
Plotting second axes outside xyplot
Hi all,
I'm trying to find a way to get xyplot to produce a second set of axes outside the right hand side of the graph. This is my progress so far:
EE <- equal.count(ethanol$E, number=9, overlap=1/4)
xyplot(NOx ~ C | EE, data = ethanol,
prepanel = function(x, y) prepanel.loess(x, y, span = 1),
xlab = "Compression Ratio", ylab = "NOx (micrograms/J)",
2010 Jul 12
2
Xyplot or Tin-R problem?
I ran the following script from xyplot Examples using Tin-R on
Windows and saw no plot produced.
EE <- equal.count(ethanol$E, number=9, overlap=1/4)
xyplot(NOx ~ C | EE, data=ethanol,
prepanel = function(x,y) prepanel.loess(x, y, span=1),
xlab="Compression Ratio", ylab="NOx (micrograms/J)",
panel = function(x,y) {
panel.grid()(h = -1, v=2)
2001 Nov 28
3
trellis plot
Hi,
I'd like to plot 4 groups of data using xyplot and panel.superpose so
that the points are overlayed on a single plot. For each group of data
I'd also like a loess smoothed function (using panel.loess). I have
tried the following:
xyplot(series ~ time | gr, data=etable,
panel = function(x,y, ...) {
panel.superpose(x,y, ...)
panel.loess(x,y,span=.15)
2008 Nov 19
1
Buggy trellis.focus() with xyplot ?
Hi:
(Tried to find a bug report about this issue, but was unable to find it, let
me know if this is a known issue)
I have been working on an interface to highlight xyplot panels on mouse
overs in JavaGD but I have stumbled with what seems to be a bug in
trellis.focus.
I am using R 2.8 with lattice 0.17-15
*** To replicate the bug:
1.- display an xyplot. For example, from the xyplot help page:
2010 Apr 09
2
lattice background grid
Hi,
I wonder if anyone knows how to control color and darkness of the
background grid generated by ' type = c("g", "p") ' in a lattice plot
(e.g., in xyplot). The documentation does not seem to offer a way to
change them.
Thanks a lot!
Wen
2006 Mar 29
2
bivariate case in Local Polynomials regression
Hi:
I am using the package "KernSmooth" to do the local polynomial regression. However, it seems the function "locpoly" can only deal with univariate covaraite. I wonder is there any kernel smoothing package in R can deal with bivariate covariates? I also checked the package "lcofit" in which function "lcofit" can indeed deal with bivariate case. The
2010 Dec 02
4
Integral of PDF
The integral of any probability density from -Inf to Inf should equal 1, correct? I don't understand last result below.
> integrate(function(x) dnorm(x, 0,1), -Inf, Inf)
1 with absolute error < 9.4e-05
> integrate(function(x) dnorm(x, 100,10), -Inf, Inf)
1 with absolute error < 0.00012
> integrate(function(x) dnorm(x, 500,50), -Inf, Inf)
8.410947e-11 with absolute error <
2010 Mar 19
3
strange behavior, adds new field by non-existent field
data:
> tmp1
Date HrMn Temp Q.4
1 19450101 0900 -37.0 1
2 19450101 1000 -35.9 2
3 19450101 1100 -35.9 3
4 19450101 1200 -36.4 4
5 19450101 1300 -36.4 5
6 19450101 1400 -36.4 6
7 19450101 1500 -36.4 7
8 19450101 1600 -37.5 9
Accidentally, I did this (I meant to write Q.4 instead of Q here)
> tmp1$Q[tmp1$Q!="1" & tmp1$Q!= "5"]<-NA
I would
2009 Dec 17
2
Wrong links in help help files? (PR#14155)
Dear R developers,
?help.search
...opens a tab in my browser (Firefox 3.5.6, Win XP):
http://127.0.0.1:31642/library/utils/html/help.search.html
When I click on the link for "?",
(http://127.0.0.1:31642/library/utils/help/?)
in the last line of "Details", I get an error message:
"Server error: invalid response from R".
Similarly,
The link to "??"
2010 Jan 11
2
Documentation: format of read.table help text (PR#14180)
Dear R developers,
I find the format of descriptions of the arguments in the read.table
help text slightly inconsistent.
For example, the logical arguments comes in seven different formats,
more or less explicit about the consequences of a TRUE (or FALSE):
1. check.names logical. If TRUE then the names...
2. blank.lines.skip logical: if TRUE blank lines...
3. flush logical: if TRUE,
2010 Sep 03
1
Help on Select.list
Hi R,
I am using select.list
names=c("Ravi", "Raj","Shubha","Nivriti")
select.list(names) provides a drop down to choose one of the 4 names.
However I would like to know if it is possible to create a
classification something like this
select.list(names) should give
Boys
Ravi
Raj
Girls
Shubha
Nivriti
I should be able to choose
2010 Jun 29
1
Use of processor by R 32bit on a 64bit machine
Dear all,
I've recently purchased a new 64bit system with an intel i7 quadcore
processor. As I understood (maybe wrongly) that to date the 32bit
version of R is more stable than the 64bit, I installed the 32bit
version and am happily using it ever since. Now I'm running a whole
lot of models, which goes smoothly, and I thought out of curiosity to
check how much processor I'm using. I
2010 Feb 10
1
R crashes when setWinProgressBar is give a numeric value for label argument
This problem can be seen by the following commands:
> pb <- winProgressBar(max=1000, label='0')
> b <- 1
> setWinProgressBar(pb, b, label=b)
This set of commands (on windows of course, XP in this case) causes R to crash.
This is not strictly a bug since the documentation states that the label argument should be a character string and using as.character(b) does work
2010 Mar 30
1
predict.kohonen for SOM returns NA?
All,
The kohonen predict function is returning NA for SOM predictions
regardless of data used... even the package example for a SOM using
wine data is returning NA's
Does anyone have a working example SOM. Also, what is the purpose of
trainY, what would be the dependent data for an unsupervised SOM?
As may be apparent to you by my questions, I am very new to kohonen
maps and am very grateful
2010 May 26
2
extracat , JGR, iWidgets install problems
[Environment: Win XP, R 2.10.1]
I'm trying to install the packages JGR and iWidgets required by the
extracat package to make the interactive plots
in the package work. I've tried various things, but nothing seems to
work. Here is my most recent attempt,
followed by my sessionInfo().
Does anyone have any suggestions how to make this work?
>
> library(extracat)
Loading
2010 Mar 18
1
probable timezone confusion with as.yearmon
It looks like a timezone issue, and it's causing confusion to me at least.
My original data:
gmt <-
c("19880101 0000", "19880101 0100", "19880101 0300", "19880101 0400",
"19880101 0500", "19880101 0600")
These were converted to local dates/times with
akst<-strptime(gmt,format="%Y%m%d %H%M")-(3600*9) # because I want
2002 Dec 06
3
smooth curves
I would like to draw smooth curves instead of polygons.
I could not find any spline function to do that : given an x and a y vectors, they all take the x in increasing order.
Is there a function to draw a smooth curve through a set of points in any order ?
with many thanks in advance
Alain Guerreau CNRS Paris
[[alternate HTML version deleted]]
2008 Jul 08
1
R crash with ATLAS precompiled Rblas.dll on Windows XP Core2 Duo
I noticed a problem using R 2.7.1 on Windows XP SP2 with the precompiled
Atlas Rblas.dll. Running the code below causes R to crash. I started R
using Rgui --vanilla and am using the precompiled Atlas Rblas.dll from
cran.fhcrc.org dated 17-Jul-2007 05:04 for Core2 Duo.
The code that causes the crash:
x <- rnorm(100)
y <- rnorm(100)
z <- rnorm(100)
loess(z ~ x * y)
loess(z ~ x) does
2010 Apr 22
2
Compare two data frames
I wonder if there is a more efficient way to do this task. Suppose I have two data frames, such as
d1 <- data.frame(x = c(1,2,3), y = c(4,5,6), z = c(7,8,9))
d2 <- d1[, c('y', 'x')]
The first dataframe d1 has more variables than d2 and the variable columns are in a different order.
So, what I want to do is compare the two frames on the variables that are common between
2010 Mar 18
1
How to read.table with “Hebrew” column names (in R)?
(I am reposting this question after a few months without a solution...)
Hi all,
I am trying to read a .txt file, with Hebrew column names, but without
success.
I uploaded an example file to: http://www.talgalili.com/files/aa.txt
And tried the command:
read.table("http://www.talgalili.com/files/aa.txt", header = T, sep = "\t")
This returns me with:
X.....ª X...ª......