Yes!!! I have overcome it!
The solution is
xyplot((V_A3/25)~year|factor(month),
strip=strip.custom(factor.levels=c(month.name,"Annual mean",
"Annual integral","Integral (Jan-Mar)","Integral
(Apr-Jun)","Integral (Jul-Sep)","Integral (Oct-Dec)",
and s
Dear r-help,
I need to draw some lattice plots and have stuck in strange problems.
I have a data frame of 405 rows:
> str(dframe)
`data.frame': 405 obs. of 4 variables:
$ year : num 1990 1990 1990 1990 1990 1990 1990 1990 1990 1990 ...
$ month: num 1 2 3 4 5 6 7 8 9 10 ...
$ V_A2 : n
How about the following:
> tstFn <- function(x){
+ xName <- deparse(substitute(x))
+ cat(xName)
+ "done"
+ }
> tstFn(Varname)
Varname[1] "done"
>
> tstF1 <- function(x){
+ tstF2 <- function(y){
+ cat(deparse(substitute(y)),
+ "doesn't work\n")
+ }
+ tstF2(x)
+
It's not so simple, but consider the following:
x=rep(1:6, each=2, length=24)
y=rep(1:6, each=4)
z=rep(0:1, length=24)
set.seed(1)
tstDF <- data.frame(x=x, y=y, z=z,
w=(x-3.5)^2+(y-3.5)+z+0.1*rnorm(24))
fit <- lm(w~x+y+z+I(x^2)+I(y^2), tstDF)
x0 <- seq(1, 5, .5)
y0 <- seq(1,
On Jul 16, 2005, at 9:58 AM, Thomas Lumley wrote:
> On Fri, 15 Jul 2005, David Groos wrote:
>
>> Thank you-all very much for your help, your responses and help has
>> been
>> very encouraging. The following doesn't close the case but it tables
>> it...
>>
>> First I copied Ken's code into my R
R2.1.1
Win 2k
I have a function, B, within a function, A.
I would like to have B print the name of the argument passed to it (not
the value of the arguments).
i.e.,
A<-function()
{
B<-function(x,y)
{
fit1<-lm(y~x,data=jo)
print(summary(fit1)
I want B to prin
Have you tried "anova(fit1, fit2)", where
fit1 <- lme(one model...)
fit2 <- lme(a submodel ... )
This "anova" does about the best that anyone knows how to do -- or at
lest did 7 years ago when it was written. If the "submodel" changes the
fixed effects, you should
Have you considered "lme" in library(nlme)? If you want to go this
route, I recommend Pinheiro and Bates (2000) Mixed-Effect Models in S
and S-Plus (Springer).
spencer graves
Brett Gordon wrote:
> Hello,
>
> I'm trying to model the entry of certain firms into a larger num
?try
spencer graves
Yimeng Lu wrote:
> Hello R-users,
>
> In a loop, if a function, such as "nls", gives an error, is it possible to
> coerce R to continue proceeding the next command with the same
> loop?
>
> Thanks so much for your advice!
>
> Hanna Lu
> [[alternati
Yimeng Lu wrote:
> Hello R-users,
>
> In a loop, if a function, such as "nls", gives an error, is it possible to
> coerce R to continue proceeding the next command with the same
> loop?
>
Yes, see the try() function. The basic usage is something like
value <- try( some calculation )
if (inheri
Hello R-users,
In a loop, if a function, such as "nls", gives an error, is it possible to
coerce R to continue proceeding the next command with the same
loop?
Thanks so much for your advice!
Hanna Lu
[[alternative HTML version deleted]]
__
R-h
Hello,
I'm trying to model the entry of certain firms into a larger number of
distinct markets over time. I have a short time series, but a large
cross section (small T, big N).
I have both time varying and non-time varying variables. Additionally,
since I'm modeling entry of firms, it seems like
Hello,
I'm trying to model the entry of certain firms into a larger number of
distinct markets over time. I have a short time series, but a large
cross section (small T, big N).
I have both time varying and non-time varying variables. Additionally,
since I'm modeling entry of firms, it seems like
Hi,
I hope this is the right list for my posting, since I've never posted
to any R list before.
I'm quite extensively using the xfig graphics device and as far as I
figured out this
device writes all the objects into xfig layer 100 (based on what I saw
in the devPS.c
file -if this is the file t
[Thomas Lumley]
> It would be nice if R could realize that you meant the cube root
> of -8, but that requires either magical powers or complicated and
> unreliable heuristics. The real solution might be a function like
> root(x,a,b) to compute x^(a/b), where a and b could then be exactly
> repres
On 7/16/05, Jeff Newmiller <[EMAIL PROTECTED]> wrote:
> I have a rather basic background in statistics, and am looking for
> assistance in solving what I expect is a common type of problem.
>
> I have measurements of physical processes, and mathematical models of
> those processes that I want to f
Hi,
Is there a way, preferably with R, to read shapefiles and transform
them in a format that I could then use with package PBSmapping?
I have been able to read such files into R with maptools' read.shape
and plot it with plot.Map, but I'd like to bring the data to
PBSmapping and plot from
On 7/16/2005 1:17 PM, Vivek Rao wrote:
> I would like to see R packages arranged by topic. CRAN
> Task Views are at
> http://lib.stat.cmu.edu/R/CRAN/src/contrib/Views/ ,
> but I'd like something more detailed. For example, the
> IMSL Fortran library, version 4 is easy to navigate
> and has procedu
I have a rather basic background in statistics, and am looking for
assistance in solving what I expect is a common type of problem.
I have measurements of physical processes, and mathematical models of
those processes that I want to feed the measurements into. A simple case
is using measurements o
When I draw a ROC curve, the xlabel is '1-spec:', and ylabel is 'sens:' . My
question is: can I specify them to '1-specificity', and 'sensitivity' ? thanks
for your help. Lu.
__
[[alternative HTML version deleted]]
__
I would like to see R packages arranged by topic. CRAN
Task Views are at
http://lib.stat.cmu.edu/R/CRAN/src/contrib/Views/ ,
but I'd like something more detailed. For example, the
IMSL Fortran library, version 4 is easy to navigate
and has procedures arranged according to following
topics:
Bas
I suspect that there are more direct ways to do this test, but it
is unclear to me just what the issue is. For example, if there
are many subjects and very few stimuli for each, you might want
to get some sort of measure of ability for each subject (many
possibilities here, then test the measure
On Sat, 16 Jul 2005, Thomas Lumley wrote:
>
> Yes, but you don't need to go via the baseline. The survival curves for
> any two covariate vectors z1 and z2 are related by
>
> S(t; z1)= S(t; z2)^(z1-z2)
>
Actually
S(t; z1)=S(t;z2) ^(beta'(z1-z2))
of course.
-thomas
___
On Sat, 16 Jul 2005, Bernardo Rangel Tura wrote:
> At 10:11 12/7/2005, [EMAIL PROTECTED] wrote:
>
>> hi all
>>
>> why does R do this:
>>
>> (-8)^(1/3)=NaN
>>
>> the answer should be : -2
>
Yes and no.
The problem is that the reciprocal of 3 is not exactly representable as a
floating point numbe
On Sat, 16 Jul 2005, Kylie-Anne Richards wrote:
> Thank you for your help.
>
>> In any case, to specify f.pom You need it to be a factor with the same set
>> of levels. You don't say what the lowest level of pom is, but if it is,
>> s
On Fri, 15 Jul 2005, David Groos wrote:
> Thank you-all very much for your help, your responses and help has been
> very encouraging. The following doesn't close the case but it tables
> it...
>
> First I copied Ken's code into my R Console and...it worked great!
> That was baffling as it looked
Thanks to all kind people who answered. Attached is a useful reply,
among many others.
--
Mauro Gasparini
Professore Straordinario di Statistica
Dipartimento di Matematica, Politecnico di Torino
Corso Duca degli Abruzzi 24 I-10129 Torino, Italy
tel: +39 011 564 7546
fax: +39 011 564 7599
ema
On Friday 15 July 2005 17:00, Deepayan Sarkar wrote:
> On 7/15/05, Federico Gherardini <[EMAIL PROTECTED]> wrote:
> > On Friday 15 July 2005 14:42, you wrote:
> > > Hi all,
> > > I've used the split argument to print four lattice plots on a single
> > > page. The problem now is that I need to reduc
> Didn't you simply try:
>
> > A<-matrix(c(1.1,1.2,1.3,1.4,1.5,1.6),ncol=3)
> > B<-matrix(c(2.1,2.2,2.3,2.4,2.5,2.6),ncol=3)
> > C<-matrix(c(3.1,3.2,3.3,3.4,3.5,3.6),ncol=3)
> > A
> [,1] [,2] [,3]
> [1,] 1.1 1.3 1.5
> [2,] 1.2 1.4 1.6
> > B
> [,1] [,2] [,3]
> [1,] 2.1 2.3 2.5
> [
Have you considered "BTm" in library(BradleyTerry)? Consider the
following example:
> cond1 <- data.frame(winner=rep(LETTERS[1:3], e=2),
+ loser=c("B","C","A","C","A","B"),
+ Freq=1:6)
> cond2 <- data.frame(winner=rep(LETTERS[1:3], e=2),
+ loser=c("B","C
Uzuner, Tolga wrote:
> Hi Folks,
> Has anyone been able to get rmpi to work under windows ?
Rmpi uses the LAM implementation of MPI.
See
http://www.lam-mpi.org/faq/category12.php3
and read FAQ 2 which implicitly tells us that there is no native port,
hence you cannot run it under Windows.
The pac
When you copied it into MS Word, did you "Copy -> "paste special" ->
"unformatted text"? This trick sometimes might expose (or eliminate)
nonprinting characters or characters with special attributes that might
have created the problem.
spencer graves
David Groos wrote:
>
> -Original Message-
> From: [EMAIL PROTECTED]
> [mailto:[EMAIL PROTECTED] On Behalf Of
> [EMAIL PROTECTED]
> Sent: 16 July 2005 02:03
> To: r-help@stat.math.ethz.ch
> Subject: [R] question from environmental statistics
>
> thanks Fran. that was useful but Im still in a fix. its a
> rea
(Ted Harding) <[EMAIL PROTECTED]> writes:
> On 16-Jul-05 Mauro Gasparini wrote:
> >
> > Dear users,
> >
> > I have a list of several matrices with the same number of columns,
> > how do I rbind them all with a vectorized command?
> >
> > A related simpler question is, how do I vectorize the ins
On 16-Jul-05 Mauro Gasparini wrote:
>
> Dear users,
>
> I have a list of several matrices with the same number of columns,
> how do I rbind them all with a vectorized command?
>
> A related simpler question is, how do I vectorize the instruction
> that rbinds together several copies of the same
At 10:11 12/7/2005, [EMAIL PROTECTED] wrote:
>hi all
>
>why does R do this:
>
>(-8)^(1/3)=NaN
>
>the answer should be : -2
Allan
In my computer:
> (-8)^(1/3)
[1] NaN
> -8^(1/3)
[1] -2
> -(8^(1/3))
[1] -2
The problem is -8 or the problem is (-8) ?
[]s
Tura
--
No virus found in this o
Dear users,
I have a list of several matrices with the same number of columns,
how do I rbind them all with a vectorized command?
A related simpler question is, how do I vectorize the instruction
that rbinds together several copies of the same matrix?
Thanks.
--
Mauro Gasparini
_
Thank you for your help.
> In any case, to specify f.pom You need it to be a factor with the same set
> of levels. You don't say what the lowest level of pom is, but if it is,
> say, -3.
>
> f.pom=factor(-3, levels=seq(-3,2.5, by=0.5))
-lf77blas is part of ATLAS, so I do suspect the RPM builder had ATLAS
installed.
lme4 needs a compatible Matrix installed.
I do think installing from the sources would solve this, but probably you
need to discuss this with the RPM maintainer as a dependency appears to
be missing.
On Fri, 15 J
39 matches
Mail list logo