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
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 statistics.
Most likely your linear model fitting is too large for the amount of
memory you have in your Windows box. However, you should check ?Memory
and in particular how your maximum RAM usage is set. (I don't know how
the `JGR JRI tools' give you access to R's startup settings.)
It may be that you
Hi again,
thank you for your previous answers. Just another question, though ...
I get the following variance components after fitting a mixed model.
Groups Name Variance Std.Dev. Corr
PlantID TreatmCtrl0.51784 0.71961
TreatmNoAccess4.77469 2.18511
AFAIK the COVTEST option just computes a Wald test! Since you want
test for variance components (which is on the boundary of the
parameter space), I'd suggest to use a LRT (i.e., anova.lme(model.1,
model.2)) and moreover consider the simulate.lme() function of the
nlme package.
Best,
Hi R users,
I use plot.circular(rad, stack=T,bins=4) and could just obtain the four
stacks in the directions of 45, 135, 225 and 315 degrees. But I want them in
0, 90, 180 and 270 degrees. Is there any parameter or any other function to
make it? Thanks for helping a R beginner.
--
Xiaohua
Hello,
i have a problem with the function nls().
This are my data in k:
V1V2
[1,]0 0.367
[2,] 85 0.296
[3,] 122 0.260
[4,] 192 0.244
[5,] 275 0.175
[6,] 421 0.140
[7,] 603 0.093
[8,] 831 0.068
[9,] 1140 0.043
With the nls()-function i want to fit following
OTOH, the Opteron Goto BLAS has no such problems on FC3, so I suspect this
is specific to only some of the architectures.
I've added a comment in the R-admin manual that says
It has been reported that on some RedHat-based Linux systems it is
necessary to set @code{OMP_NUM_THREADS=1} in
G'day Chris,
CK == Christfried Kunath [EMAIL PROTECTED] writes:
CK With the nls()-function i want to fit following formula
CK whereas a,b, and c are variables: y~1/(a*x^2+b*x+c)
CK [...]
CK The algorithm plinear give me following error:
The algorithm plinear is inappropriate
The Department of Statistics and Mathematics at the Vienna University of
Economics and Business Administration invites applications for two new
faculty positions in computational statistics and quantitative research
methodology, to begin in fall 2005. The positions will be at the
Assistant
Hey,
i'm trying to set up an cgi web-application that produces maps and other
graphics with R. for local testing i used my laptop which has an win xp OS.
there everything works fine. now i try to move the whole system over on an unix
machine, mainly because i want to establish a connection to
On 6/21/05, Christfried Kunath [EMAIL PROTECTED] wrote:
Hello,
i have a problem with the function nls().
This are my data in k:
V1V2
[1,]0 0.367
[2,] 85 0.296
[3,] 122 0.260
[4,] 192 0.244
[5,] 275 0.175
[6,] 421 0.140
[7,] 603 0.093
[8,] 831 0.068
Hello there,
I have two data sets with 14 variables each, and wish to do the test
for equality of their covariance matrices and mean vectors. Normally
these tests should be done by chi square test (box provided) and
Hotelling's T square test respectively. Which R functions could do
this kind of
Hello,
I want to create a data.frame with different number of columns per line. What I
want is something like:
example - NULL
begin - 1
while (end nrow(orig.data))
{
end - next.day(orig.data,begin) # my own function. It returns the first index
from the next day. Each day has a different
On 6/21/05, Gabor Grothendieck [EMAIL PROTECTED] wrote:
On 6/21/05, Christfried Kunath [EMAIL PROTECTED] wrote:
Hello,
i have a problem with the function nls().
This are my data in k:
V1V2
[1,]0 0.367
[2,] 85 0.296
[3,] 122 0.260
[4,] 192 0.244
[5,]
On Jun 20, 2005, at 04:58 pm, Prof Brian Ripley wrote:
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).
I have found this a useful correspondence. My own wish is to define a new
S4 class which differs from
class 'array' only in its default handling of the 'drop' argument i.e.
drop=FALSE - not an unusual wish.
My solution is the following:
I don't think conversion to a list is necessary, e.g.
f - function(x, y, ...) {
cl - sys.call()
if (missing(y)) cl$y - FALSE
cl[[1]] - as.name(cat)
cl[[2]] - rev(x)
eval(cl)
}
f(1:4) # 4 3 2 1 FALSE
On 6/21/05, [EMAIL PROTECTED] [EMAIL PROTECTED] wrote:
\
Perhaps the signature should be:
sweep(...other args go here..., warn=getOption(warn))
so that the name and value of the argument are consistent with
the R warn option.
On 6/21/05, Robin Hankin [EMAIL PROTECTED] wrote:
On Jun 20, 2005, at 04:58 pm, Prof Brian Ripley wrote:
The issue
Hi
searching in CRAN homepage for hotelling gave me this response
Thanks everyone for your help on this question. I solved the
problem by
writing a procedure to calculate Hotelling's T^2 for a one-sample,
multivariate t-test. Here's how it looks, perhaps it will be useful to
others.
data
Hi
from help page
Details:
A data frame is a list of variables of the same length with unique
^^
row names, given class 'data.frame'.
So you need to use a different sructure (e.g.list) or you have to
ensure that you rbind
hi netters
Suppose I have a factor X, with 10 elements and 3 levels: A B B C A C B A C
C .
It is easy to count the number of elements for each level:
tapply(X,X,length).
Now I have another factor Y, which formed a matrix with X:
X| A B B C A C B A C C
Y| B B C C C A A A B B
I wanna count
table() is your friend:
X - factor(scan(what=))
1: A B B C A C B A C C
11:
Read 10 items
Y - factor(scan(what=))
1: B B C C C A A A B B
11:
Read 10 items
table(X)
X
A B C
3 3 4
table(X, Y)
Y
X A B C
A 1 1 1
B 1 1 1
C 1 2 1
Andy
From: zhihua li
hi netters
Suppose I
Hi,
I just find another package shapes for solving the Hotelling T^2
test, but what about the R function for the test of covariance
equality(chi ^ 2 box's test) between two multidimensional data sets?
Shengzhe
__
R-help@stat.math.ethz.ch mailing list
try this:
X - factor(sample(c(A, B, C), 10, TRUE), levels = c(A, B,
C))
Y - factor(sample(c(A, B, C), 10, TRUE), levels = c(A, B,
C))
#
table(X, Y)
I hope it helps.
Best,
Dimitris
Dimitris Rizopoulos
Ph.D. Student
Biostatistical Centre
School of Public Health
Catholic
wu sz [EMAIL PROTECTED] writes:
Hi,
I just find another package shapes for solving the Hotelling T^2
test, but what about the R function for the test of covariance
equality(chi ^ 2 box's test) between two multidimensional data sets?
Shengzhe
Rather than looking for a package, it might
On 6/21/05, Gabor Grothendieck [EMAIL PROTECTED] wrote:
On 6/21/05, Gabor Grothendieck [EMAIL PROTECTED] wrote:
On 6/21/05, Christfried Kunath [EMAIL PROTECTED] wrote:
Hello,
i have a problem with the function nls().
This are my data in k:
V1V2
[1,]0 0.367
I think the warning condition in Robin's patch is too harsh - the following
examples seem reasonable to me, but all produce warnings
sweep(array(1:24, dim = c(4,3,2)), 1, 1:2, give.warning = TRUE)
sweep(array(1:24, dim = c(4,3,2)), 1, 1:12, give.warning = TRUE)
sweep(array(1:24, dim = c(4,3,2)),
I need to do some simulations based on a nls model. I have been able to
use nls to determine the model and I have almost gotten to the point at
which I can construct the function. When I get the function done I will
use uniroot to solve the function at a specific y value for x. I found
an example
Hi
On Jun 21, 2005, at 02:33 pm, Heather Turner wrote:
I think the warning condition in Robin's patch is too harsh - the
following examples seem reasonable to me, but all produce warnings
sweep(array(1:24, dim = c(4,3,2)), 1, 1:2, give.warning = TRUE)
sweep(array(1:24, dim = c(4,3,2)), 1,
I am trying to build a package which calls some c code. I want to call the
dyn.load function when the package is loaded. To do this I create a
.First.lib function in MyConv.R file:
.First.lib-function(/home/phineas/R_HOME/R-2.0.1/library/, MyConvolve){
### Intentionally empty body. To
Dear Petr and Shengzhe,
Also see ?manova.
John
John Fox
Department of Sociology
McMaster University
Hamilton, Ontario
Canada L8S 4M4
905-525-9140x23604
http://socserv.mcmaster.ca/jfox
-Original Message-
From: [EMAIL
On Tue, 2005-06-21 at 06:57 -0400, Gabor Grothendieck wrote:
On 6/21/05, Christfried Kunath [EMAIL PROTECTED] wrote:
Hello,
i have a problem with the function nls().
This are my data in k:
V1V2
[1,]0 0.367
[2,] 85 0.296
[3,] 122 0.260
[4,] 192 0.244
Phineas Campbell wrote:
I am trying to build a package which calls some c code. I want to call the
dyn.load function when the package is loaded. To do this I create a
.First.lib function in MyConv.R file:
.First.lib-function(/home/phineas/R_HOME/R-2.0.1/library/, MyConvolve){
Your
As far as I know, your options are (1) run Xvfb on the unix box, or
(2) use the bitmap() device instead of jpeg(). Option 2 will require
that ghostscript also be installed on the unix box.
Strictly speaking, it doesn't have to to be Xvfb. The requirement is
that an X server be running.
The R
On 6/21/05, Manuel Morales [EMAIL PROTECTED] wrote:
On Tue, 2005-06-21 at 06:57 -0400, Gabor Grothendieck wrote:
On 6/21/05, Christfried Kunath [EMAIL PROTECTED] wrote:
Hello,
i have a problem with the function nls().
This are my data in k:
V1V2
[1,]0 0.367
hi
i tried all the methods suggested above:
ave and rowsum with with function works for my situation. I think
the problem might not be due to tapply.
My data z comes from
z-y[y[[1]] %in% x[[2]], c(1,9)]
while z is supposed to have no entries for those non-matched between x and y.
however, when I
Agreed my examples may be trivial and I'm sure there are more efficient ways to
do the same thing, but I disagree that my examples are against the spirit of
sweep (). In the first example a vector is swept out from the rows, with one
value for odd rows and one value for even rows. In the second
The outcome variable in my dataset is cost. I prefer not to transform it as
readers will really be interested in pounds, not log(pounds) or
sqrt(pounds). I have fitted my model using lme and now wish to use boot on
it. I have therefore plagiarised the example in the article in RNews 2/3
What does str(z) say? I suspect the second column is a factor, which, after
the subsetting, has some empty levels. If so, just drop those levels.
Andy
From: Weiwei Shi
hi
i tried all the methods suggested above:
ave and rowsum with with function works for my situation. I think
the
We don't have the structure of your dataset. But it seems pretty clear
that your resampling is not preserving the random effects structure: you
are always fitting to the same clusters labelled by gp and so missing the
major source of variability.
You either need to resample clusters, or
To reinforce, Dimitris' excellent suggestion, I'd like to refer you
to sect. 2.4 in Pinheiro and Bates (2000) Mixed-Effects Models in S and
S-Plus (Springer, pp. 92-96). Testing for parameters on a boundary of a
parameter space, as, e.g., whether a variance component is zero,
Even before I tried, I already realize it must be true when I read
this reply! Great job! thanks, Andy.
str(z)
`data.frame': 235 obs. of 2 variables:
$ CLAIMNUM : Factor w/ 1907 levels 0,1001849,..: 1083 1083
1083 1582 1582 1084 1681 1681 1391 1391 ...
$ SIU.SAVED: int 475 3000 3000 0
Try:
(x - factor(1:2, levels=1:5))
[1] 1 2
Levels: 1 2 3 4 5
(x - x[, drop=TRUE])
[1] 1 2
Levels: 1 2
Andy
From: Weiwei Shi [mailto:[EMAIL PROTECTED]
Even before I tried, I already realize it must be true when I read
this reply! Great job! thanks, Andy.
str(z)
`data.frame': 235
XSolutions Corp (www.xlsolutions-corp.com) is proud to announce
our Advanced R/Splus programming course taught by R Development
Core Team Guru!
www.xlsolutions-corp.com/Radv.htm
*New York --- July 28th-29th, 2005
*San Francisco -- July 18th-19th, 2005
Hi R users,
I use plot.circular(rad, stack=T,bins=4) and could just obtain the four
stacks in the directions of 45, 135, 225 and 315 degrees. But I want them in
0, 90, 180 and 270 degrees. Is there any parameter or any other function to
make it? Thanks for helping a R beginner.
--
Xiaohua
The problem with simulate.lme is that it only returns logL for a given model
fitted to a simulated data set - not the simulated data set itself (which one
might have expected a function with that name to do...). It would be nice with
that functionality...
Søren
On 6/21/05, Søren Højsgaard [EMAIL PROTECTED] wrote:
The problem with simulate.lme is that it only returns logL for a given model
fitted to a simulated data set - not the simulated data set itself (which
one might have expected a function with that name to do...). It would be nice
with
On Tue, 21 Jun 2005, Douglas Bates wrote:
On 6/21/05, Søren Højsgaard [EMAIL PROTECTED] wrote:
The problem with simulate.lme is that it only returns logL for a given model
fitted to a simulated data set - not the simulated data set itself (which one
might have expected a function with that
Hi
Uwe Ligges wrote:
Jim Lemon wrote:
Thanks for the answers - I should have been more specific as I had
already tried axTicks and pretty.
The function in question is gantt.chart() in the latest plotrix
package (Thanks to Scott Waichler for the original code). I settled on
Background:
OS: Linux Mandrake 10.1
release: R 2.0.0
editor: GNU Emacs 21.3.2
front-end: ESS 5.2.3
-
Colleagues
Is there a function in R that is an equivalent of zoom in matlab? This is
very useful for being able to magnify details in a plot.
I have searched the
On 6/21/05, McClatchie, Sam (PIRSA-SARDI)
[EMAIL PROTECTED] wrote:
Background:
OS: Linux Mandrake 10.1
release: R 2.0.0
editor: GNU Emacs 21.3.2
front-end: ESS 5.2.3
-
Colleagues
Is there a function in R that is an equivalent of zoom in matlab? This is
Try RSiteSearch(zoom plot). There are some good suggestions there.
Cheers
Francisco
From: McClatchie, Sam (PIRSA-SARDI) [EMAIL PROTECTED]
To: R-Help-Request (E-mail) r-help@stat.math.ethz.ch
Subject: [R] R-help
Date: Wed, 22 Jun 2005 09:05:05 +0930
Background:
OS: Linux Mandrake 10.1
Such a wide range of backgrounds here, I thought I'd toss this out here to
get ideas.
I've lucked into some clinical trial data where schizophrenic patients were
randomly assigned to start on one of three drugs, then were followed
naturalistically over a year (or more, depending on when they
Dear All,
I am not familar with R. I want to use PCA (principal components
analysis) and MDS (multidimensional scaling). Can someone tell me
which R package I should use for PCA and MDS? I appreciate your help
in advance.
Ray
__
prcomp(stats) Principal Components Analysis
princomp(stats) Principal Components Analysis
scale(base) Scaling and Centering of Matrix-like Objects
2005/6/22, Ray Xiong [EMAIL PROTECTED]:
Dear All,
I am not familar with R. I want to use PCA (principal components
analysis) and MDS
56 matches
Mail list logo