James MacDonald wrote:
Thanks Uwe,
What is the theoretical limit for 32 bit Windows anyway?
BTW, --max-mem-size=2000M works (although maybe not all usable?)
I don't know the usable amount exactly, one has to look into Microsoft's
docs. It must be a bit smaller than 2GB (differently from the
Richard A. O'Keefe wrote:
I am puzzled by the advice to use is.na(x) - TRUE instead of x - NA.
?NA says
Function `is.na-' may provide a safer way to set missingness. It
behaves differently for factors, for example.
However, MAY provide is a bit scary, and it doesn't say WHAT the
Dear All
Is there a SIMCA (Soft Independent Modelling Class Analogy) implementation
on R or does anyone know if is it possible to replicate the SIMCA algorithm
using existing R functions?
Thanks
Mike White
__
[EMAIL PROTECTED] mailing list
Hi
In trying to fit a linear model , I use the leaps() function to determine wich
predictors I should include in my model.
I would like to plot the Mallow's Cp criteria against p with the indexes of selected
model variates as points labels
Is there already such a function? (I could not find it)
Hello,
I have been trying to write a small program to generate automatic plots. What
I want is to draw boxplots for some variables in my data frame, contrasting
with a variable called 'missing' that has value 1 if some variable in a
concrete case has at least one missing value, in order to
Note this behaviour:
a-a
a-NA
mode(a)
[1] logical
a-a
is.na(a) - T
mode(a)
[1] character
However after either way of assigning NA to a, is.na(a) is true,
and it prints as NA, so I can't see it's ever likely to matter. [Why
do I say these things? Expect usual flood of examples where it
Anne Piotet wrote:
Hi In trying to fit a linear model , I use the leaps() function to
determine wich predictors I should include in my model. I would like
to plot the Mallow's Cp criteria against p with the indexes of
selected model variates as points labels
Is there already such a function?
Xavier Fernández i Marín wrote:
...
varlist - c(var1, var2, var3, var4, ...)
Instead of a character vector with the names, it'd make life easier if
you had a list of the vectors...
# make sure you use the naming - makes life easier later.
mylist - list(var1=var1, var2=var2, var3=var3, var4=var4)
On Wed, 8 Oct 2003, Simon Fear wrote:
Note this behaviour:
a-a
a-NA
mode(a)
[1] logical
a-a
is.na(a) - T
mode(a)
[1] character
However after either way of assigning NA to a, is.na(a) is true,
and it prints as NA, so I can't see it's ever likely to matter. [Why
do I say these
Dear list,
I can not understand why the expression in
the subject does not work correct:
dcrn[which(fn == inve[2])]
numeric(0)
inve[2]
[1] 406.7
dcrn[which(fn == 406.7)]
[1] 1.3994e-07 1.3988e-07 1.3953e-07 1.3966e-07 1.3953e-07 1.3968e-07
Is this a kick self problem or an bug?
Thaks very
On Wednesday 08 October 2003 11:27, Thomas Bock wrote:
Dear list,
I can not understand why the expression in
the subject does not work correct:
dcrn[which(fn == inve[2])]
numeric(0)
inve[2]
[1] 406.7
dcrn[which(fn == 406.7)]
[1] 1.3994e-07 1.3988e-07 1.3953e-07 1.3966e-07
Well, that's a convincing argument, but maybe
it's the name that's worrying some of us. Maybe it would be
more intuitive if called set.na (sorry, I mean setNA).
Also is.na- cannot be used to create a new variable of
NAs, so is not a universal method, which is a shame for its
advocates.
I
Thomas Bock [EMAIL PROTECTED] writes:
Dear list,
I can not understand why the expression in
the subject does not work correct:
dcrn[which(fn == inve[2])]
numeric(0)
inve[2]
[1] 406.7
dcrn[which(fn == 406.7)]
[1] 1.3994e-07 1.3988e-07 1.3953e-07 1.3966e-07 1.3953e-07 1.3968e-07
Thanks, your advice worked. I don't have much experience with maths, and
therefore tried to stay away from dealing with optimization, but going
down to this level opens a lot of possibilities. For the record, the
code I used, as you suggested:
###
shape - mean(data)^2/var(data)
scale
Hi,
For a n-level factor, I'd like to specify the first contrast and have
the remaining n-2 constructed automatically so that the set is
orthogonal. I then test the contrasts with summary.lm(anova-object).
In S-Plus, the following works:
y.anova - aov( y ~ C(CO2,c(1,0,-1)) )
On Wed, 8 Oct 2003, Pascal A. Niklaus wrote:
Hi,
For a n-level factor, I'd like to specify the first contrast and have
the remaining n-2 constructed automatically so that the set is
orthogonal. I then test the contrasts with summary.lm(anova-object).
In S-Plus, the following works:
I'm sorry, but I don't have time to read all your code. However,
I saw that you tested for x alpha in your Pareto distribution
example. Have you considered reparameterizing to estimate log.del =
log(alpha-min(x))? Pass log.del as part of the vector of parameters to
estimate, then
Your question has been answered by Achim and Peter Dalgaard (at least).
Just a note:
Using
a[which(logic)]
looks like a clumsy and inefficient way of writing
a[ logic ]
and I think you shouldn't propagate its use ...
Martin Maechler [EMAIL PROTECTED]
Hi All,
I'm running R 1.7.1, and I've installed some additional packages such a
Bioconductor. Do I've to re-install all the additional packages when ugrading
to R 1.8.0 (i.e. are there compile in dependencies)?
thanks for your help,
Arne
I've rolled up R-1.8.0.tgz a short while ago. This is a new version
with major changes (see below). Notably, the Macintosh version for OS
X has been substantially improved; the old Carbon interface is no
longer being supported.
Also notice that the underscore will no longer work as an assignment
Thanks for the reply. Below is the solution and the S-Plus and R code
that does the same (for documentation).
I can't reproduce that in S-PLUS 6.1, and it is not as documented:
In S-Plus 2000, C() complements the contrast matrix with orthogonal
contrasts if only the first is given.
CO2 -
Most packages do not need to be re-installed, but R 1.8.0 has some nicer
formatting of help pages so you may want to do so, and I am reinstalling
all packages that make use of the methods package and so have saved
images.
For Bioconductor you need to ask on their list, but my understanding is
Thanks for the help, the wrapper function was very useful. I managed to
solve the problem using Spencer Graves' suggestion. I am analyzing the
interarrival times between HTTP packets on a campus network. The dataset
actually has more than 14 Million entries! It represents the traffic
generated by
I have a question regarding bootstrap coverage. I am trying to understand the
benefits of using the bootstrap for small sample sets. To do this I created a normal
population and then picked 10 from the populations and applied both traditional
statistical methods and the Bootstrap (bcanon,
Hello,
thanks for the tips on updating packages for 1.8.0. The updating is a real
problem for me, since I've to do it sort of manually using my web-browser or
wget. I'm behind a firewall that requires http/ftp authentification (username
and passwd) for every request it sends to a server outside
On Wed, 8 Oct 2003, Pascal A. Niklaus wrote:
Thanks for the reply. Below is the solution and the S-Plus and R code
that does the same (for documentation).
I can't reproduce that in S-PLUS 6.1, and it is not as documented:
In S-Plus 2000, C() complements the contrast matrix with
Hi,
my out-of-the-box installation of R-1.8.0 on Tru64 (OSFv5.1) results in
an tremendously increased verbosity of Rdinfo.
Is it really intended that Rdinfo complaints about minor mistakes like
missing empty lines at the end given the very variable quality of
Rd-pages?
Is there a way to
Also, presumably is.na- could be redefined by the user for particular
classes so if you got in the habit of setting NAs that way it would
generalize better.
---
Date: Wed, 8 Oct 2003 11:49:29 +0100 (BST)
From: Prof Brian Ripley [EMAIL PROTECTED]
I don't think it can ever `go wrong', but it
At Wednesday 03:06 PM 10/8/2003 +0200, Martin Maechler wrote:
Your question has been answered by Achim and Peter Dalgaard (at least).
Just a note:
Using
a[which(logic)]
looks like a clumsy and inefficient way of writing
a[ logic ]
and I think you shouldn't propagate its use ...
Are you interested in turning that into a monitor, processing each
day's data sequentially or even each entry as it arrived? If yes, you
may wish to evaluate the Foundations of Monitoring documents
downloadable from www.prodsyse.com. If you have any questions about
that, I might be able
Sorry, I didn' mean it the nasty way. I wouldn't have been surprised if the
R-team had told me the authentification with the firewall is my problem (i.e.
a special case that cannot be dealt with by th R-team).
Yess, and off course I should have had a much closer lookk into the docu.
Thanks again
On Wed, 8 Oct 2003 [EMAIL PROTECTED] wrote:
Hello,
thanks for the tips on updating packages for 1.8.0. The updating is a real
problem for me, since I've to do it sort of manually using my web-browser or
wget. I'm behind a firewall that requires http/ftp authentification (username
and
Dear R and sweave users
A further problem, which I couldn't resolve, using the manual: In R I
use the split.screen command to put e.g. two timecourses one above the
other into one plot:
split.screen(c(2,1))
screen(1)
plot(stick,type='h', col=red,lwd=2)
screen(2)
Cloud() function does not display anything with R1.8.0beta
in WindowsXP ...
Does any one noticed this ?
others functions from lattice seem working properly.
does it work in the final 1.8.0 for windows ?
__
[EMAIL PROTECTED] mailing list
On Wed, Oct 08, 2003 at 04:25:13PM +0200, [EMAIL PROTECTED] wrote:
Hello,
thanks for the tips on updating packages for 1.8.0. The updating is a real
problem for me, since I've to do it sort of manually using my web-browser or
wget. I'm behind a firewall that requires http/ftp
Here are some different ways of doing this. Don't know whether any
could be considered superior to the others.
# y[x==5] regarding NAs in x as not matching
x - c(5, NA, 7, 5, NA, 3)
y - c(1, 2, 3, 4, 5, 6)
subset(y,x==5)
y[x %in% 5]
y[x %in% c(5)]
y[which(x==5)]
---
Date: Wed, 08 Oct 2003
Dear all;
Installing the GLMMGibbs package to my Solaris Unix box, I got an compiling
error:
ars.c:497:10: missing terminating character
ars.c: In function `dump_arse':
ars.c:498: error: parse error before mylesj
.
The
Christoph Lehmann wrote:
Dear R and sweave users
A further problem, which I couldn't resolve, using the manual: In R I
use the split.screen command to put e.g. two timecourses one above the
other into one plot:
split.screen(c(2,1))
screen(1)
plot(stick,type='h', col=red,lwd=2)
screen(2)
Achim Zeileis wrote:
On Wednesday 08 October 2003 11:27, Thomas Bock wrote:
...
I can not understand why the expression in
the subject does not work correct:
dcrn[which(fn == inve[2])]
numeric(0)
inve[2]
[1] 406.7
...
1.) `==' comparisons have a certain tolerance
2.) the print output is
This seems to me to be a special case of the general problem of a
parameter on a boundary. Another example is the case of a variance
component that is zero. For this latter problem, Pinhiero and Bates
(2000) Mixed-Effects Models in S and S-Plus (Springer, sec. 2.4.1)
present simulation
Whoops. Hit send too quickly.
Jason Turner wrote:
tol - sqrt(.Machine$double.eps)
dcrn[(fn - inve[2]) tol]
that should be
dcrn[abs(fn - inve[2]) tol]
--
Indigo Industrial Controls Ltd.
http://www.indigoindustrial.co.nz
64-21-343-545
[EMAIL PROTECTED]
Dear listers,
I have two questions:
(1)
Is there a way in R to change the base-n of the calculations. I wnat to run
some calculations either in binary (base-2) or base-4. Is there a way to
specify that in R - to chnage from the decimal?
(2)
I also want to extract the digits from a larger number
Spencer Graves [EMAIL PROTECTED] writes:
This seems to me to be a special case of the general problem of
a parameter on a boundary.
Umm, no...
I have this problem with my data. In a GLM, I have 269 zeroes and
only 1 one:
I don't think that necessarily gets you a parameter estimate
Thanks, Peter: You are absolutely correct. Thanks again for the
correction. Spencer Graves
Peter Dalgaard BSA wrote:
Spencer Graves [EMAIL PROTECTED] writes:
This seems to me to be a special case of the general problem of
a parameter on a boundary.
Umm, no...
I have this
Have got anybody experience using winedt and R in Linux
- perhaps with wine
Or exist another editor with the ability
to parse r code into R in Linux,
because k-edit etc. using only syntax-highlithing ???
many thanks, christian
[[alternative HTML version deleted]]
On Wednesday 08 October 2003 22:02, Christian Schulz wrote:
Have got anybody experience using winedt and R in Linux
- perhaps with wine
Or exist another editor with the ability
to parse r code into R in Linux,
because k-edit etc. using only syntax-highlithing ???
I guess you have looked at
From: Andrej Kveder [mailto:[EMAIL PROTECTED]
Dear listers,
I have two questions:
(1)
Is there a way in R to change the base-n of the calculations.
I wnat to run some calculations either in binary (base-2) or
base-4. Is there a way to specify that in R - to chnage from
the decimal?
Simon Fear [EMAIL PROTECTED] suggested that
a-a
a-NA
mode(a)
[1] logical
a-a
is.na(a) - T
mode(a)
[1] character
might be a relevant difference between assigning NA and using is.na.
But the analogy is flawed: is.na(x) -
Concerning x[i] - NA vs is.na(x[i]) - TRUE
Brian Ripley wrote:
I don't think it can ever `go wrong', but it can do things other
than the user intends.
If the user writes x[i] - NA, the user has clearly indicated his intention
that the i element(s) of x should become NA. There
Good afternoon,
We currently have R installed on our HP Superdome. However, we are getting ready to
migrate from RISC to Itanium 2 chips running HP-UX (not Linux).
Does the latest version of R run on HP-UX Itanium 2?
Any information would be greatly appreciated.
Michael
Achim Zeileis [EMAIL PROTECTED] wrote:
R x - 406.7 + 1e-20
R x
[1] 406.7
R x == 406.7
[1] TRUE
that is
1.) `==' comparisons have a certain tolerance
No, all.equal() supports tolerance, == does not.
Consider
Hi everybody
I want to specify the contrasts to build a cell means model on LME (each
coefficient is the mean value of the independent variable for the specific
category of a factor variable) when there are several factors as fixed
effect in the model and also interactions between them. Can
Tongue in cheek
But surely
is.na(x) - is.na(x)
is clearer than
x[is.na(x)] - NA
(neither of which is a no-op).
/Tongue in cheek
__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
I'd just like to thank the R-Core team for the new version (all OS's), and especially
Stefano Iacus for his work on the old Carbon MacOS port and his continuing work on the
OS X port.
Cheers,
Simon.
Simon Blomberg, PhD
Depression Anxiety Consumer Research Unit
Centre for Mental Health
Hi, R-users:
Last week I send a request for help to this list. I have receive until now
two kindly responses from Spencer Graves and Renauld Lancelot. They both
point interesting things to fit an adequate model to my data but
unfortunately
it persists without a satisfactory solution.
I
On 08-Oct-03 Liaw, Andy wrote:
From: Andrej Kveder [mailto:[EMAIL PROTECTED]
I have two questions:
(1)
Is there a way in R to change the base-n of the calculations.
I wnat to run some calculations either in binary (base-2) or
base-4. Is there a way to specify that in R - to chnage from
Dear List members:
I'm using R1.7.1 (Windows 2000) and having difficulty with scoping.
I've studied the FAQ and on-line manuals and think I have identified
the
source of my difficulty, but cannot work out the solution.
For the purposes of illustration. I have three functions as defined
below:
Check out:
http://www.wiwi.uni-bielefeld.de/~wolf/software/R-wtools/decodeencode/decodeencode.rev
---
Date: Wed, 8 Oct 2003 21:39:50 +0200
From: Andrej Kveder [EMAIL PROTECTED]
Dear listers,
I have two questions:
(1)
Is there a way in R to change the base-n of the calculations. I wnat
It seems like you want in fnB
get(AA$first, envir = parent.frame(1))
but I'm entirely clear on why your original function doesn't work. My
understanding was that get() should search through the parent frames.
-roger
Peter Alspach wrote:
Dear List members:
I'm using R1.7.1 (Windows 2000) and
Hi,
I am new to R, but have extensive experience in
matlab.
I have searched on the web and in Venabels Ripley
book but I was unable to find the equivalent of the
eps function in matlab.
eps returns the Floating point relative accuracy.
thanks
__
Check out ?.Machine
-roger
christoff pale wrote:
Hi,
I am new to R, but have extensive experience in
matlab.
I have searched on the web and in Venabels Ripley
book but I was unable to find the equivalent of the
eps function in matlab.
eps returns the Floating point relative accuracy.
thanks
Debian packages of R 1.8.0 were uploaded earlier for i386.
Binaries for alpha, ia64, hppa and powerpc are already in the package pool;
the arm, mipsel and s390 architecture have their built completled and will
be added to the pool shortly while the remaining architecures should follow
over the
On Wed, 8 Oct 2003, Roger D. Peng wrote:
It seems like you want in fnB
get(AA$first, envir = parent.frame(1))
but I'm entirely clear on why your original function doesn't work. My
understanding was that get() should search through the parent frames.
No, get() searches through the
On Wed, 8 Oct 2003, Roger D. Peng wrote:
It seems like you want in fnB
get(AA$first, envir = parent.frame(1))
but I'm entirely clear on why your original function doesn't work. My
understanding was that get() should search through the parent frames.
Where did you get that idea from?
I think I misinterpreted this section of the help file for get():
If `inherits' is `FALSE', only the first frame of the specified
environment is inspected. If `inherits' is `TRUE', the search is
continued up through the parent frames until a bound value of the
right mode is
Dear List,
I downloaded R for the first time yesterday, in the hopes that I might
deal more effectively with a complex repeated measures experimental
design involving inbred strains of laboratory mice. The design below,
somewhat simplified, cannot be computed with standard ANOVA, because
Dear all,
I am looking for a function that can scatter plot the residuals obtained
from a longitudinal model i.e. plot e_{i,j} e_{i,k} for
all j k = 1,..n ( I have 7 observations for each subject). something
similar to the pairs() function.
I can do it the long way by constructing a residual
Dear
I used to use R under windows. I always save my code in a txt file. When I
want to run something, I just copy and paste to the windows.
I recently need to do a big simulation work which will cost 5 or 6 days. I
am afraid the memory of windows can not cope with this situation.
I now want
Zhen Pang wrote:
...
I now want to run the code under unix. However, I do not know how to run
this code in txt file under unix. I just log in to the unix and enter
the R, what should I do next?
One more question is: if I log off my computer when R is running under
unix (i.e., disconnect my
69 matches
Mail list logo