Hi
the emulator package of the the BACCO bundle
includes kriging as a special case.
HTH
Robin
On 4 Jan 2007, at 04:56, [EMAIL PROTECTED] wrote:
Dear R-list members,
I wish everyone a happy and successful 2007!
Does anyone know of R-based software for
optimal spatial prediction
I think I have figured out part of the answer, when I entered a debugger, the
calling environment (parent frame)
refers to the debugger program, instead of the mainfun( ). But is there
anyway to solve this problem ?
thanks
- Original Message -
From: Tong Wang [EMAIL PROTECTED]
On Thu, 4 Jan 2007, [EMAIL PROTECTED] wrote:
Remember that -2 * the difference in the likelihoods between the two
models is asymptotically chi-squared distributed, with degrees of
freedom equal to the difference in number of parameters between the
models. So you can just calculate that for
Dear helpers,
I have a question about the SVG device. It works fine, the SVG file is
indeed produced, only the graphic differs from the R window.
In the SVG file the dashed line is just a regular plain one. My toy example is:
library(RSvgDevice)
devSVG(myplot.svg, width=10, height=10)
plot(1:10)
Check out #2 in:
http://finzi.psych.upenn.edu/R/Rhelp02a/archive/85801.html
and RSiteSearch(axis(4) to find additional examples.
On 1/4/07, Arun Kumar Saha [EMAIL PROTECTED] wrote:
Dear Gabor,
Thank you very much for your letter. Actually I got partial solution from
your suggestion. Still I
Sorry for duplicating the message, the previous had an unintended
subject line...
Dear helpers,
I have a question about the SVG device. It works fine, the SVG file is
indeed produced, only the graphic differs from the R window.
In the SVG file the dashed line is just a regular plain one. My toy
Dear list members!
I've two questions concerning graphic export:
a) I want to export my graphics as PostScript files. in this way I use the
postscript() function. The tricky part is that they must have a pretended size
(7 x 7 cm) and an absoulte font size (10pt).
b) how can i (permanent)
Hi
On 4 Jan 2007 at 14:18, Arun Kumar Saha wrote:
Date sent: Thu, 4 Jan 2007 14:18:11 +0530
From: Arun Kumar Saha [EMAIL PROTECTED]
To: Gabor Grothendieck [EMAIL PROTECTED]
Copies to: r-help@stat.math.ethz.ch
On Thu, 2007-01-04 at 12:21 +0100, mirca heli wrote:
Dear list members!
I've two questions concerning graphic export:
a) I want to export my graphics as PostScript files. in this way I use
the postscript() function. The tricky part is that they must have a
pretended size (7 x 7 cm) and an
mirca heli wrote:
Dear list members!
I've two questions concerning graphic export:
a) I want to export my graphics as PostScript files. in this way I use the
postscript() function. The tricky part is that they must have a pretended
size (7 x 7 cm) and an absoulte font size (10pt).
On 1/4/2007 6:21 AM, mirca heli wrote:
Dear list members!
I've two questions concerning graphic export:
a) I want to export my graphics as PostScript files. in this way I use the
postscript() function. The tricky part is that they must have a pretended
size (7 x 7 cm) and an absoulte
Hello, and Happy New Year. My default working directory is getting very
cluttered. I know that I should be using a different working directory for
each project (I work in Windows), but do not know how to go about creating
different ones and moving back and forth between them. I have read
See ?getwd and ?setwd to set the working directory
See ?load and ?save to read the workspace.
Cheers,
Thierry
ir. Thierry Onkelinx
Instituut voor natuur- en bosonderzoek / Reseach Institute for Nature
and Forest
The basic command for this is setwd() (for set working directory)
You can save your workspace after you use the setwd command, and this
will create a .RData file in that directory. In windows you can click on
the .RData directory to open R using that particular working directory.
You can
Bill Shipley wrote:
Hello, and Happy New Year. My default working directory is getting very
cluttered. I know that I should be using a different working directory for
each project (I work in Windows), but do not know how to go about creating
different ones and moving back and forth between
use setwd() and getwd(). see ?setwd
or on the shortcut properties tab, set the target directory to your working dir
and have
as many shortcuts as working directories.
working with R commands does it for me.
all the best.
A.
- Original Message
From: Bill Shipley [EMAIL PROTECTED]
To: R
An additional note for Windows. The directory name needs to be written
as C:/Project/Working directory or C:\\Project\\Working directory as
opposed to the usual way of referencing directories in Windows.
AA wrote:
use setwd() and getwd(). see ?setwd
or on the shortcut properties tab, set the
Abhijit Dasgupta wrote:
An additional note for Windows. The directory name needs to be written
as C:/Project/Working directory or C:\\Project\\Working directory as
opposed to the usual way of referencing directories in Windows.
Which are? Those are the forms given in the Microsoft
Bill Shipley wrote:
Hello, and Happy New Year. My default working directory is getting very
cluttered. I know that I should be using a different working directory for
each project (I work in Windows), but do not know how to go about creating
different ones and moving back and forth between
On 1/3/07, Martin Henry H. Stevens [EMAIL PROTECTED] wrote:
Hi folks,
I have assumed that ratios of variance components (Fst and Qst in
population genetics) could be estimated using the output of mcmcsamp
(the series on mcmc sample estimates of variance components).
What I have started to do
Dear useRs,
I have a few hundred plots that I'd like to export to one document.
pdf() isn't an option, because the file created is prohibitively huge
(due to scatter plots with many points). So I have to use png()
instead, but then I end up with a lot of files (would prefer just
one).
1. Is
Odd behavior from get.hist.quote this AM.
get.hist.quote('sunw')
trying URL
'http://chart.yahoo.com/table.csv?s=sunwa=0b=02c=1991d=0e=03f=2007g=dq=qy=0z=sunwx=.csv'
Content type 'text/csv' length unknown
opened URL
.. .. .. .. ..
.. ..
I can think of two options:
1. Use R2HTML and save the html output as PDF
2. Use Sweave and compile the LaTeX file to PDF. Search the mailing list
archive on how to save the graphs as png or jpeg (as Sweave will
standard generate eps or pdf graphs).
Cheers,
Thierry
Here's an example illustrating a way to get a second y axis that has
a different range:
x - 1:10
y1 - 2*x
y2 - 100-3*x+rnorm(10)
par(mar=c(5.1,4.1,4.1,4.1))
plot(x,y1)
par(new=TRUE)
plot(x,y2,xaxt='n',yaxt='n',xlab='',ylab='',pch=3)
axis(4)
mtext('y2',side=4,line=2.5)
-Don
At 2:18 PM +0530
On 1/4/07, BBands [EMAIL PROTECTED] wrote:
Odd behavior from get.hist.quote this AM.
Now working, must have been a Yahoo! issue.
jab
--
John Bollinger, CFA, CMT
www.BollingerBands.com
If you advance far enough, you arrive at the beginning.
__
Bill,
I like to use Windows Explorer to find folders and then launch R with the
selected folder as the working directory.
I put some notes online about this:
http://research.stowers-institute.org/efg/R/TechNote/WindowsExplorerWorkingDirectory/index.htm
efg
Bill Shipley [EMAIL PROTECTED] wrote
Hello R List -
I have to import Excel files (either as .csv files or using RODBC) into
R (2.4.1, Windows) and operate on dates and times (e.g. find minutes
between times, change dates to days of week or analyze by weeks of
year). The help files for format.Date, strptime, as.POSIX,
On 1/4/2007 9:41 AM, Bill Shipley wrote:
Hello, and Happy New Year. My default working directory is getting very
cluttered. I know that I should be using a different working directory for
each project (I work in Windows), but do not know how to go about creating
different ones and moving
I think the best explanation of dates and times is in r-news 2.4.1 but
2.4.1 might be off so someone will hopefully correct me if I'm wrong.
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of Ben Fairbank
Sent: Thursday, January 04, 2007 12:43 PM
To:
On 1/4/2007 12:42 PM, Ben Fairbank wrote:
Hello R List -
I have to import Excel files (either as .csv files or using RODBC) into
R (2.4.1, Windows) and operate on dates and times (e.g. find minutes
between times, change dates to days of week or analyze by weeks of
year). The help
The approach I would take (possibly due to ignorance of a better option)
is to export to the multiple .png files, then use a tool like
imagemagick to combine them into a single pdf file.
For a quick test I exported 3 graphs from R and called them test1.png,
test2.png, and test3.png. The
On Jan 4, 2007, at 11:18 AM, Douglas Bates wrote:
On 1/3/07, Martin Henry H. Stevens [EMAIL PROTECTED] wrote:
Hi folks,
I have assumed that ratios of variance components (Fst and Qst in
population genetics) could be estimated using the output of mcmcsamp
(the series on mcmc sample estimates
All,
I am glad all of you have benefited from the posting of the RMySQL 5-10
zip file on my university website. I am asking for some help from the
group, I am leaving the university at the end of the semester and I need
a place to post this file until I get settled in my new position.
Anyone
On 1/4/2007 12:48 PM, Duncan Murdoch wrote:
On 1/4/2007 9:41 AM, Bill Shipley wrote:
Hello, and Happy New Year. My default working directory is getting very
cluttered. I know that I should be using a different working directory for
each project (I work in Windows), but do not know how to go
Hi, I'm trying to write a series of pipes using littler, and I get the
following behaviour: Sorry if I'm just doing something witless, I'm new to
R. I'm using the latest versions from debian testing (2.4.0 and 0.0.8).
$ r -e 'a-dget(file=stdin()); print(a)'
?list(a=2)
Segmentation fault
In R
Hi,
I have tried
for (i in 1:100) L[,i] - loess((i = =(1:100))~I(1:100), span=.5,
degree=1)$fit
to create a matrix which gives me the smoothing weights (correctly as
far as I have experienced), eg.
yhat - loess(y~I(1:100), span=.5,degree=1)$fit
yhat[30]
[1] -0.2131983
L[30,]%*%y
Dear R-wizards,
I have a population from which I want to draw a stratified sample by region.
In Venables and Ripley Modern Applied statistics with S I found some
great procedures for Simple Random Sampling (with and without replacement)
and for Systematic sampling and it works!
For
I think the question is discussed in other thread, but I don't exactly find
what I want .
I'm working in Windows XP with 2GB of memory and a Pentium 4 - 3.00Ghx.
I have the necessity of working with large dataset, generally from 300,000
records to 800,000 (according to the project), and about 300
Does anyone know a reason why, in principle, a call to randomForest
cannot accept a data frame with missing predictor values? If each
individual tree is built using CART, then it seems like this
should be possible. (I understand that one may impute missing values
using rfImpute or some other
A new version of RSQLite has been pushed to CRAN.
In this version...
* Further integration of the manifest type system available since
SQLite 3. We now obtain the column type from the DB instead of
pulling everything across as a character vector and calling
type.convert. This should
Hi All,
Weighting is the procedure to correct the distributions in the sample data to
approximate those of the population from which it is drawn. This is partly a
matter of expansion and partly a matter of correction or adjustment for both
non response and non coverage. It serves the
Hy all,
I'm plotting graphs using plot() function, they are on X axes POSIX dates:
POSIXt oldClass POSIXct POSIXlt
I can't figure out why sometimes it prints the month and days and sometimes it
prints the unix timestamp.
It appens usually when the xlim is short like only some days.
xlim is
I have a set of timestamp data that I have in a text file that I would like
to import into R for analysis.
The timestamps are formated as follows:
DT_1,DT_2
[2006/08/10 21:12:14 ],[2006/08/10 21:54:00 ]
[2006/08/10 20:42:00 ],[2006/08/10 22:48:00 ]
[2006/08/10 20:58:00 ],[2006/08/10 21:39:00 ]
I don't know about this module, but a general answer is that if you have
missing data, it may affect your model. If your data is missing at
random, then you might be lucky in your model building.
If however your data was not missing at random (e.g. censoring) , you
might build a wrong predictor.
Folks:
Motivated by the recent thread on setting working directories, below are a
couple of functions for GUI-izing saving and loading files **in Windows
only** that sort of takes care of this automatically. The simple strategy
is just to maintain a file consisting of the filenames of recently
Yes I completely agree with your statements. As far as a way around
it, I would say that CART has some facilities for dealing with
missing data. e.g. when an observation is dropped into the tree and
encounters a split at which the variable is missing, then one option
is to simply not send it
Try this:
Lines - DT_1,DT_2
[2006/08/10 21:12:14 ],[2006/08/10 21:54:00 ]
[2006/08/10 20:42:00 ],[2006/08/10 22:48:00 ]
[2006/08/10 20:58:00 ],[2006/08/10 21:39:00 ]
[2006/08/04 12:15:24 ],[2006/08/04 12:20:00 ]
[2006/08/04 12:02:00 ],[2006/08/04 14:20:00 ]
Lines2 - gsub(\\[|\\], ,
Hello,
thanks for the advice. I was aware that the version was out of date and
your message prompted me to finally upgrade to the latest version 2.4.1.
Unfortunately, running R under Emacs ESS, the problem I described earlier
persists.
I have also tried the Save as option under the standard
Hello, everyone, I'm using R to deal with a REML problem. I found lmer is
the right function for this. But I got stuck because I couldn't interpret
the result. I'm attaching a short example of my executing log. Please have a
look and give me some advice on it. Thanks a lot!
Plot Block
John Lawrence Aspden wrote:
Hi, I'm trying to write a series of pipes using littler, and I get the
following behaviour: Sorry if I'm just doing something witless, I'm new to
R. I'm using the latest versions from debian testing (2.4.0 and 0.0.8).
$ r -e 'a-dget(file=stdin()); print(a)'
You can try randomForest in Fortran codes, which has that function
doing missing replacement automatically. There are two ways of
imputations (one is fast and the other is time-consuming) to do that.
I did it long time ago.
the link is below. If you have any question, just let me know.
The good news, you don't have to shut down R. Several control-G in
the *R* buffer in emacs will recover control. A second attempt in the
same GUI graphics device did get the postscript file saved.
A workaround for this problem is to use the command line, rather than
the GUI menu, to save the
Please read the rw-FAQ Q2.9. There are ways to raise the limit, and you
have not told us that you used them (nor the version of R you used, which
matters as the limits are version-specific).
Beyond that, there are ways to use read.table more efficiently: see its
help page and the 'R Data
As the footer says
PLEASE do read the posting guide
http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
We can't help you with a problem we cannot reproduce, and random guessing
is not going to be productive.
On Thu, 4 Jan
Dear R-helpers,
I'm plotting geophysical data in the form of contours using
filled.contour. The display would be much more effective if the areas
with negative values could be color coded
by -- say -- cold colors in the range of blue to green, and conversely
the areas with positive values got
On Thu, 4 Jan 2007, Jukka Nyblom wrote:
Hi,
I have tried
for (i in 1:100) L[,i] - loess((i = =(1:100))~I(1:100), span=.5,
degree=1)$fit
to create a matrix which gives me the smoothing weights (correctly as
far as I have experienced), eg.
yhat - loess(y~I(1:100), span=.5,degree=1)$fit
Jeffrey Horner wrote:
John Lawrence Aspden wrote:
Hi, I'm trying to write a series of pipes using littler, and I get the
following behaviour:
$ r -e 'a-dget(file=stdin()); print(a)'
?list(a=2)
Segmentation fault
You've found a bug which has been fixed. Expect a new version 0.0.9 of
Hi Claudia,
It's quite easy to do this using ggplot, although you get exactly the
same appearance as filled.contour (hopefully in the next version).
Have a look at ggtile and scgradient.
Regards,
Hadley
On 1/4/07, Claudia Tebaldi [EMAIL PROTECTED] wrote:
Dear R-helpers,
I'm plotting
On 1/4/2007 5:33 PM, [EMAIL PROTECTED] wrote:
Hello,
thanks for the advice. I was aware that the version was out of date and
your message prompted me to finally upgrade to the latest version 2.4.1.
Unfortunately, running R under Emacs ESS, the problem I described earlier
persists.
I
Just in case you were not already aware, you could try Tinn-R
(http://www.sciviews.org/Tinn-R/) or JGR (http://rosuda.org/JGR/)
or one of the editors described at this site
(http://www.sciviews.org/_rgui/projects/Editors.html) as Emacs
alternatives. I love Emacs and miss the speed of using Emacs
I want to extract estimated coeffiicents of each local polynomial at given x
from loess(), locfit(), or KernSmooth(). Can some experts provide me with
suggestions? Thanks.
Delong Liu
__
R-help@stat.math.ethz.ch mailing list
[Using R 2.2.0 on Windows XP; OK, OK, I will update soon!]
I have noticed some undesirable behaviour when applying
ifelse to a data frame. Here is my code:
A - scan()
1.00 0.00 0.00 0 0.0
0.027702 0.972045 0.000253 0 0.0
A - matrix(A,nrow=2,ncol=5,byrow=T)
A == 0
Can you send ess-bugs some reproducible examples of what isn't working?
These are things that are probably easily fixed.
Rich
Original message
Date: Thu, 04 Jan 2007 21:02:17 -0500
From: [EMAIL PROTECTED]
Subject: Re: [R] R grahics: Save as hangs computer
To: [EMAIL PROTECTED],
Get the RColorBrewer package from CRAN
Description: The packages provides palettes for drawing nice maps
shaded according to a variable.
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the
Hello list
thanks for that advice Rich. You are right a few control-G allowed me to
recover and it had saved an eps file.
The dev.copy2eps appears to take exactly what is visible in the graphics
window and convert it to an eps file. This does the job for the moment.
Opening the eps file in
R is actually smarter than that. The bounding box is square.
From
plot(1:10)
dev.copy2eps(file=tmp.eps)
The bounding box (which you can see by opening tmp.eps in emacs) is
%%BoundingBox: 0 0 517 517
and LaTeX will therefore interpret it correctly. GSview is misleading
you since it puts that
An additional note.
You can see the bounding box in the GSview display by clicking
Options/Show Bounding Box
Rich
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide
another note
you can clip the GSview to just the bounding box with
Options/EPS CCip
Rich
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
It can be explained.
class(A)
[1] data.frame
length(A)
[1] 5
class(A==0)
[1] matrix
length(A==0)
[1] 10
class(-A*log(A))
[1] data.frame
length(-A*log(A))
[1] 5
as you can see, the result of A==0 is matrix with length=10, while the
result of -A*log(A) is still data.frame with length=5.
Suppose that you have the timestamps in a comma seperated text file,
times.dat then do
tid - read.csv(times.dat, header = TRUE, colClasses = character)
# Note the colClasses argument! If not used the columns of tid are factors by
default.
time.conv - function(x) as.POSIXct(strptime(x, format
70 matches
Mail list logo