On Sun, 14 Dec 2003, Douglas Trainor wrote:
Somewhere along the line, you have been confused.
You're in good company though. Factor analysis
and PCA are different entities entirely.
Not in SPSS, where the same command is used for both (although we were not
told anything like enough about
On 14 Dec 2003, Peter Dalgaard wrote:
Thomas W Blackwell [EMAIL PROTECTED] writes:
Murray -
If you could guarantee that all of the email addresses have
exactly one occurrence of the @ character in them, then
something like
snip
Otherwise, try something like this (I
Did you make sure that you loaded the library that contains the C
code on the remote hosts?
best,
-tony
Faheem Mitha [EMAIL PROTECTED] writes:
On Mon, 1 Dec 2003, A.J. Rossini wrote:
Faheem Mitha [EMAIL PROTECTED] writes:
So, can this (parallelization at the C level) be done without
Prof Brian Ripley wrote:
On Sun, 14 Dec 2003, Douglas Trainor wrote:
Somewhere along the line, you have been confused.
You're in good company though. Factor analysis
and PCA are different entities entirely.
Not in SPSS, where the same command is used for both (although we were not
told
On Sun, 14 Dec 2003, A.J. Rossini wrote:
Did you make sure that you loaded the library that contains the C
code on the remote hosts?
No, I only loaded the library on the local node (master).
However, I'm not sure how I should do so. If I simply start up R on the
remote hosts and load up
Faheem Mitha [EMAIL PROTECTED] writes:
On Sun, 14 Dec 2003, A.J. Rossini wrote:
Did you make sure that you loaded the library that contains the C
code on the remote hosts?
No, I only loaded the library on the local node (master).
However, I'm not sure how I should do so. If I simply
Whoops. Exhaustion alert. The point is that R gets started up on the
remote hosts, but you need to initialize the library on them. So you
don't have to manually start up R on the remote hosts, but you do have
to tell the remote sessions what to compute (including loading any
libraries, etc).
Hi all,
I've got a problem with the nls function.
I have an adjustment which works when I fix one of the argument of my
function (Xo=150) :
*Xo*=150
f- function (tt*,Xo*,a,b)ifelse(tt*Xo*,a*exp(-b**Xo*),a*exp(-b*tt))
Dear Sir,
I am a user of R, and I just downloaded rw1081.exe
from your site.
When I install the said file, it reports that a file
is missing. I dont know the cause. Its the only
mising file because when I press ignore, everyting
installs just fine.
Anyways, Sir, can you send me the lacking
Dear chenu,
I am not going to see your code with attention (also I do not understand the
`*' symbol you used), however it looks a changepoint-type problem.
The package segmented (on CRAN) uses a piecewise linear parameterization to
fit regression models with breakpoints.
Hope this helps,
best,
Hello there fellow R-users,
I have received some data which comes in the following format:
example1-200301
The first 4 digits correspond to the year and the remaining 2 digits
correspond to the week of the year.
I have tried to convert this to a date by using strptime as follows:
On Mon, 15 Dec 2003, Wayne Jones wrote:
Hello there fellow R-users,
I have received some data which comes in the following format:
example1-200301
The first 4 digits correspond to the year and the remaining 2 digits
correspond to the week of the year.
I have tried to convert this
Have you tried the 'sm.density' function from the sm library?
I used it for a dataset which 'only' had 13 points.
I'm new to R and am trying to perform a simple, yet
problematic task. I
have two variables for which I would like to measure the
correlation and
plot versus each other.
I can't answer all your questions right now, but I can answer the first:
Have you considered gamma: n! = gamma(n+1)
gamma(1+1:6)
[1] 1 2 6 24 120 720
Regarding the other two, have you consulted Pinhiero and Bates
(2000) Mixed-Effects Models in S and S-Plus (Springer)? Bates
Prof Brian Ripley wrote:
That is not a complete date! Which day of the week is it?
example1-200301
strptime(paste(example1, 1, format=%Y%U %d)
will work, and you need to do something like that to resolve the
ambiguity.
Doesn't seem to work:
- 6th day of 5th week:
strptime(2003 05 06,
BR == Barry Rowlingson [EMAIL PROTECTED] writes:
BR Ah ha. Use '%w' for 'day of week':
strptime(2003 05 06, format=%Y %U %w)
BR [1] 2003-02-02
BR I dont have a calendar to hand to check that 2 Feb is the 6th
BR day of the 5th week
It's not, it is the first day of the
On Mon, 15 Dec 2003, Barry Rowlingson wrote:
Prof Brian Ripley wrote:
That is not a complete date! Which day of the week is it?
example1-200301
strptime(paste(example1, 1, format=%Y%U %d)
will work, and you need to do something like that to resolve the
ambiguity.
Dear useRs:
First of all I would like to thank all the responses.
I've an error with package spdep.
I am working with a Windows XP machine (AMD-2000-XP RAM-256DDR) and 1.8.0.
R-version and when I try to load spdep appear the following error:
library(spdep)
Error in loadNamespace(i, c(lib.loc,
Berwin Turlach [EMAIL PROTECTED] writes:
BR == Barry Rowlingson [EMAIL PROTECTED] writes:
BR Ah ha. Use '%w' for 'day of week':
strptime(2003 05 06, format=%Y %U %w)
BR [1] 2003-02-02
BR I dont have a calendar to hand to check that 2 Feb is the 6th
BR day of the
On Mon, 15 Dec 2003, Perez Martin, Agustin wrote:
Dear useRs:
First of all I would like to thank all the responses.
I've an error with package spdep.
I am working with a Windows XP machine (AMD-2000-XP RAM-256DDR) and 1.8.0.
R-version and when I try to load spdep appear the following
spdep depends on package 'maptools', install it.
- Original Message -
From: Perez Martin, Agustin [EMAIL PROTECTED]
To: lista R help (E-mail) [EMAIL PROTECTED]
Sent: Monday, December 15, 2003 1:17 PM
Subject: [R] Error with spdep
Dear useRs:
First of all I would like to thank all the
Hi! I need to do a Logistic Regresion.
I installed the Design Package. I tried lrm(... ) but
the program said Object na.delete is not found.
What can I do?
Thanks!!!
Marcela
__
[EMAIL PROTECTED] mailing list
Since Spencer Graves already answered the factorial questions, I'll try to
answer one of the other two:
On Monday 15 December 2003 05:17, [EMAIL PROTECTED] wrote:
To anyone who can help,
Intelligent question (1) I keep on trying to fit a linear mixed model in R
using 'lme(y~fxd.dsgn, data =
On Mon, 15 Dec 2003 04:54:12 -0800 (PST)
Marcela XX [EMAIL PROTECTED] wrote:
Hi! I need to do a Logistic Regresion.
I installed the Design Package. I tried lrm(... ) but
the program said Object na.delete is not found.
What can I do?
Thanks!!!
Marcela
The package has a lot of
Dear R user,
I asked for packages which should be mentioned in a ecological related
presentation of R and would like to summarize the replies.
First of all E. Paradis pointed out, that he would prefer to see R presented
in a different way, please read his whole reply
Hi,
I believe the following is a very minor typo in the documentation for
merge() :
all logical; all=L is shorthand for all.x=L and all.y=L
^^ ^
Looks like 'L' should be 'FALSE'.
Someone please alert me if I am wrong. If I am not wrong, is
On Sat, 13 Dec 2003 [EMAIL PROTECTED] wrote:
On 13-Dec-03 Martin Maechler wrote:
In general, use
== for testing equality of integer numbers (of type integer
or not)
I hope this is not a suggestion to avoid usage like
which(x == max(x))
when x is a vector of
JARogers == Rogers, James A [PGRD Groton] Rogers
on Mon, 15 Dec 2003 10:08:50 -0500 writes:
JARogers I believe the following is a very minor typo in
JARogers the documentation for merge() :
JARogers all logical; all=L is shorthand for all.x=L and all.y=L
JARogers
On 15 Dec 2003, Faheem Mitha stated:
So, the question would be how to get the R slaves to load up the library
and I don't see any obvious way of doing this (using snow functions).
Say you need boot on each slave, call this on the master,
clusterEvalQ (cl, library (boot))
where cl is the
On Mon, 15 Dec 2003, A.J. Rossini wrote:
Yes, and no.
Re-read the CPH-statcomp lab, and look at the bootstrap example, which
solves the same problem.
Look carefully -- it has to initialize the library on each node.
If you are just loading the library manually, just do it on each node;
On Mon, 15 Dec 2003, Na Li wrote:
On 15 Dec 2003, Faheem Mitha stated:
So, the question would be how to get the R slaves to load up the library
and I don't see any obvious way of doing this (using snow functions).
Say you need boot on each slave, call this on the master,
Eryk Wolski wrote:
How to plot if possible flury faces in R?
/Eryk
__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
Some days ago I was looking for a function to represent multivariate data by
faces, too. For
On 15-Dec-03 Thomas Lumley wrote:
One reason that which.max() exists is that we cannot guarantee
which(x==max(x)) to work. It is possible, though rather unlikely, for
there to be no x such that x==max(x). One reason is the unpredictable
use of 10-byte wide floating point registers on Intel
Hi,
I am getting some weird results here and I think I am missing something.
I am trying to program a function that for a set of random variables
drawn from uniform distributions plots that distribution of the second
order statistic of the ordered variables. (ie I have n uniform
distributions
If you don't want to rely on the week calculations in R (which in
turn may depend on the OS?) you could try doing it yourself. The
formula is pretty simple.
Suppose:
fw = day of week of first day of each week (0=Sun, 1=Mon, etc.)
fy = day of week of first day of year (to be calculated given
The order statistics have a beta distribution, so pbeta is all you need.
On Mon, 15 Dec 2003, Eugene Salinas (R) wrote:
I am getting some weird results here and I think I am missing something.
I am trying to program a function that for a set of random variables
drawn from uniform
Dear R-user,
I started a few month ago with R and now I think that it is a great program
but in the beginning I was close to give up, because I did not manage to
read
in my files, did not know how the indexing works, where my results are or
did
not understand how R saves my data and so
On Mon, 15 Dec 2003 [EMAIL PROTECTED] wrote:
On 15-Dec-03 Thomas Lumley wrote:
One reason that which.max() exists is that we cannot guarantee
which(x==max(x)) to work. It is possible, though rather unlikely, for
there to be no x such that x==max(x). One reason is the unpredictable
use
R UseRs:
I can't figure out how to append intermediate terms inside a For loop. My
use of append(), various indexing, and use of data frames, vectors,
matrices has been fruitless. Here's a simplified example of what I'm
talking about:
i - 1
for(i in 10) {
v[i] - i/10
}
v
[1] 1 NA NA
[EMAIL PROTECTED] writes:
I can't figure out how to append intermediate terms inside a For loop. My
use of append(), various indexing, and use of data frames, vectors,
matrices has been fruitless. Here's a simplified example of what I'm
talking about:
i - 1
for(i in 10) {
v[i] -
Prof Brian Ripley [EMAIL PROTECTED] writes:
On 15 Dec 2003, Peter Dalgaard wrote:
Neither of those definitions coincide with ISO, BTW.
ISO 8601 week of the year (as distinct from ISO C where these come from)
is %V, of course.
Doesn't work with RedHat 8 though...
--
O__
Thanks (to those who responded) for the help.
Yeah, I forgot to make it a sequence. Must have overlooked that.
Thanks again.
-Jason
[[alternative HTML version deleted]]
__
[EMAIL PROTECTED] mailing list
Hi,
I'm a bit confused how julian() works. If I understand right, it returns
the number of days since the origin.
I have a vector:
SLDATX[1:10]
[1] 1986-01-06 1986-01-17 1986-02-02 1986-02-04
[5] 1986-02-04 1986-02-21 1986-03-06 1986-03-25
[9] 1986-04-06 1986-04-10
And when I did:
[diverted from R-help to R-devel; please follow up on R-devel!]
ChrisH == Christian Hennig [EMAIL PROTECTED]
on Thu, 11 Dec 2003 15:13:52 +0100 (MET) writes:
ChrisH Hi, this is rather a (presumed) bug report than a
ChrisH question because I can solve my personal statistical
I hope I am not telling you things you already know. If so, I apologize
in advance.
There are several C-library addons available to try to deal with the
problem that comparisons of floating point numbers can be
unpredictable. I think your example with the greater than sign would
not be a
Hi!
Old story again. I put it away for a while because there are always other thinks to do.
But I cant deny that I still like to like to comment assignment functions.
Under linux I get.
Cannot handle Rd file names containing ''.
These are not legal file names on all R platforms.
Please rename
Ko-Kang Kevin Wang wrote:
Hi,
I'm a bit confused how julian() works. If I understand right, it returns
the number of days since the origin.
I have a vector:
SLDATX[1:10]
[1] 1986-01-06 1986-01-17 1986-02-02 1986-02-04
[5] 1986-02-04 1986-02-21 1986-03-06 1986-03-25
[9]
Dr Murray Jorgensen ([EMAIL PROTECTED]) wrote:
I have some email addresses that I would like to sort in reverse
lexicographic order so that addresses from the same domain will be
grouped together. How might that be done?
Because he wants addresses from the same domain
What you can do to handle this timezone problem is either to use
POSIXt with GMT or use chron (which does not use timezones so
can't cause problems like this):
Suppose:
SLDATX - c( 1986-01-06, 1986-01-17, 1986-02-02, 1986-02-04,
,1986-02-04, 1986-02-21, 1986-03-06, 1986-03-25,
,1986-04-06,
Thanks! chron() is very useful indeed.
Just out of interest, is it possible to do, say in this case, the number
of months (or quarters) after January 1986? i.e. use a different time
interval?
On Mon, 15 Dec 2003, Gabor Grothendieck wrote:
What you can do to handle this timezone problem is
On Tue, 16 Dec 2003, Richard A. O'Keefe wrote:
Here's a strrev() I knocked together quickly:
strrev -
function (s) paste(rev(strsplit(s, character(0))[[1]]), collapse=)
If anyone can tell me how to vectorise this, I would be glad of the lesson.
strrev- function(ss) {
Assuming you want to find a whole number months and quarters:
require(chron)
z - chron( SLDATX, format=y-m-d )
months - with( month.day.year(z), 12*(year-1986)+month-1 )
quarters - months %/% 3
---
Date: Tue, 16 Dec 2003 13:58:13 +1300 (NZDT)
From: Ko-Kang Kevin Wang [EMAIL
Hi all,
I am trying to install add-on in R (rw1070) to work with WinEdit. Libraries
Swinregistry and Rwinedt have been installed via Install Package(s) from local zip
files from R. However, when I run library(Rwinedt) I get the following messages:
library(Rwinedt)
Loading required package:
I wrote:
If anyone can tell me how to vectorise this, I would be glad of the lesson.
where this was
strrev -
function (s) paste(rev(strsplit(s, character(0))[[1]]), collapse=)
Thomas Lumley [EMAIL PROTECTED] suggested
strrev- function(ss) {
sapply(lapply(
Hi,
I am new to R (I have most of my experience in SAS and SPSS). I was
wondering if anyone has used both Resampling Stats and R, and could comment
on strengths/relationships. Also, I have no clue on how to do the various
examples from the book Resampling: The New Statistics in R. Can anyone
Here is a way to do it without using apply. sep must be set to
a character not in any of the strings. Below we show its much
faster than using apply yet gives the same answer.
strRev - function(x, sep = \10) {
z - unlist( strsplit( paste( x, sep, sep= ), ) )
z - unlist(
Brandon Vaughn wrote:
...
I am new to R (I have most of my experience in SAS and SPSS). I was
wondering if anyone has used both Resampling Stats and R, and could comment
on strengths/relationships.
There are a few add-on packages for resampling with R. boot is the
one I've used, and can
57 matches
Mail list logo