You might wish to have a look at the 'its' package for irregular
time-series on CRAN. If your prices are in an its called price,
then the following will get you on your way. Since it is not efficient
either in storage or computation, I offer it because it might be
convenient for display,
I am trying to fit a rank-frequency distribution with 3 unknowns (a, b
and k) to a set of data.
This is my data set:
y - c(37047647,27083970,23944887,22536157,20133224,
20088720,18774883,18415648,17103717,13580739,12350767,
8682289,7496355,7248810,7022120,6396495,6262477,6005496,
Is there any way to change the y axis range of values in a plot.gam()? I
need that two different GAM plots to be of the same scale.
- Sorry, I don't think there is a simple way of doing it, except by
editing plot.gam to fix `ylim' by hand It's an obvious thing to want
to do and I'll fix it
Dr Andrew Wilson [EMAIL PROTECTED] writes:
I am trying to fit a rank-frequency distribution with 3 unknowns (a, b
and k) to a set of data.
This is my data set:
y - c(37047647,27083970,23944887,22536157,20133224,
20088720,18774883,18415648,17103717,13580739,12350767,
Martin Maechler has pointed out to me that I omitted to include my values
for x when circulating my query. They are a simple rank list from 1 to
26:
x - 1:26
Thanks,
Andrew
__
[EMAIL PROTECTED] mailing list
Thanks for this clarification.
I have learned in the meantime that it is necessary to be very careful when
using all these POSIX things.
As another example, here is something that made me scratch my head just
yesterday:
When I create a sequence of days that happens to start before and ends in
Pascal A. Niklaus wrote:
Is there an shortcut (compared to manually swap row vectors) for the
random permutation of rows of a matrix?
sample does not act on the element row vectors as a whole.
True, but sample(nrow(foo)) gives you a permutation of 1:nrow(foo)
that you can subscript with:
Hi,
i getting following error and don't know doing
something wrong.
mysqlWriteTable(con,model1,model1,overwrite=T)
Error in [.data.frame(value, from:to, drop = FALSE) :
undefined columns selected
In addition: Warning message:
drop argument will be ignored in: [.data.frame(value, from:to,
Hi,
is it correct that i need ~ 2GB RAM that it's
possible to work with the default setting
ntree=500 and a data.frame with 100.000 rows
and max. 10 columns for training and testing?
P.S.
It's possible calculate approximate the
memory demand for different settings with RF?
Many thanks
Hello List
does anyone have an R function for the Lambert W function? I need
complex arguments.
[the Lamert W function W(z) satisfies
W(z)*exp(W(z)) = z
but I could'nt even figure out how to use uniroot() for complex z]
--
Robin Hankin
Uncertainty Analyst
Southampton Oceanography Centre
Hi,
is it correct that i need ~ 2GB RAM that it's
possible to work with the default setting
ntree=500 and a data.frame with 100.000 rows
and max. 10 columns for training and testing?
no. You may parallelize the computations: perform 5 runs of RF with `ntree
= 100' (or less) and save the
Hello R-Users!
Perhaps you can help me!
I am a newbie at R and I am searching for a propability to realize the
Coxian-Distribution.
http://www.owlnet.rice.edu/~elec428/handouts/Coxian.pdf
Thanks in advance
Tobias
[[alternative HTML version deleted]]
Updating my 1.8.0 R installation (update.packages() ) I obtain the following
(SORRY FOR THE LENGTH OF THE LOG BUT IT HELPS!!!):
downloaded 135Kb
KernSmooth :
Version 2.22-11 in /usr/lib/R/library
Version 2.22-12 on CRAN
Update (y/N)? y
mgcv :
Version 0.9-3.1 in
On Tue, 25 Nov 2003, [EMAIL PROTECTED] wrote:
I think you need to tell is specifically which platform you are using
(output of version), and whether you installed R from source or a binary
distribution. Most likely you are on a Unix or Linux platform, and
installed a binary distribution, but
Hi
I am using R as a back-end to some CGI scripts, written in Perl. My platform is Suse
Linux 8.2, Apache 1.3.7. So the CGI script takes some form parameters, opens a pipe
to an R process, loads up some Bioconductor libraries, executes some R commands and
takes the ouput and creates a web
Hi, there!
I finally became a disciple of 'R', after having lost years of my life
handling data with a popular, rather wide-spread spreadsheet-software.
Now I want to plot the results of many runs of my simulation over time,
so that the means +/- Standard error are on the y-axis, and time on
From: Christian Schulz
Hi,
is it correct that i need ~ 2GB RAM that it's
possible to work with the default setting
ntree=500 and a data.frame with 100.000 rows
and max. 10 columns for training and testing?
If you have the test set, and don't need the forest for predicting other
data,
Is the contributed package `deldir' on CRAN what you're looking for?
HTH,
andy
From: Pascal A. Niklaus [mailto:[EMAIL PROTECTED]
Hi all,
Given a number of points (x,y) in a plane, I'd like to plot a map of
polygons, so that
1) each polygon contains exactly one point
2) the
or tripack perhaps...
url:www.econ.uiuc.edu/~roger/my.htmlRoger Koenker
email [EMAIL PROTECTED] Department of Economics
vox:217-333-4558University of Illinois
fax:217-244-6678Champaign, IL 61820
Luis,
Thanks for your thoughtful comments. Indeed you've uncovered
a bug/problem in that there's no way for users to control the
allow.keywords= argument in calls to dbWriteTable() -- this
needs to be fixed. Regarding the default value for allow.keywords,
I'm not sure it is wise to set it to
With a custom compiled kernel, I've run R processes that used more than 5GB
of RAM on a Linux box with 8GB RAM and dual Xeons. So it seems to work on
32-bit Linux with big memory kernel.
Andy
From: Duncan Murdoch [mailto:[EMAIL PROTECTED]
[snip]
Normally the maximum memory allowed for any
Starting up R and loading libraries can be very time consuming. For my
RSOAP system (http://www.analytics.washington.edu/Zope/projects/RSOAP/) I
took the step of pre-starting the R process, including the loading of some
libraries, and then handing work of to the pre-started process. You should
See the function plotCI() in package gregmisc on CRAN.
EP
At 15:03 25/11/2003 +0100, vous avez écrit:
Hi, there!
I finally became a disciple of 'R', after having lost years of my life
handling data with a popular, rather wide-spread spreadsheet-software.
Now I want to plot the results of many
Hi Angel and * !
I also use kate and in my opinion it is very convenient. I don't know which
features you don't like, but if you mean some syntax highlighting features,
you might try my R syntax highlighting (XML) file for kate. It is based on
Egon Willighagen's one and contains a few bug
Since x - 1:26 and your y's are all positive, your model,
ignoring the error term, is mathematically isomorphic to the following:
x - 1:26
(fit - lm(y~x+log(x)))
Call:
lm(formula = y ~ x + log(x))
Coefficients:
(Intercept)x log(x)
35802074 -371008 -8222922
If the numbers are letter frequencies, I would suggest Poisson
regression using glm the default link is logarithms, and that should
work quite well for you.
hope this helps. spencer graves
###
Very many thanks for your help.
What do these numbers
Hello. I want to specify a diagonal structure for the covariance matrix
of random effects in the lme() function.
Here is the call before I specify a diagonal structure:
fit2-lme(Ln.rgr~I(Ln.nar-log(0.0011)),data=meta.analysis,
+ random=~1+I(Ln.nar-log(0.0011)|STUDY.CODE,na.action=na.omit)
Pascal == Pascal A Niklaus [EMAIL PROTECTED]
on Tue, 25 Nov 2003 16:10:56 +0100 writes:
Pascal Hi all, I am playing around with latin squares, and
Pascal wrote a recursive function that searches for valid
Pascal combinations. Apart from the fact that there are
Pascal very
On Sun, Nov 23, 2003 at 11:38:05AM +, Florian Roedhammer wrote:
I would need information concerning the features of R such as
static/dynamic typing, method overlapping, object oriontation,
My question is whether you know any site where such features are
discussed on or
On 25 Nov 2003 15:03:49 +0100
Jan Wantia [EMAIL PROTECTED] wrote:
Hi, there!
I finally became a disciple of 'R', after having lost years of my life
handling data with a popular, rather wide-spread spreadsheet-software.
Now I want to plot the results of many runs of my simulation over
Your starting values for the parameters are no even in the general
ballpark. Here is what I got for the final fit:
Parameters:
Estimate Std. Error t value Pr(|t|)
a 3.806e+07 1.732e+06 21.971 2e-16 ***
k -3.391e-02 6.903e-02 -0.4910.628
b 9.000e-01 1.240e-02 72.612 2e-16 ***
On Tue, 25 Nov 2003, michael watson (IAH-C) wrote:
Hi
I am using R as a back-end to some CGI scripts, written in Perl. My platform is
Suse Linux 8.2, Apache 1.3.7. So the CGI script takes some form parameters, opens a
pipe to an R process, loads up some Bioconductor libraries, executes
Bill Shipley [EMAIL PROTECTED] writes:
Hello. I want to specify a diagonal structure for the covariance matrix
of random effects in the lme() function.
Here is the call before I specify a diagonal structure:
fit2-lme(Ln.rgr~I(Ln.nar-log(0.0011)),data=meta.analysis,
+
It looks like image_and_label has only 2 columns, so when you take
img_and_label[,2] you have a vector left. Even if that weren't the case,
you're going to need to pass in both the gray scale points and labels,
presumably in a data frame. You've created a character matrix below, so
you're just
It's not clear to me what you want to do, but if I understand your problem
somewhat, I don't see how randomForest would be relevant.
Sounds like you are doing the following:
o Read in a 512x512 image with pixel intensities.
o You somehow fit a 3-component normal mixture model to the intensity
How to control line width ?
if I do:
postscript(IC50-density.eps, width = 4.0, height = 3.0, horizontal =
FALSE, onefile = FALSE, paper = special, title = IC50 distribution)
plot(d$x, d$y, xlab = -log10(IC50), ylab = density)
lines(d$x, d$y, lwd = 0.1)
dev.off()
but whatever value I give
Sorry. I need to retract my claim. There seems to be a 3G limit, even
though the OS could handle nearly 8G. (I can have two simultaneous R
processes each using near 3G.)
On another note, on our dual Opteron box R (compiled as 64-bit) could easily
use nearly all the 16G in that box (that's one
When I use
plot(..., type = line) then ldw parameter makes a difference...
Because of large number of points overlapping I simply han an impression
before that I am getting thick line...
R
Ryszard Czerminski/PH/[EMAIL PROTECTED]
Sent by: [EMAIL PROTECTED]
11/25/2003 01:34 PM
To:
Dear MacR users
During the weekend I did a clean install of MacOSX 10.3.1, of Apples X11,
of Apples development tools, and of the basic Fink 0.6.2 package.
Now I have just downloaded Raqua.dmg from CRAN and installed the
RAqua, libreadline and tcltk packages.
Sorrowly, double-clicking on StartR
From: Douglas Bates [mailto:[EMAIL PROTECTED]
Liaw, Andy [EMAIL PROTECTED] writes:
Sorry. I need to retract my claim. There seems to be a 3G
limit, even
though the OS could handle nearly 8G. (I can have two
simultaneous R
processes each using near 3G.)
On another note, on
Dear MacR users
Sorry for the earlier mail, now everything works really great.
For some reason I could not start R immediately after installation,
I had to log out first and then login again, then R did start from
a really great R Console.
Best regards
Christian
cstrato wrote:
Dear MacR users
michael watson (IAH-C) wrote:
I am trying to make my cgi scripts quicker and it turns out that the
bottle-neck is the loading of the libraries into R - for example
loading up marrayPlots into R takes 10-20 seconds, which although not
long, is long enough for users to imagine it is not working and
Hi,
I have some difficulties to figure how to set a range to my histogram bands.
I have values that are [0,50] and they appear once in my list. How
can I do a histogram that include all the values between a range of
1 together? [0,1],[10001,2],[21,3], ...
Thanks,
Mathieu
How do I go about generating a WEIGHTED mean (and standard error) of a
variable (e.g., expenditures) for each level of a categorical variable
(e.g., geographic region)? I'm looking for something comparable to PROC
MEANS in SAS with both a class and weight statement.
Thanks.
Marc
Mathieu Drapeau wrote:
Hi,
I have some difficulties to figure how to set a range to my histogram
bands.
I have values that are [0,50] and they appear once in my list. How
can I do a histogram that include all the values between a range of
1 together?
check out plotCI and plotmeans in the gregmisc library.
Duncan
*
Dr. Duncan Mackay
School of Biological Sciences
Flinders University
GPO Box 2100
Adelaide
S.A.5001
AUSTRALIA
Ph (08) 8201 2627FAX (08) 8201 3015
Christian Reilly [EMAIL PROTECTED] writes:
I'm having trouble plotting a curve (basically a dose-response
function) through a set of data.
I have created a dataframe (df) of Stimulus Intensities (xstim) and
Normalized Responses (yresp), and I've used nls() to calculate a nonlinear
[EMAIL PROTECTED] wrote:
How do I go about generating a WEIGHTED mean (and standard error) of a
variable (e.g., expenditures) for each level of a categorical variable
(e.g., geographic region)? I'm looking for something comparable to PROC
MEANS in SAS with both a class and weight statement.
Dear Marc,
For the weighted mean, one possible solution is as follows
and will hopefully give you the general idea:
tmp - data.frame(x=sample(1:5, 100, replace=TRUE),
y=sample(1:100, 100, replace=TRUE),
w=runif(100))
lapply(split(tmp[, 2:3], tmp[, x]),
Hi,
I have a column in a dataframe in the form of:
as.vector(SLDATX[1:20])
[1] 1/6/1986 1/17/1986 2/2/1986 2/4/1986 2/4/1986
[6] 2/21/1986 3/6/1986 3/25/1986 4/6/1986 4/10/1986
[11] 4/23/1986 4/30/1986 5/8/1986 5/29/1986 6/15/1986
[16] 6/18/1986 6/23/1986 6/29/1986 7/16/1986 7/25/1986
Ko-Kang Kevin Wang wrote:
Hi,
I have a column in a dataframe in the form of:
as.vector(SLDATX[1:20])
[1] 1/6/1986 1/17/1986 2/2/1986 2/4/1986 2/4/1986
[6] 2/21/1986 3/6/1986 3/25/1986 4/6/1986 4/10/1986
[11] 4/23/1986 4/30/1986 5/8/1986 5/29/1986 6/15/1986
[16] 6/18/1986 6/23/1986
Can I use nlm through the R API from C?
It's an internal, so I believe it's not supported yet.
If not, any suggestions on a workaround.
__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
I would like to announce that my DataLoad utility which can be found at:
http://www.vsn-intl.com/genstat/downloads/datald.htm
has been updated to support R data.frames. DataLoad reads a
large variety of data formats (listed on the web page) and
can convert these to ASCII or XDR data.frames. Only
I'm using Visual Studio 6.
I set my project to get headers from C:\Program Files\R\rw1081\src\include and get
libraries from C:\Program Files\R\rw1081\src\gnuwin32(though i'm not running
gnuwin at all, just Visual Studio 6)
And I added the rdll.lib from C:\Program
What is the correct syntax for running a batch program (test.r) from the a:
drive. I have tried no quotes, single and double quotes around a:test.r to
no avail.
Thanks in anticipation.
John Byrne.
__
[EMAIL PROTECTED] mailing list
Hi,
Has anyone got any R code (or are there any packages) that calculates
the great circle distance between two geographical (lat, lon) positions?
Cheers
Toby Patterson
Pelagic Ecosystems Research Group
CSIRO Marine Research
Email: [EMAIL PROTECTED]
Have you seen Online calculations Downloadable spreadsheets to perform
Geodetic Calculations.
at http://www.ga.gov.au/nmd/geodesy/datums/calcs.jsp ?
Duncan
*
Dr. Duncan Mackay
School of Biological Sciences
Flinders University
GPO Box 2100
Adelaide
S.A.
strptime takes a character input and produces a POSIXlt output so
the format you specify to strptime is the format of the input,
not the output:
format( strptime(10/22/1986, %m/%d/%Y), %Y-%m )
---
Date: Wed, 26 Nov 2003 13:23:45 +1300 (NZDT)
From: Ko-Kang Kevin Wang [EMAIL PROTECTED]
To:
58 matches
Mail list logo