What's the reommended way, in R, to determine the strength of
association between a discrete variable and a continuous variable?
Yes, I have read the manuals, trawled the archives, c.
__
[EMAIL PROTECTED] mailing list
Gabor Grothendieck wrote:
Duncan Murdoch dmurdoch at pair.com writes:
: On Sun, 11 Jul 2004 21:28:44 + (UTC), Gabor Grothendieck
: ggrothendieck at myway.com wrote:
:
: 1. when I run skeleton.package realize that I must use the arg
:path = library
:The example that is shown there
On 07/12/04 18:23, Richard A. O'Keefe wrote:
What's the reommended way, in R, to determine the strength of
association between a discrete variable and a continuous variable?
Analysis of variance? aov()? R^2?
It seems to me, though, that there are many possible answers, and
this really doesn't
How does one implement a likelihood-ratio test, to test whether the
variances of the random effects differ between two groups of subjects?
Suppose your data consist of repeated measures on subjects belonging to
two groups, say boys and girls, and you are fitting a linear mixed-effects
model for
I'm wondering if mutual information al la Cover Thomas (1991, Ch 2) is
not the killer association measure for all types of random variables?
Murray Jorgensen
PS Yes, this is probably OT!
Richard A. O'Keefe wrote:
What's the reommended way, in R, to determine the strength of
association between
I had an error message while running a macro from Yudi Pawitan's web site:
source(ex2-13.r)
Error in parse(file, n, text, prompt) : syntax error on line 2
Inspecting ex2-13.r I found that the error was generated by another
source() command.
Clearly R does not like nested source()s, which is
Hi
I use pixmapIndexed
tmp.vimp - array(0,c(x.dim,y.dim))
tmp.vimp - pixmapIndexed(tmp.vimp, col=rainbow)
to plot values of a 2D matrix. I 'fill' the pixmapIndexed like:
for (x in 1:x.dim) {
for (y in 1:y.dim) {
[EMAIL PROTECTED],y] - my.matrix[x,y]
}}
how can I
tmp.vimp - matrix(NA, nrow = x.dim, ncol = y.dim)
tmp.vimp - image(tmp.vimp, col=rainbow)
gives:
Error in image.default(tmp.vimp, col = rainbow) :
invalid z limits
In addition: Warning messages:
1: no finite arguments to min; returning Inf
2: no finite arguments to max; returning -Inf
even
Uwe Ligges ligges at statistik.uni-dortmund.de writes:
:
: Gabor Grothendieck wrote:
:
: Duncan Murdoch dmurdoch at pair.com writes:
:
: : On Sun, 11 Jul 2004 21:28:44 + (UTC), Gabor Grothendieck
: : ggrothendieck at myway.com wrote:
: :
:
: : 1. when I run skeleton.package
Murray Jorgensen wrote:
I had an error message while running a macro from Yudi Pawitan's web site:
source(ex2-13.r)
Error in parse(file, n, text, prompt) : syntax error on line 2
Inspecting ex2-13.r I found that the error was generated by another
source() command.
Clearly R does
It can't calculate the range of the x/y axes so there's nothing to plot.
-roger
Christoph Lehmann wrote:
tmp.vimp - matrix(NA, nrow = x.dim, ncol = y.dim)
tmp.vimp - image(tmp.vimp, col=rainbow)
gives:
Error in image.default(tmp.vimp, col = rainbow) :
invalid z limits
In addition: Warning
Yes, I have source()-ed recursively without problems before. Try
sourcing the second script on its own and see what errors it has.
I have download lkpack.zip
(http://www.mep.ki.se/~yudpaw/likelihood/lkpack.zip)
and source(EX2-13.R) and source(li.r) which is called without
problem. The timestamp
Gabor Grothendieck wrote:
Uwe Ligges ligges at statistik.uni-dortmund.de writes:
:
: Gabor Grothendieck wrote:
:
: Duncan Murdoch dmurdoch at pair.com writes:
:
: : On Sun, 11 Jul 2004 21:28:44 + (UTC), Gabor Grothendieck
: : ggrothendieck at myway.com wrote:
: :
:
: : 1. when I
__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
From: Gabor Grothendieck
[snip]
The objective should be that creating a package is as easy as this:
f - function()1; g - function()2; d - 3; e - 4:5
package.skeleton(list=c(f,g,d,e), name=AnExample)
library(AnExample)
f()
which means that the package needs to be
On Mon, 12 Jul 2004 04:32:55 + (UTC), Gabor Grothendieck
[EMAIL PROTECTED] wrote :
Duncan Murdoch dmurdoch at pair.com writes:
: On Sun, 11 Jul 2004 21:28:44 + (UTC), Gabor Grothendieck
: ggrothendieck at myway.com wrote:
:
: 1. when I run skeleton.package realize that I must use the
On Mon, 12 Jul 2004 15:14:27 +0200, Uwe Ligges
[EMAIL PROTECTED] wrote :
Gabor Grothendieck wrote:
The objective should be that creating a package is as easy as this:
f - function()1; g - function()2; d - 3; e - 4:5
package.skeleton(list=c(f,g,d,e), name=AnExample)
Hello All:
Suppose the following little data frame:
x - data.frame(dog = c(3,4,6,2,8), cat = c(8,2,3,6,1))
x$cat
[1] 8 2 3 6 1
How can I get the paste() function to do the same thing. The command
below is obviously wrong
paste(x, cat, sep = $)
Thanks,
ANDREW
Do you mean something like:
paste(x$cat, collapse= )
[1] 8 2 3 6 1
??
Andy
From: Andrew Criswell
Hello All:
Suppose the following little data frame:
x - data.frame(dog = c(3,4,6,2,8), cat = c(8,2,3,6,1))
x$cat
[1] 8 2 3 6 1
How can I get the paste() function to do the
On Mon, 12 Jul 2004 13:16:06 +0700, Andrew Criswell
[EMAIL PROTECTED] wrote :
Hello All:
Suppose the following little data frame:
x - data.frame(dog = c(3,4,6,2,8), cat = c(8,2,3,6,1))
x$cat
[1] 8 2 3 6 1
How can I get the paste() function to do the same thing. The command
below is
Andrew Criswell wrote:
Suppose the following little data frame:
x - data.frame(dog = c(3,4,6,2,8), cat = c(8,2,3,6,1))
x$cat
[1] 8 2 3 6 1
How can I get the paste() function to do the same thing. The command
^^
On Mon, 2004-07-12 at 01:16, Andrew Criswell wrote:
Hello All:
Suppose the following little data frame:
x - data.frame(dog = c(3,4,6,2,8), cat = c(8,2,3,6,1))
x$cat
[1] 8 2 3 6 1
How can I get the paste() function to do the same thing. The command
below is obviously wrong
Now that I see Duncan's reply, I believe that's what Andrew wanted. I
really should read messages more carefully...
Andy
From: Liaw, Andy
Do you mean something like:
paste(x$cat, collapse= )
[1] 8 2 3 6 1
??
Andy
From: Andrew Criswell
Hello All:
Suppose the
Marc Schwartz wrote:
On Mon, 2004-07-12 at 01:16, Andrew Criswell wrote:
Hello All:
Suppose the following little data frame:
x - data.frame(dog = c(3,4,6,2,8), cat = c(8,2,3,6,1))
x$cat
[1] 8 2 3 6 1
How can I get the paste() function to do the same thing. The command
below is obviously wrong
Hello,
My name is Matthias and I do look for syntax regarding hierarchal anova
models in R. How can I express that a factor is nested within the
combination of two other factors A(B,C), e.g. for aov(...)? I did not find
the corresponding expression. Furthermore, I wanted to ask whether block
Hi!
I have a grid based system where I can only install R via RPM. Is there a
working RPM spec file that I can use to create RPM package? Or is there
already RPMs for different types of platform?
Thanks!
Soichi H.
**
Hi,
How do I open an excel file in R? I have save the excel file in unicode text
format, but it is not possible to open the file in R.
/Pernilla
__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
PLEASE do read the
-Ursprüngliche Nachricht-
Von: TEMPL Matthias
Gesendet: Montag, 12. Juli 2004 17:29
An: 'Pernilla Karlsson'
Betreff: Re: [R] Excel file
It is possible to open the file in R.
?read.table
excelfile - read.table(../file.txt, sep=,, header=T)
when you have a komma as
On Mon, 12 Jul 2004 17:17:00 +0200, Pernilla Karlsson
[EMAIL PROTECTED] wrote :
Hi,
How do I open an excel file in R? I have save the excel file in unicode text
format, but it is not possible to open the file in R.
See the R Data Import/Export manual (via Help|Manuals). You probably
want to
Pernilla Karlsson wrote:
Hi,
How do I open an excel file in R? I have save the excel file in unicode text
format, but it is not possible to open the file in R.
/Pernilla
__
[EMAIL PROTECTED] mailing list
Hello,
For a simple problem with 1 predictor (x) and 2 classes (0 and 1), the
linear discriminant function should be something like
2(mu_0 - mu_1)/var x+x-independent-terms
Dear R users
this may be a simple question - but i would appreciate any thoughts
does anyone know how you would get one lower and one upper confidence
interval for a set of data that consists of proportions. i.e. taking a
usual confidence interval for normal data would result in the lower
Christoph Lehmann wrote:
Hi
I use pixmapIndexed
tmp.vimp - array(0,c(x.dim,y.dim))
tmp.vimp - pixmapIndexed(tmp.vimp, col=rainbow)
to plot values of a 2D matrix. I 'fill' the pixmapIndexed like:
for (x in 1:x.dim) {
for (y in 1:y.dim) {
[EMAIL PROTECTED],y] -
From: Paul Roebuck
On Sun, 11 Jul 2004, Gabor Grothendieck wrote:
search.path -
function(fn,
paths = strsplit(Sys.getenv(PATH), split = ;)[[1]],
fsep = \\) {
for(d in paths) {
f - file.path(d, fn, fsep = fsep)
if (file.exists(f))
On Mon, 12 Jul 2004 11:01:49 -0500 (CDT), Paul Roebuck
[EMAIL PROTECTED] wrote :
I glanced this and thought this might be handy to keep for
possible use. To make it less Windows-specific, I was going
to replace Gabor's fsep default value with '.Platform$file.sep'
when I noticed that .Platform
On Mon, 12 Jul 2004 12:50:30 -0400, Duncan Murdoch [EMAIL PROTECTED]
wrote :
I'm not sure I follow. What would be the difference between
.Platform$file.sep and .Platform$path.sep ?
Stupid me: path.sep would be the separator between paths in the PATH
variable.
Duncan Murdoch
You should consider prop.test or binom.test.
The problem you will find is that these functions are
not intended to do what you want but to give you one
confidence interval at a time.
However a starting point could be :
I-sample(1:50) #the numerator
N-sample(50:200, 50) #the denominator
Have you tried something like the following:
mylmeMF0 - lme( y ~ time, data=DAT, random=~time | id, method=ML)
mylmeMF1 - lme( y ~ time*sex, data=DAT, random=~time | id, method=ML)
anova(mylemFM0, mylmeMF1)
Please check Pinheiro and Bates
My recommendations would depend on the nature of the discrete and
continuous variables. What is the application? Is one variable random
and the other a control variable that might be used to predict? If the
continuous variable is random, what about using the R^2 produced by
lm? For
Hi Pernilla
Under menu go to help manuals R Data Import/Export and read the section
on spreadsheet-like data import.
Also you can try ?read.table
Francisco
From: Pernilla Karlsson [EMAIL PROTECTED]
To: [EMAIL PROTECTED]
Subject: [R] Excel file
Date: Mon, 12 Jul 2004 17:17:00 +0200
Hi,
How do
Folks:
Without entering the debate, just a little addendum to Andy's suggestions and
Duncan's and Gabor's comments about the difficulty of building (simple,
containing no binaries and just a few functions) packages for R.
I have for some time now built such R-code only .rda files. It is a
Dear Mr. Graves:
Thank you. You're right -- I should have said method=ML.
Putting sex into the model as a fixed effect enables us to see if the relationship
between time and y differs between the sexes. But it doesn't do anything to the
covariance of the random effects.
Jake
On Mon, 12 Jul
Please see:
Brown, Cai and DasGupta (2001) Statistical Science, 16: 101-133
and (2002) Annals of Statistics, 30: 160-2001
They show that the actual coverage probability of the standard
approximate confidence intervals for a binomial proportion are quite
poor, while the
Running R scripts via 'R_exe BATCH CMD inpufile outputfile' works fine
with jpeg commands in Windows, but the jpeg commands give an error under
Linux because GUI is set to none. Is there a way to use jpeg commands in
BATCH CMD in Linux.
Thanks, Moises
In The American Statistician vol 42 (1988) pages 257 - 280, Yoav
Benjamini investigates some variations on the box plot, including
vaseplots, which maek the width of each box vary proportionally to he
estimated density at a particular point.
Has anyone implemented these in R ?
Thanks as always
Darren also might consider binconf() in library(Hmisc).
library(Hmisc)
binconf(1, 10, method=all)
PointEstLower Upper
Exact 0.1 0.002528579 0.4450161
Wilson 0.1 0.005129329 0.4041500
Asymptotic 0.1 -0.085938510 0.2859385
binconf(10, 100,
It seems to me that a transformation is in order since [0,1] can't
possibly contain a normal distribution without cutting off both tails.
On Mon, 12 Jul 2004, Rolf Turner wrote:
Darren Shaw wrote:
this may be a simple question - but i would appreciate any thoughts
does anyone know
See ?bitmap.
Andy
From: Moises Hassan
Running R scripts via 'R_exe BATCH CMD inpufile outputfile' works fine
with jpeg commands in Windows, but the jpeg commands give an
error under
Linux because GUI is set to none. Is there a way to use jpeg
commands in
BATCH CMD in Linux.
Thanks,
FWIW, if the exact intervals are what is desired here, as another poster
has already suggested, binom.test() will get you there:
binom.test(1, 10)$conf.int
[1] 0.002528579 0.445016117
attr(,conf.level)
[1] 0.95
binom.test(10, 100)$conf.int
[1] 0.04900469 0.17622260
attr(,conf.level)
[1] 0.95
According to Brown, Cai and DasGupta (cited below), the exact
confidence intervals are hyperconservative, as they are designed to
produce actual coverage probabilities at least the nominal. Thus for a
95% confidence interval, the actual coverage could be 98% or more,
depending on the
Short course: Statistical Learning and Data Mining
Trevor Hastie and Robert Tibshirani, Stanford University
Georgetown University Conference Center
Washington DC
September 20-21, 2004
This two-day course gives a detailed overview of statistical models
for data mining, inference and
There has been a plethora of responses over the past hour or so to a
question posed by Darren Shaw about how to estimate (get a confidence
interval for) a proportion based on a data set consisting of a number
of proportions. These responses have been all off the point. I
would suggest to the
I don't have that article handy, but it sounds like `violin plots' that I've
read (also in Am. Stat.) quite a while ago. I wrote a S-PLUS function to do
something like it back in the dark ages. It probably isn't too hard to
modify for use with R. See if it does something similar:
Hi all,
I'm looking for smooth monotone estimation packages, preferably using splines.
I downloaded the 'cobs' package and intend to use it, but since it offers only
quadratic splines based on L1 minimization, I'd like to compare its performance to
that of a more 'mainstream' cubic-spline,
At 08:45 AM 7/12/2004, marzban wrote:
Hello,
For a simple problem with 1 predictor (x) and 2 classes (0 and 1), the
linear discriminant function should be something like
2(mu_0 - mu_1)/var x+x-independent-terms
where var is the common variance.
Question 1: Why does lda() report only
I believe the `mgcv' package has facility for monotone splines.
HTH,
Andy
From: Assaf P Oron
Hi all,
I'm looking for smooth monotone estimation packages,
preferably using splines.
I downloaded the 'cobs' package and intend to use it, but
since it offers only quadratic splines based
Dear all, there is a R function to construct
simultaneous confidence intervals for multinomial
proportions?
Best regards
jl
__
Acabe com aquelas janelinhas que pulam na sua tela.
AntiPop-up UOL - É grátis!
Perhaps some reading would be helpful. I suggest you look first at the
help file for lda(). Second, I suggest you read Venables and Ripley, MASS,
4th Edition, where lda() is discussed extensively. Third, I suggest you
read Ripley's Pattern Recognition and Neural Networks, where the
Hi Group,
I have a set of number (~300) and I want to find the
unique numbers in it. Basically I want to get rid of
duplicate numbers. Is there a function in R to do
this.
Please help.
thank you
SP
--- F Z [EMAIL PROTECTED] wrote:
Hi Pernilla
Under menu go to help manuals R Data
JEdit has an R syntax highlighting file at http://community.jedit.org. It
is listed as a file pretty far down on the right side menus as R Edit Mode
- Extensive Version.
Ben Stabler
-Original Message-
From: Yi-Xiong Sean Zhou [mailto:[EMAIL PROTECTED]
Sent: Wednesday, July 07, 2004
I implemented Uwe's R-WinEdt package as a closer to true WinEdt mode so that
WinEdt could be used with other modes as well. I sent Uwe the
implementation and we talked a little about cleaning it up, but I haven't
gotten around to it and I don't think he has either. It is still a little
messy as
I haven't done this in years but I think the `scaling' element in the
list returned by lda is the original data matrix multiplied by the
rotation matrix from the SVD. Taking a look at
getAnywhere(lda.default)
will probably answer your question.
-roger
marzban wrote:
Perhaps some reading would
?unique
S Peri wrote:
Hi Group,
I have a set of number (~300) and I want to find the
unique numbers in it. Basically I want to get rid of
duplicate numbers. Is there a function in R to do
this.
Please help.
thank you
SP
--- F Z [EMAIL PROTECTED] wrote:
Hi Pernilla
Under menu go to help
Hi,
Is there a way to use regular expressions to capture two or more words in a
sentence? For example, I wish to to find all the lines that have the words thomas,
perl, and program, such as thomas uses a program called perl, or perl is a
program that thomas uses, etc.
I'm sure this is a
Hi, all;
Basically, the subject contains the all information I need to know.
In e1071 library, there are functions to tune parameters.
They provide several values one of which is the performance.
Does any body know the definition of performance here?
Is it percentage of error or just the error
I'd suggest doing it with multiple regular expressions -- you could
construct a single regular expression for this, but I expect it would get
quite complicated and possibly very slow.
The expression for y in the example below tabulates how many words
matched for each line (i.e., line 2 matched
Paul Roebuck roebuck at odin.mdacc.tmc.edu writes:
On Sun, 11 Jul 2004, Gabor Grothendieck wrote:
search.path -
function(fn,
paths = strsplit(Sys.getenv(PATH), split = ;)[[1]],
fsep = \\) {
for(d in paths) {
f - file.path(d, fn, fsep = fsep)
help.search()
on my machine turns up only:
mono.con(mgcv) Monotonicity constraints for a cubic
regression spline.
smooth.monotone(fda)Monotone Smoothing of Data
pmreg(ftnonpar) Piecewise monotone regression with taut
strings
Looking at the body of tune(), it has:
...
repeat.errors[reps] - if (is.factor(true.y))
1 - classAgreement(table(pred, true.y))
else crossprod(pred - true.y)/length(pred)
...
where classAgreement() is a function defined inside tune() that
Sangick Jeon sijeon at ucdavis.edu writes:
Is there a way to use regular expressions to capture two or more words in a
sentence? For example, I wish to to find all the lines that have the
words thomas,
perl, and program, such as thomas uses a program called perl, or perl
is a
program
Hi everyone,
I'm running R1.9.1 on RedHat Linux. I'm trying to read in a matrix
file with 13956 by 858 dimensions. I realize this is pretty huge, though
I think the amount of memory I have should be able to handle it. R reads
the entire file and tells me Read in 11974247 values. This is
On Mon, 12 Jul 2004 21:53:34 -0400 (EDT), Tianyu Tom Wang
[EMAIL PROTECTED] wrote:
Hi everyone,
I'm running R1.9.1 on RedHat Linux. I'm trying to read in a matrix
file with 13956 by 858 dimensions. I realize this is pretty huge, though
I think the amount of memory I have should be able to
I doubt this is a memory problem, considering that R reported that it
read in the data! What exactly were the commands that you used to read
in the data?
-roger
Tianyu Tom Wang wrote:
Hi everyone,
I'm running R1.9.1 on RedHat Linux. I'm trying to read in a matrix
file with 13956 by 858
I can't get confint.glm to work from within a function. Consider
the following (using R 1.9.1, Windows 2000):
# FIRST: SOMETHING THAT WORKS FROM A COMMAND PROMPT
DF - data.frame(y=.1, N=100)
(fit - glm(y~1, family=binomial, data=DF,
weights=DF[,N]))
Call: glm(formula = y ~ 1,
If you use the RPerl interface
(http://www.omegahat.org/RSPerl/), you can do this simply
for any number of expressions by letting Perl do the
matching. If strings is the array of strings that you
want to test, all you need is
for (@strings) {
if (/thomas/ /perl/ /prog/) {do
__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
76 matches
Mail list logo