Dear helper.
I am a beginer.
I have difficulties to handle axis. I want to draw axis label such that axis has range
of [-0.4,0.4] with intervel 0.2 for x and y axis.
Some part of range do not have data points. Thus, plot does not show whole range. How
can I enforce plot to depict the whole rang
On Tue, Oct 21, 2003 at 08:31:16PM -0400, Gabor Grothendieck wrote:
>
>
> I have Excel files containing data that I would like to move to R.
> They are in the standard form of a one row header followed by
> rows of data, one record per row EXCEPT that there are a few
> rows of comments before th
> From: Jason Turner [mailto:[EMAIL PROTECTED]
>
> Zhen Pang wrote:
>
> > We are not allowed to submit job directly, so I never type
> R to use R,
> > just make a batch. How can I use try() to correct my codes? In the
> > interactive mode, I know how to continue, but now I never
> enter the R
I have Excel files containing data that I would like to move to R.
They are in the standard form of a one row header followed by
rows of data, one record per row EXCEPT that there are a few
rows of comments before the header. The number of rows of comments
varies. For Excel files of this form
Prof. Bates may be able to give us more recent references on this, but
the best literature I know on this is Pinhiero and Bates (2000)
Mixed-Effects Models in S and S-Plus (Springer, sec. 2.4). This
includes description of a "simulate.lme" function, which you can use to
generate random numbers
Thanks to Thomas, Roger and Spencer on your reply!
I am trying to replicate/validate the results from the "black box"
Econometrics software Eviews, compared to R and EViews uses the Almon
approach. I will try the code kindly offered by Thomas.
I would be interested in the way to apply the B-spl
Hi
Arend P. van der Veen wrote:
Hi,
Can anyone point in the direction of how to write a custom output device
in R. I currently generate output using PS and JPEG but need to produce
output in our own vector graphics language.
Look at ...
R/src/include/R_ext/GraphicsDevice.h
... for the in-pro
Hi all,
I'm trying to access my account at
https://www.stat.math.ethz.ch/mailman/options/r-help
but the following appears:
*Error: */Authentication failed./
Yes, I've already checked my password. Is it a problem
with mailman or something?
Thanks in advance.
Hector
__
Hi,
Can anyone point in the direction of how to write a custom output device
in R. I currently generate output using PS and JPEG but need to produce
output in our own vector graphics language.
Thanks in advance,
Arend van der Veen
__
[EMAIL PROTECTED]
David Brahm <[EMAIL PROTECTED]> writes:
> Language differences:
> - Scoping rules differ. In R, functions see the functions they're in. Try:
> f1 <- function() {x <- 1; f2 <- function() print(x); f2()}; f1()
It may be more accurate to say "the functions they're defined in".
It's lexical s
Contributions of code to provide alternative calculations of
denominator degrees of freedom are welcome :-)
I think it would be good to bear in mind that the use of the t and F
distributions for models with mixed effects is already an
approximation. If your design is such that you end up with a v
According David Jame's response to my earlier question, there is a
problem with setGeneric.setMethod in R 1.8.0 that affects DBI and RMySQL.
Is there a fix for this? David Jame's refers to an 'R-patched version'
but I haven't seen anything like this on CRAN. Would going back to an
older ver
Zhen Pang wrote:
We are not allowed to submit job directly, so I never type R to use R,
just make a batch. How can I use try() to correct my codes? In the
interactive mode, I know how to continue, but now I never enter the R
window, where to find my results and save seed to continue?
Like you'
It's actually easier than I thought (and perhaps than you described). What
Prof. Ripley suggested, combined with the options to the configure script,
should make this fairly straight-forward, even with "make install".
Cheers,
Andy
> From: David Brahm [mailto:[EMAIL PROTECTED]
>
> Andy Liaw <[E
On Tue, 21 Oct 2003, Robin Hankin wrote:
> Hi everyone
>
> I've been playing with do.call() but I'm having problems understanding it.
>
> I have a list of "n" elements, each one of which is "d" dimensional
> [actually an n-by-n-by ... by-n array]. Neither n nor d is known in
> advance. I want to
Andy Liaw <[EMAIL PROTECTED]> wrote:
> What's not clear to me is a good way of keeping two versions of R
> simultaneously (for ease of transition). Can anyone suggest a good strategy
> for doing that on *nix?
I'm not sure what you mean, but I'll tell you what we do. We have built
/res/R/R-1.6.2
Purvis Bedenbaugh <[EMAIL PROTECTED]> wrote:
> I started looking for an R-S compatibility table but didn't find it.
> Examples:
> 'stdev' is now 'sd' - is it exactly the same computation ?
> couldn't find a built-in for error.bar()
> syntax that is an error in R: param(thisframe,"b") <- valu
Hello all.
I was wondering if there is any way to adjust the denominator degrees of
freedom in lme(). It seems to me that there is only one method that can be
used. As has been pointed out previously on the list, the denominator
degrees of freedom given by lme() do not match those given by SAS
Hi
Vittorio wrote:
Paul Murrell [r-help] <20/10/03 09:13 +1300>:
Hi
.
The "nasty rectangles" are the output of the layout.show() function.
This function draws a simple diagram (consisting of nasty rectangles) to
indicate the regions that a ca
Thank you for your responses, Sean, Professor Ripley, Liaw.
(B
(B Sorry for my thanks.
(B
(B I'll try to use "SOM" in both "class" and "(Gene)som".
(B
(BHowever, I couldn't know how to draw 2D hexagon maps using these packages.
(B
(B Could you give any suggestion?
(B
(B Thanks.
(B
(B___
Earlier, I had posted the following question to the group :
> Hello. I have come across a curious result that I cannot explain.
> Hopefully, someone can explain this. I am doing a 1-way ANOVA with 6
> groups (example: summary(aov(y~A)) with A having 6 levels). I get an
> F of 0.899 with 5
Thanks, I appreciate knowing that.
abind() can currently take a fractional value for along, and behaves as per
your description of 'catenation' in APL.
Does APL supply any hints as to what sort of value to give 'along' to tell
abind() to perform 'lamination'?
-- Tony Plate
At Tuesday 01:22 PM
Please ignore this email. I just reread the abind
documentation and it already does this!
--- On Tue 10/21, Gabor Grothendieck < [EMAIL PROTECTED] > wrote:
From: Gabor Grothendieck [mailto: [EMAIL PROTECTED]
To: [EMAIL PROTECTED], [EMAIL PROTECTED], [EMAIL PROTECTED]
Date: Tue, 21 Oct 2003 1
I suggest following APL as that is a well thought out system.
In APL terms there are two operations here called:
- catenation. In abind, this occurs when along = 1,2,...,length(dim)
- lamination. In abind, this occurs when along = length(dim) + 1
however, the latter is really only one case of
On 21-Oct-03 Bill Shipley wrote:
> Hello. I have come across a curious result that I cannot explain.
> Hopefully, someone can explain this. I am doing a 1-way ANOVA with 6
> groups (example: summary(aov(y~A)) with A having 6 levels). I get an F
> of 0.899 with 5 and 15 df (p=0.51). I then do th
Dear Anders and Paulo,
Thanks very much for your recommendations!
I did it this way:
test2 <- unlist(lapply(test, function(x)
cbind(x$par[1], x$par[2], x$value, x$conv)))
m<- as.data.frame(matrix(test2, nrow = dim(grid)[1], ncol = 4, byrow = T))
names(m) <- c("s2", "rho", "value
Lucas Gonzalez Santa Cruz wrote:
Hi
I'd like to use R in epidemiology and disease surveillance.
In EpiInfo you can have a script (.pgm) which calls a predefined report
(.rpt), where a table is calculated and values picked from that table
and placed where the author of the report wants them, with
"Bill Shipley" <[EMAIL PROTECTED]> writes:
> Hello. I have come across a curious result that I cannot explain.
> Hopefully, someone can explain this. I am doing a 1-way ANOVA with 6
> groups (example: summary(aov(y~A)) with A having 6 levels). I get an F
> of 0.899 with 5 and 15 df (p=0.51). I
The ANOVA assumes equal variances in the groups. Suppose groups 5 and 6
had much lower variances than groups 1 to 4, and group 6 had a different
mean from the other 5 (which were about equal)?
Given how small the groups apperat to be, this could happen.
On Tue, 21 Oct 2003, Bill Shipley wrote:
> do.call("abind" c(list.of.arrays, list(along=4)))
This reminds me that I had been meaning to submit an enhancement of abind()
that allows the first argument to be a list of arrays so that you could
simply do abind(list.of.arrays, along=4), as I find this is a very common
pattern.
-- Tony Pla
I've also been thinking about how to specify that 'along' should be
length(dim)+1. At the moment one can specify any number from 0 up to
length(dim)+1, but as you point out you have to spell out length(dim)+1 as
the value for the along argument. It would possible to make abind()
automatically
Try using lapply()
For instance like:
val<-unlist(lapply(test, function(x)x$value))
You can also extend this by having your function return
everything you need from the list.
Cheers,
Anders.
On Tue, 21 Oct 2003, Tord Snall wrote:
> Dear all,
> I try extracting information from a list with
Hi everyone
I've been playing with do.call() but I'm having problems understanding it.
I have a list of "n" elements, each one of which is "d" dimensional
[actually an n-by-n-by ... by-n array]. Neither n nor d is known in
advance. I want to bind the elements together in a higher-dimensional
ar
Hello. I have come across a curious result that I cannot explain.
Hopefully, someone can explain this. I am doing a 1-way ANOVA with 6
groups (example: summary(aov(y~A)) with A having 6 levels). I get an F
of 0.899 with 5 and 15 df (p=0.51). I then do the same analysis but
using data only corre
Hi,
The solution to my problem is to use
lattice:::calculateAxisComponents to calculate appropriate labels
for the time axis in trellis plots.
# For example, given
x <- seq.POSIXt(strptime("2003/01/01", format = "%Y/%m/%d"),
strptime("2003/10/01", format = "%Y/%m/%d"), by = "mon
Dear all,
I try extracting information from a list with several levels, but I would
be happy for recommendation on writing more efficient code:
> h0<- seq(0,100, by = 20); expo<- seq(0.1, 0.5, l = 5)
> grid<- expand.grid(h0, expo)
> test<- apply(grid, 1, pcp, point.data = as.points(dat[,c("x","y"
For what it is worth, I would have thought that expressing
the lag coefficients in a B-spline expansion would be preferable
to going back to Almon approach. This would give a relatively
simple lm() application.
url:www.econ.uiuc.edu/~roger/my.htmlRoger Koenker
email [EMAIL PROTECTED
On Tue, 21 Oct 2003, Spencer Graves wrote:
> Have you checked "www.r-project.org" -> search -> "R site search"? I
> just got 15 hits for "polynomial lag". If you haven't already tried
> this, I'd guess that some of these hits (though certainly not all) might
> help you.
>
Only one of these is a
I suspect you are looking for segments().
Andy
> From: GWIGGNER Claus-Peter (EXT)
>
> Hello,
>
> Given x1, ..., xn and y1, ..., yn I'd like to draw n lines
> between xi,yi. The xi, yi shoulfd be 2-D coordinates.
>
> What is an elegant solution?
> Thanks.
>
>
>
> This message and any
Hello,
Given x1, ..., xn and y1, ..., yn I'd like to draw n lines between xi,yi.
The xi, yi shoulfd be 2-D coordinates.
What is an elegant solution?
Thanks.
This message and any files transmitted with it are legally p...{{dropped}}
__
[EMAIL PR
Hi R-Users
Is there any package in R that work Indicator Kriging?
Or somebody is working in a link between GSLIB library in
FORTRAN for R?
Thank you very much for your help.
--
__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listi
Christophe:
I suggest moving this discussion from r-help to r-sig-geo, there will be a
higher density of people there who may have had successful experiences.
There are so many different configuration issues that might impact this,
that r-help is too broad a forum. I have also CC-ed this to the
Dear Karth,
The Anova function in the car package can calculate "type-III" sums of
squares, though it doesn't do so by default. Be careful with the contrast
coding or you will get nonsense (and you might want to think about whether
you really want type-III SSs).
I hope that this helps,
John
A
> From: Rajarshi Guha [mailto:[EMAIL PROTECTED]
>
> Hi,
> I was trying to use the nnet library and am not sure of
> whats going on. I am calling the nnet function as:
>
> n <- nnet(x,y,size=3,subset=sets[[1]], maxit=200)
Please give us output of something like:
str(x)
summary(y)
Also, I b
On 10/21/03 14:38, Christoph Bier wrote:
>Jonathan Baron schrieb:
>> A very elementary overview like this is in our "Notes on R for
>> psychology experiments and questionnaires," in CRAN "contributed
>> documents" and in my R page below. We expanded it a bit from the
>> even-more elementary versio
Peter Dalgaard wrote:
On Tue, 21 Oct 2003, Christoph Bier wrote:
[...]
RTMFs are welcome =/. But I read 'help(plot)' (plot is
what I actually use for the graphic above¹) and 'help(par)',
(Read the muckin' *what*?? ;-) )
Oops :-D
[...]
(Odd, BTW, www.springer.de says it ships within 3 days
Dear R experts,
I am trying to use R with Arcgis Desktop 8.1.
When i try to add a "StatConnectorGraphicsDevice" control to my form, VBA returns an
specified error and nothing else.
Does anybody has a successful experience with Arcgis and R ?
Thanks,
Christophe Saint-Jean.
__
Jonathan Baron schrieb:
On 10/21/03 12:22, Christoph Bier wrote:
Hi,
is there an graphics overview, where the graphic capabitlities
of R are shown with the corresponding code?
A very elementary overview like this is in our "Notes on R for
psychology experiments and questionnaires," in CRAN "co
Prof Brian Ripley <[EMAIL PROTECTED]> writes:
> On Tue, 21 Oct 2003, Christoph Bier wrote:
>
> > Hi,
> >
> > is there an graphics overview, where the graphic capabitlities
> > of R are shown with the corresponding code? I already tested
> > 'demo(graphics)', that isn't that comprehensive,
> >
Prof Brian Ripley wrote:
Chapter 4 of MASS (the book) is a pretty comprehensive set of examples,
but given that there are lots of plots associated with e.g. multivariate
analysis (try chapter 11 of MASS) and time series (try chapter 14 of MASS)
the scope is enormous.
I ordered it right now in o
Hi, Lucas
You should try Sweave in the 'tools' package (see
http://www.ci.tuwien.ac.at/~leisch/Sweave/).
You will have to get a TeX/LaTeX distribution and learn a little of LaTeX
but it is worth the effort.
I frequently use R with Sweave on EpiData files (http://www.epidata.dk/)
with great succe
On Tue, 21 Oct 2003, Vittorio wrote:
> Paul Murrell [r-help] <20/10/03 09:13 +1300>:
> > Hi
> >.
> > The "nasty rectangles" are the output of the layout.show() function.
> > This function draws a simple diagram (consisting of nasty rectangles)
HI All:
Can anyone give me the formulae/steps for calculating the type III sum of
squares for an unbalanced 2-way ANOVA design? Eg. we are looking at 8
treatments x 4 doses, with unequal numbers of replications within the
groups. I really need the stepwise calculation, as I would try to put it
i
Paul Murrell [r-help] <20/10/03 09:13 +1300>:
> Hi
>.
> The "nasty rectangles" are the output of the layout.show() function.
> This function draws a simple diagram (consisting of nasty rectangles) to
> indicate the regions that a call to layout
We are not allowed to submit job directly, so I never type R to use R, just
make a batch. How can I use try() to correct my codes? In the interactive
mode, I know how to continue, but now I never enter the R window, where to
find my results and save seed to continue?
From: Jason Turner <[EMAI
It is probably a bug: does it happen when you copy from the screen in png?
I would expect not, hence that may be a workaround for you.
When I have both time and access to a Windows machine I may be able it take
a closer look: meanwhile you do have access and have the source code so
please investig
On 10/21/03 12:22, Christoph Bier wrote:
>Hi,
>
>is there an graphics overview, where the graphic capabitlities
>of R are shown with the corresponding code?
A very elementary overview like this is in our "Notes on R for
psychology experiments and questionnaires," in CRAN "contributed
documents" a
On 21-Oct-03 Michele Grassi wrote:
> Hi.
> 1)I have two variables: call a<-c(e.g.0,3,6,7...)
>b<-c(e.g.6,8,3,4...)
> I want to create a third vector z wich contain the
> pairs values z<-c(0,6,3,8,6,3,7,4and so on for each
> pairs (a,b)).
> There is a specific funct
Dear R list,
I do encounter the following problem by generating either a png-file
(example below) or a jpeg-file:
By employing 'box(which="outer")' a box is drawn, except for the right line.
If I generate the plot without the 'box(which="outer")', a line at the
bottom in the graphics file still a
Hi
I'd like to use R in epidemiology and disease surveillance.
In EpiInfo you can have a script (.pgm) which calls a predefined report
(.rpt), where a table is calculated and values picked from that table
and placed where the author of the report wants them, with text around
those values. (Please
On Tue, 21 Oct 2003, Christoph Bier wrote:
> Hi,
>
> is there an graphics overview, where the graphic capabitlities
> of R are shown with the corresponding code? I already tested
> 'demo(graphics)', that isn't that comprehensive,
> 'demo(image)', 'demo(lattice)', searched the Mailarchive,
> g
Hi,
is there an graphics overview, where the graphic capabitlities
of R are shown with the corresponding code? I already tested
'demo(graphics)', that isn't that comprehensive,
'demo(image)', 'demo(lattice)', searched the Mailarchive,
googled and the FAQ keeps silent, too.
For example, I kn
Michele Grassi wrote:
Hi.
1)I have two variables: call a<-c(e.g.0,3,6,7...)
b<-c(e.g.6,8,3,4...)
I want to create a third vector z wich contain the
pairs values z<-c(0,6,3,8,6,3,7,4and so on for each
pairs (a,b)).
There is a specific function?
How can i write my own
Michele Grassi wrote:
Hi.
1)I have two variables: call a<-c(e.g.0,3,6,7...)
b<-c(e.g.6,8,3,4...)
I want to create a third vector z wich contain the
pairs values z<-c(0,6,3,8,6,3,7,4and so on for each
pairs (a,b)).
There is a specific function?
How can i write my ow
Hi.
1)I have two variables: call a<-c(e.g.0,3,6,7...)
b<-c(e.g.6,8,3,4...)
I want to create a third vector z wich contain the
pairs values z<-c(0,6,3,8,6,3,7,4and so on for each
pairs (a,b)).
There is a specific function?
How can i write my own function?
2)When i t
Zhen Pang wrote:
I want to do 200 simulations. I found during the fisrt 128 simulations,
some parameters may be NAs, since I use if (abs(aold-anew)<1e-5) {print
(anew) break} to break the one estimation.
...
> I want to know how to resume my program with the seeds saved, and do
like continuein
Have you checked "www.r-project.org" -> search -> "R site search"? I
just got 15 hits for "polynomial lag". If you haven't already tried
this, I'd guess that some of these hits (though certainly not all) might
help you.
hope this helps. spencer graves
Francisco Vergara wrote:
Does anybody
Barnet Wagman wrote:
Do you know where I find a patched version, or where I find the patch an
instructions on how to install it? (I didn't see anything about this on
CRAN.) I'm running R-base-1.8.0-1.i386.rpm (the most recent binary
available for SuSE) and it appears to have the 'valueClass'
Dear list,
Is there a way to obtain a matrix of weighted
correlations in the presence of missing values ?
cor() can deal with NAs but cov.wt() apparently
can't. Is there any package that offers such a
function, e.g. one that uses all complete pairs
of observations ?
Thanks in advance,
Tobias
> Dear all,
>
> i have the question about the using of Jarque-Bera Test by
> using R. The question is that I do not have in my package
> "ts" this test and can not obtain any information in the
> help-file. Could you help my? Where could I download the
> package and which one, to use the Jarq
Dear all,
i have the question about the using of Jarque-Bera Test by using R. The question is
that I do not have in my package "ts" this test and can not obtain any information in
the help-file. Could you help my? Where could I download the package and which one, to
use the Jarque-Bera Test?
> "Deepayan" == Deepayan Sarkar <[EMAIL PROTECTED]>
> on Mon, 20 Oct 2003 12:58:36 -0500 writes:
Deepayan> On Monday 20 October 2003 12:35,
Deepayan> [EMAIL PROTECTED] wrote:
<.>
>> Previously, I used 'calculateAxisComponents' to massage
>> the labels manually
72 matches
Mail list logo
