Hello all...
I am really new to statistics and I am trying to figure out a way to
apply Chauvenet's criterion using the t-distribution on a set of numbers
in perl. I was unable to find a TDIST and TINV function for perl. I am
getting these functions from Excel. So, I figured that I
I get the second set each time, on Windows, using the build from CRAN.
Which BLAS are you using?
On Tue, 14 Sep 2004, Francisco Chamu wrote:
I start a clean session of R 1.9.1 on Windows and I run the following code:
library(MASS)
data(painters)
pca.painters - princomp(painters[ ,1:4])
On Mon, 13 Sep 2004, Tamas K Papp wrote:
The fact that as.integer(TRUE) is 1 (and that for FALSE, it gives
zero) is a really nice feature, eg when constructing piecewise
functions (for example, as in -x*(x0)+x*(x=0)) and for many other
things.
Since I haven't found a reference about this,
Tamas K Papp [EMAIL PROTECTED] writes:
On Mon, Sep 13, 2004 at 08:40:15AM -0400, Frank Samuelson wrote:
If anyone has a few extra CPU cycles to spare,
I'd appreciate it if you could verify a problem that I
have encountered. Run the code
below and tell me if it crashes your R before
We are hardly likely to know what those are in Excel. Possibly pt and qt,
but see help.search(Student t distribution) for where to look for what R
provides.
I also do not know what Chauvenet's criterion has to do with Student's t,
and
Tamas K Papp [EMAIL PROTECTED] writes:
The fact that as.integer(TRUE) is 1 (and that for FALSE, it gives
zero) is a really nice feature, eg when constructing piecewise
functions (for example, as in -x*(x0)+x*(x=0)) and for many other
things.
Since I haven't found a reference about this, I
Dear all
I use the library(netCDF) to read in NCEP data. The file I want to read has size 113
Mb.
When i try to read it I get the following message:
Error: cannot allocate vector of size 221080 Kb
In addition: Warning message:
Reached total allocation of 255Mb: see help(memory.size)
I get a
Prodromos Zanis [EMAIL PROTECTED] writes:
Dear all
I use the library(netCDF) to read in NCEP data. The file I want to
read has size 113 Mb. When i try to read it I get the following
message:
Error: cannot allocate vector of size 221080 Kb
In addition: Warning message:
Reached total
Prodromos Zanis wrote:
Dear all
I use the library(netCDF) to read in NCEP data. The file I want to read has size 113
Mb.
When i try to read it I get the following message:
Error: cannot allocate vector of size 221080 Kb
In addition: Warning message:
Reached total allocation of 255Mb: see
[EMAIL PROTECTED] wrote:
I'm using v 1.9.1 under Windoz XP.
Can I do the equivalent of ls -l on my R objects? R's ls() lists
only the names.
For example ll() in package gregmisc.
Uwe Ligges
Thanks!
David L. Reiner
Rho Trading
440 S. LaSalle St -- Suite 620
Chicago IL 60605
312-362-4963
On 14 Sep 2004, Peter Dalgaard wrote:
Prodromos Zanis [EMAIL PROTECTED] writes:
Dear all
I use the library(netCDF) to read in NCEP data. The file I want to
read has size 113 Mb. When i try to read it I get the following
message:
Error: cannot allocate vector of size 221080 Kb
On Tue, 14 Sep 2004, Uwe Ligges wrote:
[EMAIL PROTECTED] wrote:
I'm using v 1.9.1 under Windoz XP.
Can I do the equivalent of ls -l on my R objects? R's ls() lists
only the names.
For example ll() in package gregmisc.
But that does not give datetimes, since they are not recorded.
Prof Brian Ripley wrote:
On Tue, 14 Sep 2004, Uwe Ligges wrote:
[EMAIL PROTECTED] wrote:
I'm using v 1.9.1 under Windoz XP.
Can I do the equivalent of ls -l on my R objects? R's ls() lists
only the names.
For example ll() in package gregmisc.
But that does not give datetimes, since they are not
Can somebody remember me which is the command to erase
columns from a data frame?
Thanks Michele
___
http://it.seriea.fantasysports.yahoo.com/
__
[EMAIL PROTECTED] mailing list
Hi,
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of michele lux
Sent: Dienstag, 14. September 2004 10:44
To: [EMAIL PROTECTED]
Subject: [R] erase columns
Can somebody remember me which is the command to erase
columns from a data frame?
Thanks Michele
I
On Mon, 13 Sep 2004, Shuangge Ma wrote:
constrOptim() will do this, but it isn't a particularly efficient
algorithm when the number of constraints is large.
-thomas
Hello:
I have been trying to program the following maximization problem and would
definitely welcome some help.
the target
?subset
/E
*** REPLY SEPARATOR ***
On 9/14/2004 at 10:44 AM michele lux wrote:
Can somebody remember me which is the command to erase
columns from a data frame?
Thanks Michele
___
http://it.seriea.fantasysports.yahoo.com/
Perhaps the package Rcmdr is a compromise for the
people don't like command-line software.
christian
Am Dienstag, 14. September 2004 12:02 schrieb Paolo Ariano:
Hi *
i've done my anova anlysis but now i need post-hoc test, are these
included in R ?
I've a Big problem, working with people
On Tue, 14 Sep 2004, Paolo Ariano wrote:
i've done my anova anlysis but now i need post-hoc test, are these
included in R ?
Yes, in function TukeyHSD and in package multcomp for example. There are
worked examples in the MASS scripts.
I've a Big problem, working with people that don't like
Hello,
Frank Samuelson schrieb:
If anyone has a few extra CPU cycles to spare,
I'd appreciate it if you could verify a problem that I
have encountered. Run the code
below and tell me if it crashes your R before
completion.
library(lme4)
data(bdf)
dump-sapply( 1:5, function(i) {
fm -
Hi,
I am fairly new to GNU R !
At the moment I am doing an intensive learning on the basics of GNU
R-1.91, especially graphics like plots and alike, by reading the
introductory docs!
Well, except some occasional glitches (X11 output errors) everything
seems to be fine, thanks to developers for
Hi Thomas!
Try ?par at the R prompt, there you will get all the necessary
information to change the appearance of graphics.
Daniel
--
Daniel Hoppe
Department of Marketing
University of Vienna
Bruenner Strasse 72
1210 Vienna
Austria
-Original Message-
From: [EMAIL PROTECTED]
Thomas Schönhoff wrote:
Hi,
I am fairly new to GNU R !
At the moment I am doing an intensive learning on the basics of GNU
R-1.91, especially graphics like plots and alike, by reading the
introductory docs!
Well, except some occasional glitches (X11 output errors) everything
seems to be fine,
The X11 device has an argument `pointsize' (which may well mean pixel size
in a particular implementation of X11, as my laptop for example has width,
height and pointsize all much smaller than specified): just increase it.
The advice to look at ?par is incorrect if you want to scale everything,
Brian Ripley wrote:
We are hardly likely to know what those are in Excel. Possibly pt
and qt, but see help.search(Student t distribution) for where to
look for what R provides.
I also do not know what Chauvenet's criterion has to do with
Student's t, and
As others have said, this needs tools not CPU cycles: gctorture or valgrind.
Valgrind found (after a few seconds and on the first pass)
==23057== Invalid read of size 4
==23057==at 0x3CF4E645: ssc_symbolic_permute (Mutils.c:373)
==23057==by 0x3CF5BF75: ssclme_create (ssclme.c:168)
Il mar, 2004-09-14 alle 12:26, Prof Brian Ripley ha scritto:
BTW, I guess you are working on Windows, but have not told us what OS,
even. Please read the posting guide.
sorry, i use DebianGNU/Linux and my collegues windows
thanks
paolo
--
Paolo Ariano
Neuroscience PhD Student @ UniTo
Una
Prof Brian Ripley [EMAIL PROTECTED] writes:
The X11 device has an argument `pointsize' (which may well mean pixel size
in a particular implementation of X11, as my laptop for example has width,
height and pointsize all much smaller than specified): just increase it.
The advice to look at
Dear R-users,
I have a data matrix with 20 rows and 10 columns which is stored in the hard drive as
.csv file called c:\DataFile.csv and a 10 elements vector called xVec.
I would like to be able to copy and paste the information contained in xVec into (say)
the 2nd row of DataFile.csv
Prof Brian D Ripley wrote:
As others have said, this needs tools not CPU cycles: gctorture or valgrind.
Valgrind found (after a few seconds and on the first pass)
==23057== Invalid read of size 4
==23057==at 0x3CF4E645: ssc_symbolic_permute (Mutils.c:373)
==23057==by 0x3CF5BF75:
Many thanks to Prof. Ripley. The problem is that memory.profile does not
exist in *nix environment and there is probably a very good reason why.
I was reading help(Memory) and in the Details section :
You can find out the current memory consumption (the heap and cons
cells used as
On Tue, 14 Sep 2004, Adaikalavan Ramasamy wrote:
Many thanks to Prof. Ripley. The problem is that memory.profile does not
exist in *nix environment and there is probably a very good reason why.
memory.size?
I was reading help(Memory) and in the Details section :
You can find out the
If you only have simple objects in your function, you might be able to
use a function like
totalMem - function() {
sum(sapply(ls(all = TRUE), function(x) object.size(get(x / 2^20
}
which should give you a rough idea of the memory usage (in MB) in the
current environment.
-roger
Hi,
The following is a cut/paste from http://developer.r-project.org/200update.txt:
...
3) When a package is installed, all the data sets are loaded to see
what they produce. If this is undesirable (because they are
enormous, or depend on other packages that need to be installed
Is there a way of getting the argument list of a function from within
that function? For example, something like
f - function(x, y = 3) {
fargs - getFunctionArgList()
print(fargs) ## Should be `alist(x, y = 3)'
}
Thanks,
-roger
__
Ergh, yes, that's exactly it. I didn't realize you could use it in
that way.
Thanks,
-roger
Jeff Gentry wrote:
Is there a way of getting the argument list of a function from within
that function? For example, something like
Would formals() work for you here?
Paul, Thank you very much! They all works!
Sean
On Tue, 14 Sep 2004, Paul Murrell wrote:
Hi
xiang li wrote:
Also, I am wondering if there is any source where the expressions of
many symbols are collected.
Thanks you very much!!!
(Assuming you mean draw the angstrom symbol on a
Hello,
Thomas Schönhoff wrote:
Well, thanks, I'll have a look at your advices.
regards
Thomas
__
[EMAIL PROTECTED] mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Hello,
I have been trying to use the date library (mdy.date) to create notched boxplots, but
have not been successful. Can anybody help with the code for this command?
Thank you,
Stephanie
__
[EMAIL PROTECTED] mailing list
Hello-
I'm trying to do some repeated measures ANOVAs. In the past, using SAS,
I have used the framework outlined in Littell et al.'s SAS System for
Mixed Models, using the REPEATED statement in PROC MIXED to model
variation across time within an experimental unit. SAS allows you to
specify
Hi all,
I have a data.frame with the following colnames pattern:
x1 y11 x2 y21 y22 y23 x3 y31 y32 ...
I.e. I have an x followed by a few y's. What I would like to do is turn
this wide format into a tall format with two columns: x, y. The
structure is that xi needs to be associated with yij
I have run this on both Windows 2000 and XP. All I did was install
the binaries from CRAN so I think I am using the standard Rblas.dll.
To reproduce what I see you must run the code at the beginning of the
R session. After the second run, all subsequent runs give the same
result as the second
Yes. Try something akin to
fm1- lme(y~time, data, random=~time|ID)
fm2-update(fm1, correlation=corAR1(form~time|ID)
You can then use anova(fm1,fm2) to compare.
Harold
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of Chris Solomon
Sent: Tuesday,
On Tue, 14 Sep 2004, Francisco Chamu wrote:
I have run this on both Windows 2000 and XP. All I did was install
the binaries from CRAN so I think I am using the standard Rblas.dll.
To reproduce what I see you must run the code at the beginning of the
R session.
We did, as you said `start
On Tue, 14 Sep 2004, Chris Solomon wrote:
Hello-
I'm trying to do some repeated measures ANOVAs. In the past, using SAS,
I have used the framework outlined in Littell et al.'s SAS System for
Mixed Models, using the REPEATED statement in PROC MIXED to model
variation across time within an
Hello Chris,
as far as I know from the Pinheiro and Bates book Mixed-Effects Models in S
and S-PLUS, the autoregressive parameter has to be specified only as
initial value for the estimation. That is, the parameter will be estimated,
but the result may depend on the prespecified value. You do not
Hi all,
I was able to replicate Francisco's observation. I'm using R-1.9.1
installed from binaries on Windows 2000 Pro.
[Previously saved workspace restored]
library(MASS)
data(painters)
pca.painters - princomp(painters[ ,1:4])
loadings(pca.painters)
Loadings:
Comp.1 Comp.2
Sundar:
As I understand it, you can easily create an index variable (a pointer,
actually) that will pick out the y columns in order:
z-yourdataframe
y-as.vector(z[,indexvar])
So if you could cbind() the x's, you'd be all set.
Again, assuming I understand correctly, the x column you want is:
I can't find a way to control the distance between tick marks and tick
labels
in lattice. In a stripplot I'm making these components are too close.
I don't
see anything like base graphics mgp in the scales list.
Thanks for your help,
Scott Waichler
Pacific Northwest National Laboratory
Ditto here, although not from a fresh session. Also 1.9.1 binary from CRAN,
on WinXPPro:
library(MASS)
data(painters)
pca.painters - princomp(painters[ ,1:4])
loadings(pca.painters)
Loadings:
Comp.1 Comp.2 Comp.3 Comp.4
Composition 0.484 -0.376 0.784 -0.101
Drawing 0.424
FWIW, I see the same behavior as Francisco on my Windows machine (also an
installation of the windows binary without trying to install any special
BLAS libraries):
library(MASS)
data(painters)
pca.painters - princomp(painters[ ,1:4])
loadings(pca.painters)
Loadings:
Comp.1
Try:
x - data.frame(x1 = 1: 5, y11 = 1: 5,
x2 = 6:10, y21 = 6:10, y22 = 11:15,
x3 = 11:15, y31 = 16:20,
x4 = 16:20, y41 = 21:25, y42 = 26:30, y43 = 31:35)
df.names-names(x)
ynames-df.names[grep(y,df.names)]
xnames-substring(sub(y,x,ynames),1,2)
Bert,
Coming up with nvec was what I was missing. Modifying your solution
slightly, here's what I ended up with:
z - x[, nmx]
nvec - seq(length(nmx))
z - unlist(z[, rep(nvec, repy)])
z2 - data.frame(x = z, y = y, row.names = NULL)
Thanks again,
--sundar
Berton Gunter wrote:
Sundar:
As I
Try this:
is.x - substr(colnames(x),1,1) == x # TRUE if col name starts with x
x. - unlist(rep(x[,is.x], diff(which(c(is.x,TRUE)))-1)) # repeat x cols
names(x.) - NULL
y. - unlist(x[,!is.x])
DF - data.frame(x = x., y = y., row.names = NULL)
Sundar Dorai-Raj sundar.dorai-raj at PDF.COM
On Tuesday 14 September 2004 11:44, Waichler, Scott R wrote:
I can't find a way to control the distance between tick marks and tick
labels
in lattice. In a stripplot I'm making these components are too close.
I don't
see anything like base graphics mgp in the scales list.
There's none in
Hi Does R have a proceedure/software for capture recapture? Thank you.
Lawrence Lessner
Please respond to [EMAIL PROTECTED]
__
[EMAIL PROTECTED] mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide!
At www.r-project.org - search - R site search, I just got 11
hits for capture-recapture. Did you try this? They may not all be
relevant, but I suspect that some of them might be.
hope this helps. spencer graves
Lawrence Lessner wrote:
Hi Does R have a proceedure/software for
Lawrence-
Are you familiar with the computer software MARK
(http://www.cnr.colostate.edu/~gwhite/mark/mark.htm)?
It is a very complete package for analysis of capture-recapture data. I
have written an interface to MARK in R but it isn't ready for
distribution. I'm presently in the process of
Here is another variation. It uses LOCF which is last
observation carried forward -- a function which takes a
logical vector and for each element provides the index of
the last TRUE value. The version of LOCF here assumes
that the first element of the argument is TRUE which
happens to be the
I have (still) some memory problems, when trying to allocate a huge array:
WinXP pro, with 2G RAM
I start R by calling:
Rgui.exe --max-mem-size=2Gb (as pointed out in R for windows FAQ)
R.Version(): i386-pc-mingw32, 9.1, 21.6.2004
## and here the problem
x.dim - 46
y.dim - 58
slices - 40
Christoph Lehmann wrote:
I have (still) some memory problems, when trying to allocate a huge array:
WinXP pro, with 2G RAM
I start R by calling:
Rgui.exe --max-mem-size=2Gb (as pointed out in R for windows FAQ)
Not sure that it actually says to use 2Gb there. You might try
Did you read the *rest* of what the rw-FAQ says?
Be aware though that Windows has (in most versions) a maximum amount of
user virtual memory of 2Gb, and parts of this can be reserved by
processes but not used. The version of the memory manager used from R
1.9.0 allocates large objects in
Hello:
I would like to generate rolling, multiperiod forecasts from an
estimated ARIMA model, but the function predict.Arima seems
only to generate forecasts from the last observation in the data
set. To implement this, I was looking for an argument like
'newdata=' in predict.lm.
I can
I have the following problem.
I want to use pairs function and get a matrix of scatterplots with the
correlations in the upper panel and the ordinary scatterplots in the
lower panel.
Moreover, I want to have points colored in five differet ways in the
lower panel, because I have five
On Sep 13, 2004, at 8:47 PM, Tamas K Papp wrote:
It ran smoothly on my installation.
version
_
platform powerpc-apple-darwin6.8
arch powerpc
os darwin6.8
system powerpc, darwin6.8
status
major1
minor9.1
year 2004
month06
day 21
language R
Typical lines
I am an economist who decided it's high time that I learned some
Bayesian statistics. I am following An Introduction to Modern
Bayesian Econometrics by T. Lancaster.
The book recommends using BUGS, but I wonder if there are any
alternatives which are free software and fully integrated to R
I've found
Bayesian Methods: A Social and Behavioral Sciences Approach
by Jeff Gill
useful as an introduction. The examples are written in R and S with generalized scripts for doing
a variety of problems. (Though I never got change-point analysis to successfully in R.)
Best, Mark Hall
I have just realised that I sent this to Per only. For those interested on
the list:
-Original Message-
From: Gygax Lorenz FAT
Sent: Tuesday, September 14, 2004 4:35 PM
To: 'Per Toräng'
Subject: RE: [R] glmmPQL and random factors
Hi Per,
Hi,
I am wondering if a Monte Carlo method (or equivalent) exist permitting to test the
randomness of a cluster analysis (eg got by
hclust(). I went through the package fpc (maybe too superficially) but dit not find
such method.
Thanks for any hint,
Patrick Giraudoux
69 matches
Mail list logo