On 20 Oct 2003, Rajarshi Guha wrote:
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)
Where x is a 272x4 matrix of observations (examples) and y is a 272x1
matrix of target
On Tue, 21 Oct 2003, Paul Livingstone wrote:
I am using R 1.6.1 on Windows. I usually write/edit code in the R console,
and when I get it working right, copy and paste it to a text file for later
use. In Splus for Windows, this could be done easily using just the
keyboard, with up-arrow,
Farouk Nathoo [EMAIL PROTECTED] writes:
parameter - 0 # initial value
for (i in 1:1000)
{
parameter - one.step(parameter, data)
mem - memory.size()
cat(parameter, , mem,\n)
}
I output the
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 but that function
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?
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
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,
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'
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 know
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 continueing
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 try
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 own
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
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,
googled
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.
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
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 function?
How
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 and in
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
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
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() has
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
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) to
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
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
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,
'demo(image)',
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
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.
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
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 version that was
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 believe
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
At
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
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
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
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 files
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 about
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
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)]),
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 = month)
y -
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
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
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
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()
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 Plate
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
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 then
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,
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, conv)
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 the
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
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
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
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
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
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 call
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
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) - value
It's a
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,
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 bind the
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 [EMAIL
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
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
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
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 scoping,
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
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
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
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
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
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
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
window,
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 the
74 matches
Mail list logo