wu sz wrote:
> Hello there,
>
> how to do the test for a significant difference between two groups
> mean vectors or across several group means vectors?
This is not an R question.
You cannot test on mean differences without assumptions re. the
variance. Please read a basic textbook on statisti
Hi Doug and Spencer,
Many thanks - Excellent!
All worked out nicely
Regards
Stephen
From: Spencer Graves [mailto:[EMAIL PROTECTED]
Sent: Mon 20/06/2005 17:54
To: Stephen
Cc: r-help@stat.math.ethz.ch
Subject: Re: [R] Mixed model
(comments in line)
On 6/20/05, Weiwei Shi <[EMAIL PROTECTED]> wrote:
> hi,
> i have another question on tapply:
> i have a dataset z like this:
> 5540 389100307391 2600
> 5541 389100307391 2600
> 5542 389100307391 2600
> 5543 389100307391 2600
> 5544 389100307391 2600
> 5546 381300302513
Hi
Barry Rowlingson wrote:
> Pierre Kleiber wrote:
>
>> I want to use contourLines() to get contour line coordinate vectors,
>> but I don't want to make a plot. However contourLines() insists on
>> opening a graphics device. Is there a way tell it not to do this?
>
>
> contourLines() calls
Dear Ben,
To get big factor loadings like this, I'd guess that you have large
communalities (not a bad thing, of course, and probably less rather than
more likely with crudely measured variables) and strong correlations between
factors. I suppose that if the latter get too large that might suggest
On 20 June 2005 at 08:29, bogdan romocea wrote:
| I timed the same code (simulation with for loops) on the same box
| (dual Xeon EM64T, 1.5 Gb RAM) under 3 OSs and was surprised by the
| results:
| Windows XP Pro (32-bit): Time difference of 5.97 mins
| 64-bit GNU/Linux (Fedora Core 4): Ti
On 6/20/05, William H. Asquith <[EMAIL PROTECTED]> wrote:
> I am porting some FORTRAN to R in which an Inf triggers an if(). The
> trigger is infinite on exp(). What is the canonical R style of
> determining OVER when exp(OVER)== Inf? The code structure that I am
> porting is best left intact--s
On 6/20/05, Weiwei Shi <[EMAIL PROTECTED]> wrote:
> hi,
> i have another question on tapply:
> i have a dataset z like this:
> 5540 389100307391 2600
> 5541 389100307391 2600
> 5542 389100307391 2600
> 5543 389100307391 2600
> 5544 389100307391 2600
> 5546 381300302513
On Mon, 2005-06-20 at 18:15 -0500, Weiwei Shi wrote:
> hi,
> i have another question on tapply:
> i have a dataset z like this:
> 5540 389100307391 2600
> 5541 389100307391 2600
> 5542 389100307391 2600
> 5543 389100307391 2600
> 5544 389100307391 2600
> 5546 381300302513
This may help
R>wei
V1 V2 V3
1 5540 389100307391 2600
2 5541 389100307391 2600
3 5542 389100307391 2600
4 5543 389100307391 2600
5 5544 389100307391 2600
6 5546 381300302513 NA
7 5547 387000307470 NA
8 5548 387000307470 NA
9 5549 387000307470 NA
10 5550 38700030747
I am porting some FORTRAN to R in which an Inf triggers an if(). The
trigger is infinite on exp(). What is the canonical R style of
determining OVER when exp(OVER)== Inf? The code structure that I am
porting is best left intact--so I need to query R somehow to the value
of OFL that causes ex
Thanks to Robert and Gabor for their replies, but neither was what I was
looking for, undoubtedly because of the poor phrasing of my question (+ a
typo -- however, even if correctly typed it doesn't work). I finally
realized that the "elegant" approach I sought can easily be done without
setGener
hi,
i have another question on tapply:
i have a dataset z like this:
5540 389100307391 2600
5541 389100307391 2600
5542 389100307391 2600
5543 389100307391 2600
5544 389100307391 2600
5546 381300302513NA
5547 387000307470NA
5548 387000307470NA
5549 3
R>a<-c(1,2,3)
R> b<-c(1,1,2,4)
R>1*b%in%a*b
[1] 1 1 2 0
R> b<-c(1,1,2,4,5,1,2,3,4)
R>1*b%in%a*b
[1] 1 1 2 0 0 1 2 3 0
I think this works
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of Weiwei Shi
Sent: June 20, 2005 6:35 PM
To: R-help@stat.math.ethz.ch
Su
On Mon, 2005-06-20 at 17:34 -0500, Weiwei Shi wrote:
> Hi,
>
> I have a question, explained by the following example:
> > a<-c(1,2,3)
> > b<-c(1,1,2,4)
> > b[match(a,b, nomatch=0)]
> [1] 1 2
>
> which means it returns "the first match", but I want to get
> 1 1 2 instead of 1 2
>
> In a word, how
On Mon, 2005-06-20 at 16:44 -0400, Hikel, Jerry wrote:
> I am currently attempting to work on some graphs whose plotted lines
> have simple polygons at regular intervals, such as a triangle or a
> square.
>
> I haven't been able to find anything in the base R plotting packages, or
> in any extensi
Hi,
I have a question, explained by the following example:
> a<-c(1,2,3)
> b<-c(1,1,2,4)
> b[match(a,b, nomatch=0)]
[1] 1 2
which means it returns "the first match", but I want to get
1 1 2 instead of 1 2
In a word, how to do multiple matching?
thanks,
--
Weiwei Shi, Ph.D
"Did you always kno
Hi,
I am running R version rw2010 on a Windows 2000 desktop. I am invoking R
from Java via the JGR JRI tools. My process consists of repeated calls to R
in order to create linear models and process the resulting statistics. I
find, however, that the process often dies due to memory allocatio
Here's an example of how to replicate rows according to a count that
is provided in one of the variables.
> foo <- data.frame(id=letters[1:3],cl=LETTERS[1:3],n.yes=c(3,5,2))
> foo
id cl n.yes
1 a A 3
2 b B 5
3 c C 2
> cbind(foo[rep(1:nrow(foo),foo$n.yes),c('id','cl')],res
I am performing a large (105 variable) factor analysis with factanal,
specifying promax rotation. I kow that some loadings over 1.0 are not
unsual with that rotation, but I have some as large as 1.2, which seems
extreme. I am skirting the assumptions of the model by using responses
on a 7-point r
On Thu, Jun 16, 2005 at 10:49:05AM -0500, David Forrest wrote:
> I posted a link to it on the GraphGallery page of the R Wiki:
> http://fawn.unibw-hamburg.de/cgi-bin/Rwiki.pl?GraphGallery
Thanks.
> Why not make it a package?
This is high on the todo list, but first I need to define more option
setMethod("[","foo",function(x,i, j,
.,drop=TRUE)callGeneric([EMAIL PROTECTED],i,j,drop=drop) )
^^^
You have a typo. Use ... instead of . (note ?"[").
Best,
Robert
Robert McGehee
Quantitative Analyst
Geode Capital Management, LLC
53 State Street, 5th Floor | Bo
I am currently attempting to work on some graphs whose plotted lines
have simple polygons at regular intervals, such as a triangle or a
square.
I haven't been able to find anything in the base R plotting packages, or
in any extensions that would allow me to do this easily. I am familiar
with polyg
Federico -
"match" will give you the (first) index of each element of its first
argument in its second argument. So match(vector.1, vector.2) tells you
where each element of vector.1 appears in vector.2. So if you use "table"
on that vector, you'll see how many times each element of vector.2 a
Hello there,
how to do the test for a significant difference between two groups
mean vectors or across several group means vectors?
Thank you,
Shengzhe
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read
v3 <- numeric()
v3[v1] <- table(v2)[v1]
Jim
__
James Holtman"What is the problem you are trying to solve?"
Executive Technical Consultant -- Convergys Labs
[EMAIL PROTECTED]
+1 (513) 723-2929
On 6/20/05, Berton Gunter <[EMAIL PROTECTED]> wrote:
> Folks:
>
> This is a question about the S4 formal class system.
>
> Suppose I have a class, 'foo', defined by:
>
> setClass('foo',representation(dat='matrix', id='character') )
>
> I wish to define a '[' method for foo that will extract fro
On 20 Jun 2005, at 21:24, Erin Hodgess wrote:
> Hello, Federico!
>
> I'm a bit confused about your question, please:
>
> What sorts of things are in Vector1, please?
numbers (as in "numeric") that code individuals
>
> Why are you counting NAs in Vector3, please?
I am counting how many times th
Folks:
This is a question about the S4 formal class system.
Suppose I have a class, 'foo', defined by:
setClass('foo',representation(dat='matrix', id='character') )
I wish to define a '[' method for foo that will extract from the 'dat' slot.
I would have thought that the following would work, b
Hi All,
I am counting the number of occurrences of the terms listed in one
vector in another vector.
My code runs:
for( i in 1:length(vector3)){
vector3[i] = sum(1*is.element(vector2, vector1[i]))
}
where
vector1 = vector containing the terms whose occurrences I want to count
vector2 =
Kevin J Emerson a écrit :
> R-masters,
>
> I have a problem that I have been working on for a while and it seems
> that there may be a simple solution that I have yet to figure out, so I
> thought that I would venture to post to the help list.
>
> Let's say there was a data.frame with three vecto
On Sun, 2005-06-19 at 22:29 +0200, [EMAIL PROTECTED] wrote:
> The Fedora upgrade process should normally install backward compatibility
> libraries when it finds an RPM linked to a library (or library version) that
> isn't in the new release. In this case "compat-libf2c-32" provides libg2c on
> FC
On 6/20/05, Søren Højsgaard <[EMAIL PROTECTED]> wrote:
> Thanks. I wonder if there is a general way of extracting var(u) and var(e),
> which would be needed to simulate u and e. Clearly, one can get the estimated
> parameters, but is there a clever way of 'setting up' the matrices??
> Best
> Sør
R-masters,
I have a problem that I have been working on for a while and it seems
that there may be a simple solution that I have yet to figure out, so I
thought that I would venture to post to the help list.
Let's say there was a data.frame with three vectors, two that are
factors identifying the
I have a question about weights in glm for binomial model. Can I use any
arbitrary number between 0 and 1 as weights instead of integers? I kept getting
warning message when I tried that.
Thanks.
__
R-help@stat.math.ethz.ch mailing list
https://stat.e
I have a question about weights in glm for binomial model. Can I use any
arbitrary number between 0 and 1 as weights instead of integer? I kept getting
warning message when I did that.
Thanks.
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz
Hello there,
Which function can be used to do the test of equal covariance matrix
between two groups or among several groups?
Thank you,
Shengzhe
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the p
'rle' might be your friend. This will find the 'run of a sequence'
Here is some code working off the 'visit' data that you created.
# $Log$
x.1 <- matrix(visit, ncol=4) # your data
x.rle <- apply(x.1, 1, rle) # compute 'rle' for each row
Passed <- lapply(x.rle, function(x){ # now process
On 20 Jun 2005 17:41:13 +0200, Peter Dalgaard <[EMAIL PROTECTED]> wrote:
> Joshua Gilbert <[EMAIL PROTECTED]> writes:
>
> > > > As I understand it, Lapack __can__ computer them
> > > > (http://www.netlib.org/lapack/lawn41/node111.html) and R can use
> > > > Lapack. If there is no function already,
Dear all,
Finally -- after one week of trial & error -- the fog lifted thanks to
the hints of Kazushige Goto.
Therefore, I would like to report back the workaround (for Fedora Core
based systems, including FC4) to the list and a pre-mature analysis,
inspired
by K. Goto:
A brief howto (applicable
On Mon, 2005-06-20 at 08:29 -0400, bogdan romocea wrote:
> Dear useRs,
>
> I timed the same code (simulation with for loops) on the same box
> (dual Xeon EM64T, 1.5 Gb RAM) under 3 OSs and was surprised by the
> results:
> Windows XP Pro (32-bit): Time difference of 5.97 mins
> 64-bit GNU/
Have you considered "simulate.lme"?
Also, have you read Pinheiro and Bates (2000) Mixed-Effects Models in
S and S-Plus (Springer). It is excellent.
spencer graves
Søren Højsgaard wrote:
> Thanks. I wonder if there is a general way of extracting var(u) and var(e),
The problem is that "size" in "dbinom" must be an integer:
> dbinom(x=2, size=4.5, prob=.5)
[1] NaN
Warning message:
NaNs produced in: dbinom(x, size, prob, log)
Consider the following:
> library(MASS)
> set.seed(123)
> x <- rbinom(100, 10, 0.4)
> nx <- max(x)
> fit <- f
R friends,
I am using R 2.1.0 in a Win XP . I have a problem working with lists, probably I
do not understand how to use them.
Lets suppose that a set of patients visit a clinic once a year for 4 years
on each visit a test, say 'eib' is performed with results 0 or 1
The patients do not all visit
The issue here is that the equivalent command array(1:5, c(6,6)) (to
matrix(1:5,6,6)) gives no warning, and sweep uses array().
I am not sure either should: fractional recycling was normally allowed in
S3 (S4 tightened up a bit).
Perhaps someone who thinks sweep() should warn could contribute a
On 6/20/05, Gabor Grothendieck <[EMAIL PROTECTED]> wrote:
> On 6/20/05, Paulo Brando <[EMAIL PROTECTED]> wrote:
> >
> > Dear All,
> >
> > I created a script to calculate averages - two groups: "parcel" and
> > "date" - and, based on these averages, make a graph. The problem is that
> > 'R' does not
On 6/20/05, Paulo Brando <[EMAIL PROTECTED]> wrote:
>
> Dear All,
>
> I created a script to calculate averages - two groups: "parcel" and
> "date" - and, based on these averages, make a graph. The problem is that
> 'R' does not recognize the first column even if I try to insert one.
>
> A brief
Joshua Gilbert <[EMAIL PROTECTED]> writes:
> > > As I understand it, Lapack __can__ computer them
> > > (http://www.netlib.org/lapack/lawn41/node111.html) and R can use
> > > Lapack. If there is no function already, can I access Lapack from R
> > > and use those routines directly?
> >
> > Yes, yo
On Mon, 20 Jun 2005, Joshua Gilbert wrote:
> On 6/17/05, Prof Brian Ripley <[EMAIL PROTECTED]> wrote:
>> On Thu, 16 Jun 2005, Joshua Gilbert wrote:
>>
>>> I need to compute generalized eigenvalues. The eigen function in base
>>> doesn't do it and I can't find a package that does.
>>
>> They are ve
try this:
results.pas <- results[results$Score != "Failed", ]
I hope it helps.
Best,
Dimitris
Dimitris Rizopoulos
Ph.D. Student
Biostatistical Centre
School of Public Health
Catholic University of Leuven
Address: Kapucijnenvoer 35, Leuven, Belgium
Tel: +32/16/336899
Fax: +32/16/337015
Web
You can use 'subset()' if you have a data frame.
d <- read.table(...)
subset(d, Score %in% c("Passed", "Borderline"))
-roger
Bernard L. Dillard wrote:
> Hello. I have looked at R Site Search for this problem, and it didn't
> give me exactly what I needed.
>
> Consider this dataset called "resu
Hello. I have looked at R Site Search for this problem, and it didn't
give me exactly what I needed.
Consider this dataset called "results". It has the following information:
Student DaySubjectScore
Mary 1 Math Failed
David 2 Sci
On 6/20/05, Spencer Graves <[EMAIL PROTECTED]> wrote:
> (comments in line)
>
> Stephen wrote:
> > Dear Fellow R users,
> >
> >
> >
> > I am fairly new to R and am currently conducting a mixed model.
> >
> >
> >
> > I have 7 repeated measures on a simulated clinical trial
> >
> >
> >
> > If I under
Thanks. I wonder if there is a general way of extracting var(u) and var(e),
which would be needed to simulate u and e. Clearly, one can get the estimated
parameters, but is there a clever way of 'setting up' the matrices??
Best
Søren
On 6/19/05, Søren Højsgaard <[EMAIL PROTECTED]> wrote:
> I wou
(comments in line)
Stephen wrote:
> Dear Fellow R users,
>
>
>
> I am fairly new to R and am currently conducting a mixed model.
>
>
>
> I have 7 repeated measures on a simulated clinical trial
>
>
>
> If I understand the model correctly, the outcome is the measure (as a
> factor) the
If I read your question correctly, you want to integrate the
indicated expression over x = 0 to Inf. If I substitute z = x/sqrt(b),
your integral becomes one constant plus another times the expected value
of (y-z/sqrt(b))^2, where z follows Student's t with v degrees of
freedom. Exp
On 6/17/05, Prof Brian Ripley <[EMAIL PROTECTED]> wrote:
> On Thu, 16 Jun 2005, Joshua Gilbert wrote:
>
> > I need to compute generalized eigenvalues. The eigen function in base
> > doesn't do it and I can't find a package that does.
>
> They are very rarely used in statistics, so this is not sur
James
The main reason for the adjusted R^2 (Fisher) is that it is less biased than
the ordinary R^2. The ordinary R^2 has a positive bias that is a function of
the true Rho^2, the number of predictors p, and the sample size n. The maximum
bias occurs at Rho^2 = 0, where the expected R^2 is
Dear All,
I created a script to calculate averages - two groups: "parcel" and
"date" - and, based on these averages, make a graph. The problem is that
'R' does not recognize the first column even if I try to insert one.
A brief example
Raw data:
Data <- sample(1:100, 30, replace = FALSE, pro
Hi Heather
thanks for this. (I have messed up quote level somewhere. Sorry for
that).
On Jun 20, 2005, at 02:01 pm, Heather Turner wrote:
>
Robin Hankin <[EMAIL PROTECTED]> 06/20/05 10:32am >>>
>> Hi
>>
>> [snip]
>> I would expect sweep() to give
>> me a warning if the elements don't rec
?unique
On Mon, 20 Jun 2005 13:50:22 +0100
"Uzuner, Tolga" <[EMAIL PROTECTED]> wrote:
> Hi,
>
> Am sure this is a trivial question but for some reason, haven't been able to
> figure it out.
>
> I want to bucket data in a vector, and then iterate over the buckets.
>
> Say the data set is:
>
Dr H Turner
Research Assistant
Dept. of Statistics
The University of Warwick
Coventry
CV4 7AL
Tel: 024 76575870
>>> Robin Hankin <[EMAIL PROTECTED]> 06/20/05 10:32am >>>
Hi
I had a hard-to-find bug in some of my code the other day, which I
eventually
traced to my misusing of sweep().
I woul
Many thanks all, that was it.
Regards,
Tolga
-Original Message-
From: Marc Schwartz [mailto:[EMAIL PROTECTED]
Sent: 20 June 2005 14:02
To: Uzuner, Tolga
Cc: 'r-help@stat.math.ethz.ch'
Subject: Re: [R] Bucketting data
On Mon, 2005-06-20 at 13:50 +0100, Uzuner, Tolga wrote:
> Hi,
>
> Am
On Mon, 2005-06-20 at 13:50 +0100, Uzuner, Tolga wrote:
> Hi,
>
> Am sure this is a trivial question but for some reason, haven't been
> able to figure it out.
>
> I want to bucket data in a vector, and then iterate over the buckets.
>
> Say the data set is:
>
> > cleandata[,4]
> [1] 26 26
try this:
unique(cleandata[, 4])
length(unique(cleandata[, 4]))
I hope it helps.
Best,
Dimitris
Dimitris Rizopoulos
Ph.D. Student
Biostatistical Centre
School of Public Health
Catholic University of Leuven
Address: Kapucijnenvoer 35, Leuven, Belgium
Tel: +32/16/336899
Fax: +32/16/337015
?table maybe is what you want
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of Uzuner, Tolga
Sent: June 20, 2005 8:50 AM
To: 'r-help@stat.math.ethz.ch'
Subject: [R] Bucketting data
Hi,
Am sure this is a trivial question but for some reason, haven't been a
I'm not certain what "buckets" are, but based on your example try
?unique
On 6/20/05, Uzuner, Tolga <[EMAIL PROTECTED]> wrote:
> Hi,
>
> Am sure this is a trivial question but for some reason, haven't been able to
> figure it out.
>
> I want to bucket data in a vector, and then iterate over the
It is a nevv category called smart luxury goods. They are from Rolexes,
Cartiers, Bvlgaries, FrankMullers, Harry Winstons, Breguets,
Jaeger-lecoultre, Brietilings, TagHeuers and Tudors. With one closer l00k,
you will be convinced at these beauties. They also have their SerialNumbers
and logo. We pa
Hi,
Am sure this is a trivial question but for some reason, haven't been able to
figure it out.
I want to bucket data in a vector, and then iterate over the buckets.
Say the data set is:
> cleandata[,4]
[1] 26 26 26 26 26 26 26 26 26 26 26 26 61 61 61 61 61 61
61 61 61
Dear useRs,
I timed the same code (simulation with for loops) on the same box
(dual Xeon EM64T, 1.5 Gb RAM) under 3 OSs and was surprised by the
results:
Windows XP Pro (32-bit): Time difference of 5.97 mins
64-bit GNU/Linux (Fedora Core 4): Time difference of 6.97 mins
32-bit GNU/Li
I've rolled up R-2.1.1.tar.gz a short while ago. This is a maintenance
release containing mainly bugfixes.
See the full list of changes below.
You can get it from
http://cran.r-project.org/src/base/R-2/R-2.1.1.tar.gz
(give it some time to arrive there) or wait for it to be mirrored at a
CRAN s
There is now an Epi-package on CRAN.
It is intended for epidemiological analysis in R.
It has its own homepage, http://www.pubhealth.ku.dk/~bxc/Epi
The package has been used at the course "Statistical practise
in Epidemiology with R", see http://www.pubhealth.ku.dk/~bxc/SPE.
A mailing list for
Hi,
how can i run SJava from an IDE like Netbeans? If i try to run a example
class from RJava i get this message:
Exception in thread "main" java.lang.UnsatisfiedLinkError: no
RInterpreter in java.library.path
at java.lang.ClassLoader.loadLibrary(ClassLoader.java:1682)
at java.la
hi all
at the outset i must APOLOGIZE for sending the following mail. it is not
R related but since there are many stats and maths buffs that use the
list i decided to send the following question.
integrate ((1+((y-bx)^2)/(av))*(1+(x^2)/(bv)))^(-0.5*(v+1))
over the interval 0 to inf
a>0, b>0 a
On Mon, 20 Jun 2005, Guillaume Allain wrote:
Either download the source code and unpack the file you need, or use the
online anonymous SVN facility:
https://svn.r-project.org/R/trunk/src/library/stats/src/distance.c
> Hi,
>
> Is it possible to access to the C source code of the "dist" functio
you can access "distance.c" in R-sources in the directory
\src\library\stats\src
I hope it helps.
Best,
Dimitris
Dimitris Rizopoulos
Ph.D. Student
Biostatistical Centre
School of Public Health
Catholic University of Leuven
Address: Kapucijnenvoer 35, Leuven, Belgium
Tel: +32/16/336899
Fax
Hi
I had a hard-to-find bug in some of my code the other day, which I
eventually
traced to my misusing of sweep().
I would expect sweep() to give
me a warning if the elements don't recycle nicely, but
X <- matrix(1:36,6,6)
sweep(X,1,1:5,"+")
[,1] [,2] [,3] [,4] [,5] [,6]
[1,]2
Hi,
Is it possible to access to the C source code of the "dist" function? I
mean when u do read the dist code, it points almost directly to a C
function :
d <- .C("R_distance", x = as.double(x), nr = N, nc = ncol(x),
d = double(N * (N - 1)/2), diag = as.integer(FALSE),
method
Dear Fellow R users,
I am fairly new to R and am currently conducting a mixed model.
I have 7 repeated measures on a simulated clinical trial
If I understand the model correctly, the outcome is the measure (as a
factor) the predictors are clinical group and trial (1-7). The fixed
facto
79 matches
Mail list logo