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
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.
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]
_
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
https://www.stat.
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 Files\R\rw1081\src\gnuwin3
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 t
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
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"
[
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
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"]),
f
[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.
That'
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
> r
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
http://www.scieng.flinders.edu.au/
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? [0,1],[10001,2],[21,3000
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
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
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
regression, like so:
> f <- yresp ~
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 s
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
Du
> 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
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 ha
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
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 o
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 v
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 d
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 pa
"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,
>
> + random
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, execu
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 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
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
> "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
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)
Perhaps you could take the trouble to read the error message, which is
Error in inherits(x, "factor") : evaluation is nested too deeply:
infinite recursion?
The evaluation depth is controlled by options(expressions=). Increasing
it allows your code to run, albeit very slowly.
On Tue, 25
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 rep
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 -82229
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 fix
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
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
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
Hi all,
I am playing around with latin squares, and wrote a recursive function
that searches for valid combinations.
Apart from the fact that there are very many, I run into troubles
beginning with size 10x10 because the recursion depth becomes too large
(max of 10x9-1=89 in this case).
Why is
Hi,
I think this is caused by changes to subscripting for data.frames
made in 1.8.0 (grep "Subscripting" in the 1.8.0 NEWS file).
I'll fix this in an RMySQL update I'm currently working on
(and soon to be available).
Regards,
--
David
Christian Schulz wrote:
> Hi,
>
> i getting following err
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 T
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
O
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)
> 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 othe
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 th
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 pag
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 ar
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 /usr/lib/R/l
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]]
_
> 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 sav
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?
P.S.
It's possible calculate approximate the
memory demand for different settings with RF?
Many thanks & regar
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, fr
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:
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.
Thanks
Pascal
__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/
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
day
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
https://www.stat.math.ethz.ch/mailman/listi
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,
> 8682289,74
> 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
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,
5065887,4594147
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, furt
63 matches
Mail list logo