Hi All
I had to assure the randomness of backgammon game's two dices.
I took 8,000,000 observations and found that one of the dices is not uniform
distributed.
I took 80,000,000 observations and found that the dice is uniform
distributed.
In both tests I used Pearson Chi Square
On Mon, 15 May 2006, Satoshi Takahama wrote:
Hello everyone,
If I want to convert or combine a (large) character matrix into a data
frame without having any of its columns convert into a factor class,
is there a simple solution? I() says it will operate on 'an
object' but it seems that
Nancy Lo wrote:
Hi,
I am using a note book computer to do simulation. Even I added a commend
line in the shortcut on my desktop for R. I still get a message saying
reached total memory allocation of 502 Mb. Any suggestions will be welcome.
Nancy
Please specify more details, as the
For those of you who want to learn how to integrate AD Model Builder with R,
and are going to the JSM in Seattle, the following course may be of interest:
Statistical and ecological modeling in AD Model Builder
Date: August 5, 2006 (the day before the Joint Statistical Meetings open)
For more
Dear All:
I tried to fit negative binomial distribution to data in terms of mean and
mean is also a quadratic function of another variable. The likelihood function
is:
function (parameters, y1,x11)
{
p-parameters[1]
alpha1-parameters[1]
beta1-parameters[2]
delta1-parameters[3]
Rogerio Porto wrote:
Hello all.
I have a large .txt file whose variables are fixed-columns,
ie, variable V1 goes from columns 1 to 7, V2 from 8 to 23 etc.
This is a 60GB file with 90 variables and 60 million observations.
I'm working with a Pentium 4, 1GB RAM, Windows XP Pro.
I tried
Nameeta Lobo wrote:
Hello all
thank you very much for all your suggestions. I actually need binary
representations. I tried all the methods that Marc,Jim and Charles have
suggested and they ran fine(thanks a lot). I tried doing it then with 26 and
13
and that's when the computer gave
Hello,
I am new to R and want to use it for some statistical tests. Before the
tests, I want to visualize my data with pairwise scatterplots using
pairs(myData). To analyse it by eye, I would like to draw the function
f(x)=y into each of the scatterplot panels. Unfortunately I get lost
when I
Hello,
I am new to R and want to use it for some statistical tests. Before the
tests, I want to visualize my data with pairwise scatterplots using
pairs(myData). To analyse it by eye, I would like to draw the function
f(x)=y into each of the scatterplot panels. Unfortunately I get lost
when I
rpart allows user-specified splits, so yes, it can.
On Mon, 15 May 2006, Mazerolle, Dan wrote:
Hi all,
I am using the rpart library to construct wildlife distribution models.
As many of the sites that we have sampled have been sampled over many
years, we would like to evaluate binomial
Harald,
for your specific problem (add main diagonal) for example:
abline.points - function (x,y,...) { points(x, y, ...); abline(0,1) }
pairs(USJudgeRatings[1:5], panel = abline.points, cex = 1.5 )
example(pairs) has nicer examples.
Regards
Michael
Am Tuesday 16 May 2006 10:46 schrieb
Dear all,
Is there anybody who can help me to avoid scientific number in the summary of
an lm model?
Here there is an example of a usual output of the lm model.
Thank you!
Guillermo
Example
summary(lm(promiscuity.Index~allK))
Call:
lm(formula = promiscuity.Index ~ allK)
Residuals:
Min
Manuel López-Ibáñez wrote:
Jan T. Kim wrote:
That's an idea I like very much too -- much better than the currently
popular idea of protecting users from the unfriendliness of
programming, anyway...
It is just my opinion that the amount of mail in R-help speaks volumes
about the current
Hello!
I want to calculate the intra- and interrater reliability of my study. The
design is very simple, 5 raters rated a diagnostic score 3 times for 19
patients.
Are there methods/funtions in R? I only found packages to calculate interrater
variability and intraclass correlation
jivan parab wrote:
can u please give me the code written in matlab for partial least square
regression
This is an email list about R, not about Matlab.
will you pls provide me the code written in matlab and also the exlanatiion
to each step
You not only want someone to do the work for you,
From: Erin Hodgess [EMAIL PROTECTED]
To: r-help@stat.math.ethz.ch
Date: Mon, 15 May 2006 05:49:11 -0500
Subject: [R] Rorca zip?
Dear R People:
Is there a Windows binary for ROrca, please?
When I did the search on the R site, there were several links
to a zip file, but all of the links
hi
does anyone have knowledge on how to connected to StatConnector from PHP
I have tried:
new COM(STATCONNECTORSRVLib.StatConnector) or die(Unable to
instantiate StatConnector);
with no succes
regards soren
__
R-help@stat.math.ethz.ch mailing list
Hi Karl,
this sounds like a variance components problem - I suspect that you
will find that you can tackle it using the lme() or lmer() functions
of the nlme and lme4 packages respectively.
The key reference for this tool is:
@Book{PB00:mixed,
author = {J. C. Pinheiro and D. M. Bates},
?options
and see the scipen argument.
For example:
ctl - c(4.17,5.58,5.18,6.11,4.50,4.61,5.17,4.53,5.33,5.14)
trt - c(4.81,4.17,4.41,3.59,5.87,3.83,6.03,4.89,4.32,4.69)
group - gl(2,10,20, labels=c(Ctl,Trt))
weight - c(ctl, trt)
lm.D9 - lm(weight ~ group)
summary(lm.D9)
options(scipen=20)
It sounds as thought you are interested in Hoyt's Anova which is a form
of generalizability theory. This is usually estimated using by getting
the variance components from ANOVA.
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of Karl Knoblick
Sent:
You can read fixed-width-files with read.fwf(). But my rough calculation says
that your dataset will require 40GB of RAM. I don't think you'll be able to
read the entire thing into R. Maybe look at a subset?
-roger
Rogerio Porto wrote:
Hello all.
I have a large .txt file whose variables
On 5/16/2006 5:46 AM, Joerg van den Hoff wrote:
Manuel López-Ibáñez wrote:
Jan T. Kim wrote:
That's an idea I like very much too -- much better than the currently
popular idea of protecting users from the unfriendliness of
programming, anyway...
It is just my opinion that the amount of
Dear R users,
I would like to output the results from the image() function to a specific
grid viewport (to be able to create an .svg afterwards using the gridToSVG()
function).
Is there a function comparable to the image() function but specifically for
the grid system?
The trial code I used
Dear Andrew and Karl,
See the irr, condord, and psy packages.
Regards,
John
On Tue, 16 May 2006 20:37:00 +1000
Andrew Robinson [EMAIL PROTECTED] wrote:
Hi Karl,
this sounds like a variance components problem - I suspect that you
will find that you can tackle it using the lme() or lmer()
Try to open your db with MySQL and use RMySQL
- Message d'origine
De : Roger D. Peng [EMAIL PROTECTED]
Ã⬠: Rogerio Porto [EMAIL PROTECTED]
Cc : r-help@stat.math.ethz.ch
Envoyé le : Mardi, 16 Mai 2006, 1h55mn 41s
Objetà : Re: [R] Large database help
You can read fixed-width-files
Dear R Users
I have difficulties accessing the variance components for an lme fit when the
variance covariance matrix of the random effects is not positive definite.
Can anyone inform me on how to get by this ?
Thanks in advance
Pryseley
I wish to reproduce a locally-weighted fit from (old-fashioned) lowess with
(new-fangled) loess.
Although ?lowess advises that loess has different defaults, does anyone know
how to set up loess so it does exactly the same as lowess did?
Fiddling about suggests setting loess(x, span = 0.8,
On Tue, 2006-05-16 at 09:45 +0200, Uwe Ligges wrote:
Nameeta Lobo wrote:
Hello all
thank you very much for all your suggestions. I actually need binary
representations. I tried all the methods that Marc,Jim and Charles have
suggested and they ran fine(thanks a lot). I tried doing it
Hy all,
I wish to use a variable as rownames for a set of datas.
By example :
nom-toto
prenom-tutu
res-c(1,2)
res-t(res)
res
[,1] [,2]
[1,]12
colnames(res)-c(nom,prenom)
res
toto tutu
[1,]12
nom
[1] toto
I wish to call the rowname
Xin jasonshi510 at hotmail.com writes:
Dear All:
Then error messga there: initial value in 'vmmin' is not finite
In addition: There were 38 warnings (use warnings() to see them).
Could you give some advice please?
Thanks a lot!
Xin Shi
[[alternative HTML version
On Tue, 16 May 2006, Stuart Leask wrote:
I wish to reproduce a locally-weighted fit from (old-fashioned) lowess with
(new-fangled) loess.
`new-fangled': from 1991! (That is, well before R.)
Although ?lowess advises that loess has different defaults, does anyone know
how to set up loess so
Dear R users,
I'm trying to create images of 2 symmetrical matrices using heatmap.2 of the
'gplots' library. Both matrices have the same row and column names
For the first matrix I'm using the default clustering and column denrogram
options:
heatmap.2(Matrix1, symm=T, breaks=20,
Dear Folks:
My recent question concerning split data.frame has been solved using both
following approaches. (Question is listed at the end.)
solution 1: Try melt from the reshape package (by Gabor)
library(reshape)
melt(DF, 1:2)
You may may need to resort it if the order is important.
Hi all
I am actually trying to create a permuatation matrix whose number of columns
equals the number of subjects I have. This matrix is actually supposed to have
-1s and 1s in a unique combination and hence the expand grid question before.
The reason I need this was because I was going to
Many thanks.
No, I confess I didn't look up the loess ref in Statistical Models in S.
Since loess has been around for so long, I assumed that someone else had
encountered this challenge, and found a solution, some time ago.
In the general case this leads to an informative posting in the
Hi,
has anybody an example of an Engel curve analysis in R
or does there exist a package to estimate and plot
Engel curves from expenditure / income data in R?
Thanks a million for your hints,
Werner
__
R-help@stat.math.ethz.ch mailing list
RSiteSearch('Engel')
This search capability can also be used directly via a browser in CRAN.
-- Bert Gunter
Genentech Non-Clinical Statistics
South San Francisco, CA
The business of the statistician is to catalyze the scientific learning
process. - George E. P. Box
-Original
Thks for your answer,
I wish to access colnames by using vars like in my example nom=toto , if i
type toto in R it will display the column named toto with the values that the
column keep, i wish the same by typing nom, but in my case it shows the content
of the var nom not the values of the
If you are SURE you need the exact distribution of the t-statistic (and
asymptotics ala Pitman, 1938, just will not do), then
You can get it by hacking the
dec.index - dec.index + something.else
line in the for loop and making several passes through the data to
accumulate the
On May 16, 2006, at 8:15 AM, justin bem wrote:
Try to open your db with MySQL and use RMySQL
I've seen this offered up as a suggestion a few times but with little
detail. In my experience, even using SQL to pull in data from a
MySQL DB, R would need to load the entire data set into RAM
On Tue, 16 May 2006, Robert Citek wrote:
On May 16, 2006, at 8:15 AM, justin bem wrote:
Try to open your db with MySQL and use RMySQL
I've seen this offered up as a suggestion a few times but with little
detail. In my experience, even using SQL to pull in data from a
MySQL DB, R would
This isn't an R question is it?
PLEASE do read the posting guide!
http://www.R-project.org/posting-guide.html
But I know why this happened, so I will answer your question anyway.
During the first few months of rolling your die, it really was not
uniform and this is what accounts for
Hi
I was confused by your example, too and still am a little bit
confused, so please excuse me if my example does not answer your
question:
You have a data.frame dat
dat - data.frame(colname1 = 1:4, colname2 = 5:8)
You have a vector with the colnames of dat
namen - colnames(dat)
You want to
help
Hello,
I´m a statistic student in Austria and I have to do a survival analysis in R
by using psplines as regressor. My problem is that I sometimes (I think it
depends on the choose of the parameters) get a error message, but I do not
know what it means. After that I tried the procedure with
Hello,
Thanks for the hint.
grid.segments seemed the closest I got.
I did manage to draw (well fake) a line with it. I can only address the
whole drawing frame, which means I can only adjust the position and length
of the line
by trial and error. I see no way to address the y axis scale of my
To whom it may concern,
I seem to be encountering RAM problems using R. I am currently running R
2.2.0 in Windows XP version 2002. My computer has 4 GB of RAM. The
problem I am encountering is that R does not seem to allow me to obtain
more than 1.5GB of RAM as I construct various
On Tue, 16 May 2006, Birgit Rauchenschwandtner wrote:
help
Hello,
I´m a statistic student in Austria and I have to do a survival analysis in R
by using psplines as regressor. My problem is that I sometimes (I think it
depends on the choose of the parameters) get a error message, but I do not
Need to compute the numerical value of the solution of the quadratic form:
X'AX=I
where all the matrices are (n*n) [n may be very large]; X is the matrix
containing n*n unknowns (those I need) and A contains the coefficients
of the problem. I is the identity matrix of order n.
thanks in advance,
Hi folks!
I'm trying to update a model using:
form - substitute(.~.-bad.var, list(bad.var=bad.var))
model.update-update(model, form)
where bad.var is a character string like lag(X1, 10)
The problem is that substitute puts lag(X1, 10) in quotes - something I
don't want it to do.
On Tue, 16 May 2006, Kerpel, John wrote:
I'm trying to update a model using:
form - substitute(.~.-bad.var, list(bad.var=bad.var))
model.update-update(model, form)
where bad.var is a character string like lag(X1, 10)
The problem is that substitute puts lag(X1, 10) in quotes - something I
I'm trying to label axes using 'expression', but I would also like to change
the font to Times New Roman (a requirement of the journal). I'd appreciate any
advice.
Here's one example of several things that I've tried, including changes to
par() parameters 'family', 'font', 'font.axis', and
Jean-Pierre Dube wrote:
To whom it may concern,
I seem to be encountering RAM problems using R. I am currently running R
2.2.0 in Windows XP version 2002. My computer has 4 GB of RAM. The
problem I am encountering is that R does not seem to allow me to obtain
more than 1.5GB of RAM
Look at
?eigen
In your example, you are looking for
eigen(A)
__
R-help@stat.math.ethz.ch 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 everyone,
I have a large dataset (x) with some rows that have duplicate variables
that I would like to remove. I find which rows are the duplicates with
X1-which(duplicated(x)). That gives me the rows with duplicated
variables. Now, how can I remove just those rose from the original data
On Tue, 16 May 2006 17:42:22 +0100 Viktor Tron wrote:
Hello,
Thanks for the hint.
grid.segments seemed the closest I got.
I did manage to draw (well fake) a line with it. I can only address
the whole drawing frame, which means I can only adjust the position
and length of the line
by trial
On Tue, 2006-05-16 at 14:37 -0400, Guenther, Cameron wrote:
Hello everyone,
I have a large dataset (x) with some rows that have duplicate variables
that I would like to remove. I find which rows are the duplicates with
X1-which(duplicated(x)). That gives me the rows with duplicated
Marc,
I have tried unique but unique looks at the entire row. I have a data
set with a variable TRIPID. The dataset has 469,000 rows. In most
cases TRIPID is a unique value. However, in some cases I have the same
TRIPID value but different values for other variables. What this
amounts to is
Thanks Phil
That worked pergectly.
Cameron Guenther, Ph.D.
Associate Research Scientist
FWC/FWRI, Marine Fisheries Research
100 8th Avenue S.E.
St. Petersburg, FL 33701
(727)896-8626 Ext. 4305
[EMAIL PROTECTED]
-Original Message-
From: Phil Spector [mailto:[EMAIL PROTECTED]
Sent:
Dear R-users,
I am working with some grid-search optimization for 13 values of an
object function.
At glance one may compute the object function for each possible
combination of set of parameters, but in my case would not be feasible:
taking for example 13 possible values for each parameters
On 5/16/06, Pryseley Assam [EMAIL PROTECTED] wrote:
Dear R Users
I have difficulties accessing the variance components for an lme fit when
the variance covariance matrix of the random effects is not positive definite.
Well, that shouldn't happen. A variance-covariance matrix is, by
--- percy tiglao [EMAIL PROTECTED] wrote:
Thanks alot :-)
I'm already reading through it.
Yes. A very useful piece. Thanks
On 5/15/06, Owen, Jason [EMAIL PROTECTED] wrote:
R Community:
In case you haven't seen it, I have written a
manuscript
called The R Guide that is
On Tue, 2006-05-16 at 14:54 -0400, Guenther, Cameron wrote:
Marc,
I have tried unique but unique looks at the entire row.
Right, as I noted in the last line of my reply.
I have a data
set with a variable TRIPID. The dataset has 469,000 rows. In most
cases TRIPID is a unique value.
Hi
Scott Rollins wrote:
I'm trying to label axes using 'expression', but I would also like to change
the font to Times New Roman (a requirement of the journal). I'd appreciate
any advice.
Here's one example of several things that I've tried, including changes to
par() parameters
I am trying to apply nmds to a data matrix but I receive the
following error message:
Error in isoMDS(dis, y = y, k = k, maxit = maxit) :
zero or negative distance between objects 5 and 7
The data are in a vegetation cover-class matrix (species in columns,
plots in rows, classes 1-8
The irr package seems to install correctly:
install.packages(irr)
trying URL 'http://cran.us.r-project.org/src/contrib/irr_0.61.tar.gz'
Content type 'application/x-tar' length 13848 bytes
opened URL
==
downloaded 13Kb
* Installing *source* package
You have not told us the device you are using, nor your version of R.
The way to select 'Times New Roman' is via the family, which will work if
it is supported on the device and your R is current (2.3.0, preferably
patched).
font=6 is only supported on a few devices, and not on X11 on Linux.
Duncan Murdoch wrote (amongst other things):
Statistical computing is not easy, so how could R be? Who has ever
claimed it is? Any package that makes statistical computing appear to
be easy is probably giving you wrong answers half the time, or is
extremely limited in scope.
Rick -
You didn't say what operating system you're using, but I suppose
that it's UNIX-based, because you're using the tar.gz files.
If this is the case, the problem can be resolved by setting the
environmental variable LANG to an empty string. If you're currently
just typing R at the shell
On May 16, 2006, at 11:19 AM, Prof Brian Ripley wrote:
Well, there *is* a manual about R Data Import/Export, and this does
discuss using R with DBMSs with examples. How about reading it?
Thanks for the pointer:
http://cran.r-project.org/doc/manuals/R-data.html#Relational-
databases
On Tue, 2006-05-16 at 13:51 -0700, Phil Spector wrote:
Rick -
You didn't say what operating system you're using, but I suppose
that it's UNIX-based, because you're using the tar.gz files.
If this is the case, the problem can be resolved by setting the
environmental variable LANG to an
In ancient times, 1999 or so, Alvaro Novo and I experimented with an
interface to mysql that brought chunks of data into R and accumulated
results.
This is still described and available on the web in its original form at
http://www.econ.uiuc.edu/~roger/research/rq/LM.html
Despite
Rolf Turner wrote:
Duncan Murdoch wrote (amongst other things):
Statistical computing is not easy, so how could R be? Who has ever
claimed it is? Any package that makes statistical computing appear to
be easy is probably giving you wrong answers half the time, or is
extremely limited in
On Tue, 16 May 2006, roger koenker wrote:
In ancient times, 1999 or so, Alvaro Novo and I experimented with an
interface to mysql that brought chunks of data into R and accumulated
results.
This is still described and available on the web in its original form at
Can the packages that I've installed in R 2.2.1 be automatically imported to
R 2.3.0?
Otherwise it is a hassle everytime I update a new version of R...
[[alternative HTML version deleted]]
__
R-help@stat.math.ethz.ch mailing list
On Tue, 2006-05-16 at 17:49 -0300, Rolf Turner wrote:
Duncan Murdoch wrote (amongst other things):
Statistical computing is not easy, so how could R be? Who has ever
claimed it is? Any package that makes statistical computing appear to
be easy is probably giving you wrong answers half
On Tue, 16 May 2006, Rick Bilonick wrote:
On Tue, 2006-05-16 at 13:51 -0700, Phil Spector wrote:
Rick -
You didn't say what operating system you're using, but I suppose
that it's UNIX-based, because you're using the tar.gz files.
If this is the case, the problem can be resolved by
On Tue, 16 May 2006 16:50:36 -0500 Marc Schwartz (via MN) wrote:
On Tue, 2006-05-16 at 17:49 -0300, Rolf Turner wrote:
Duncan Murdoch wrote (amongst other things):
Statistical computing is not easy, so how could R be? Who has
ever claimed it is? Any package that makes statistical
I am fairly new to R so I don;t know much about manipulating it yet.
I would like to make some sparklines http://www.edwardtufte.com/
bboard/q-and-a-fetch-msg?msg_id=0001ORtopic_id=1 using R. I have th
pre-compiled R for MacOS on intel. Any help would be greatly
appreciated.
Thanks
Robert
Press ESC is supposed to be able to stop the current computation... but
often times it does not really work... I had to kill the R process to stop
and restart the session...
What is a good way to stop/cancel the current computation?
[[alternative HTML version deleted]]
I think you did not answer my question... I now upgraded my main R program
from 2.2.1 to 2.3.0 and I removed the 2.2.1 installation, but all the
libraries that I've installed myself previously in 2.2.1 are now gone... I
have to reinstall them one by one... in 2.3.0
On 5/16/06, Berton Gunter
I am using Windows XP... Ctrl-C or Break does not stop the computation in
R...
On 5/16/06, Berton Gunter [EMAIL PROTECTED] wrote:
This is usually an OS issue (you did not tell us what OS you are on), but
it's usually ctrl-C or Break not esc that does it, I think.
-- Bert
-Original
Well,
if you removed your installed packages together with R then of course you
need to reinstall them consider installing the packages in a separate
directory, not in the R tree. (I install packages into ~/.R/library on my
Linux box for example.) This way at the next R update you don't need
Also this may depend on the package you're using. For some functions in some
packages ctrl+c (or anything else) does not work. This is because C code
cannot be interrupted unless it is written with some care, and there are
some packages which do not support interruption. At least my igraph
I think you did not answer my question... I now upgraded my
main R program
from 2.2.1 to 2.3.0 and I removed the 2.2.1 installation, but all the
Wait until after you use update.packages() to remove your previous
installation. You can keep multiple versions of R simultaneously, so this is
no
I'm not familiar with Pairwise likelihood. RSiteSearch(pairwise
likelihood) produced 22 hits, none of which looked to me like they were
relevant, though some may be. Google returned some articles that might
be relevant, but that was not clear to me.
However, some
The process that you outlined is described in the Windows R FAQ 2.8 What's
the best way to upgrade?.
The same process should apply to other platforms as well, but I don't think
it is in the general R FAQ. Perhaps it should.
-Christos
-Original Message-
From: [EMAIL PROTECTED]
Thanks loads, extremely useful.
So just for the record again.
You have to get into the viewport that contains the actual graphplot with
the axes.
So, how you do it:
1. viewports have names.
In case of histograms the relevant viewport is something like
plotA.panel.X.Y.vp
where A is the plot
an upgrade: from the flintstones -- to the michelin man...
On May 16, 2006, at 4:40 PM, Thomas Lumley wrote:
On Tue, 16 May 2006, roger koenker wrote:
In ancient times, 1999 or so, Alvaro Novo and I experimented with an
interface to mysql that brought chunks of data into R and
two ARIMA models, both have several bars signicant in ACF and PACF plots of
their residuals,
but when run Ljung Box tests,
both don't show any significant correlations...
however, one model has p-value that is larger than the other model,
based on the p-values,
can I say the model with larger
If first fit my data column V1 to column V2 using normal lm fitting,
call it fit1,
then I used acf(fit1$residuals, type='cov', 40) function to obtain the
autocovariance of the residuals,
and then constructed a autocovariance matrix, I chose it to be 40x40.
Call this autocovariance matrix B,
Last chance to register for Sydney, Australia Data Mining Training and
Workshop.
Training: June 5-7, 2006
Workshop: June 8, 2006
Register for 1, 2, 3 or 4 days.
Location:
Australian Graduate School of Management, Lecture Theatre
Number 1 O'Connell Street
Sydney, Australia
Objective: The
Thanks! But I could not find a corresponding example in the mentioned book - at
least the calculation of the intra- and interrater variability or the
corresponding intraclass coefficients.
Karl
- Ursprüngliche Mail
Von: Andrew Robinson [EMAIL PROTECTED]
An: Karl Knoblick [EMAIL
Daer John,
Thanks! My problem is that in the mentioned packages (irr, condord, and psy)
only a n*m matrix for subjects and raters can be given as data. So I can not
give the repetitions as data and therefore it will be impossible to get results
for the INTRArater reliability.
Regards,
Karl
93 matches
Mail list logo