[R] NNet value and convergence

2003-09-30 Thread Hannah Wood
Hi, I'm using the R nnet package and have a few simple (?) questions. What is the value that is output after every 10 iterations during the training of the network and how is it calculated? # weights: 177 initial value 506.134586 iter 10 value 128.222774 iter 20 value 95.399782 iter

[R] FW: error predicting values from the LME

2003-09-30 Thread Andrej Kveder
HI all, I might add some more information in order to possibly solve my problem. I'm really stuck and no obvious solutions do the trick. I'm using R 1.7.1 on Windows 2000 with the packages regurarly updated. I'm using hypothetical data constructed as a pseudo population conforming to a certain

Re: [R] shared object

2003-09-30 Thread Uwe Ligges
p hu wrote: Dear R-helper, I have one question for creating shared objects (C code) on Unix. Since there is a set of softwares (such as Rtools, Perl, etc.) to be dowloaded for generating shared objects on windows, I am wondering whether I need to install a set of these kind of softwares to

Re: [R] shared object

2003-09-30 Thread Martin Maechler
Hu == p hu [EMAIL PROTECTED] on Mon, 29 Sep 2003 23:38:24 -0400 (EDT) writes: Hu Dear R-helper, I have one question for creating shared Hu objects (C code) on Unix. Since there is a set of Hu softwares (such as Rtools, Perl, etc.) to be dowloaded Hu for generating shared

[R] using identify() together with plot () and pixmap()

2003-09-30 Thread Christoph Lehmann
Dear R users I have a two-dimensional array, whose values I want to plot, using the pixmapGrey class. Plotting works fine, and now I would like to be able to identify some of the points in the plot using identify(). But I get the following message while pressing the left mouse button:

[R] Biclassification

2003-09-30 Thread Dowkiw, Arnaud
Is there any R package/function for biclassification (hierarchical reciprocal clustering of 2 factors) ? Thanks, Arnaud DISCLAIMER**...{{dropped}} __ [EMAIL PROTECTED] mailing list

Re: [R] using identify() together with plot () and pixmap()

2003-09-30 Thread Barry Rowlingson
Christoph Lehmann wrote: plot(pixmapGrey(fmri.vtc[,,slice,volume])) identify(fmri.vtc[,,slice,volume]) warning: no point with 0.25 inches pressing the right mouse button I get: numeric(0) what is the problem here and how can I solve it? The problem is that there is no method for identifying

Re: [R] using identify() together with plot () and pixmap()

2003-09-30 Thread Friedrich . Leisch
On Tue, 30 Sep 2003 10:16:11 +0200, Christoph Lehmann (CL) wrote: Dear R users I have a two-dimensional array, whose values I want to plot, using the pixmapGrey class. Plotting works fine, and now I would like to be able to identify some of the points in the plot using identify().

Re: [R] BIC or AIC from nnet

2003-09-30 Thread Uwe Ligges
Paul Green wrote: Is AIC or BIC available when using the nnet package? Given you are talking about multinom(): Yes, it returns the AIC for the fit, but that's already mentioned in the help file... Uwe Ligges Thank you Paul Green __ [EMAIL

Re: [R] using identify() together with plot () and pixmap()

2003-09-30 Thread Friedrich . Leisch
On Tue, 30 Sep 2003 09:38:31 +0100, Barry Rowlingson (BR) wrote: Christoph Lehmann wrote: plot(pixmapGrey(fmri.vtc[,,slice,volume])) identify(fmri.vtc[,,slice,volume]) warning: no point with 0.25 inches pressing the right mouse button I get: numeric(0) what

Re: [R] NNet value and convergence

2003-09-30 Thread Uwe Ligges
Hannah Wood wrote: Hi, I'm using the R nnet package and have a few simple (?) questions. What is the value that is output after every 10 iterations during the training of the network and how is it calculated? # weights: 177 initial value 506.134586 iter 10 value 128.222774 iter 20

Re: [R] CP for rpart

2003-09-30 Thread Uwe Ligges
weidong zhang wrote: Hi All, I have some questions on using library rpart. Given my data below, the plotcp gives me increasing 'xerrors' across different cp's with huge xstd (plot attached). What causes the problem or it's not a problem at all? I am thinking 'xerror's should be decreasing

[R] cluster mgcv update

2003-09-30 Thread Martin Wegmann
Hello, After reinstalling the whole OS and R as well, I tried to update.packages() and get the follwing error message: concerning the mgcv update: atlas2-base is installed and blas as well (on debian). I haven't found lf77blas, I assume it's a library or something similar associated with

Re: [R] cluster mgcv update

2003-09-30 Thread Uwe Ligges
Martin Wegmann wrote: Hello, After reinstalling the whole OS and R as well, I tried to update.packages() and get the follwing error message: concerning the mgcv update: atlas2-base is installed and blas as well (on debian). I haven't found lf77blas, I assume it's a library or something

Re: [R] cluster mgcv update

2003-09-30 Thread Martin Wegmann
You need to add atlas2-base-dev: $ apt-get install atlas2-base-dev I installed atlas2-base-dev and g77 but know I get the error messages pasted below. Both (cluster and mgcv) requires lfrtbegin, but that does not seem to be programm which I can install via apt-get. Martin * Installing

Re: [R] overlay two pixmap

2003-09-30 Thread Christoph Lehmann
when I try to overlay a completely transparent pixmap on another pixmap, I get an error. For reproduction: just the transparent pixmap itself gives an error: tmp - array(0,c(x.dim,y.dim)) tmp - pixmapIndexed(tmp[,]) for (x in 1:x.dim) { for (y in 1:y.dim) { [EMAIL PROTECTED],y] - NA } }

Re: [R] checking generic/method consistency

2003-09-30 Thread Arne Henningsen
Hi *, thanks for all your answers and discussions. And additionally special thanks to Henrik Bengtsson for writing the first draft of the R Coding Conventions. I think that this document contains a lot of good ideas to make the code more readable. Since my package is *new* code, I adjusted it

Re: [R] cluster mgcv update

2003-09-30 Thread Dirk Eddelbuettel
On Tue, Sep 30, 2003 at 02:04:23PM +0200, Martin Wegmann wrote: You need to add atlas2-base-dev: $ apt-get install atlas2-base-dev I installed atlas2-base-dev and g77 but know I get the error messages pasted below. Both (cluster and mgcv) requires lfrtbegin, but that does not seem to

[R] non-linear trends in kriging model

2003-09-30 Thread David Pleydell
Hi I am struggling to fit a non-linear trend using the likfit function in geoR. Specifically I want a sigmoidal function, something like SSfpl in the nls package to fit the trend. But it seems trend.spatial in geoR only works with lm or glm type models. Any ideas how I can specify the model to

Re: [R] non-linear trends in kriging model

2003-09-30 Thread Paulo Justiniano Ribeiro Jr
David Indeed, the trend term in geoR must be a linear one. You can: 1) fit you model by alternating between using nls() and passing the residuals to the functions in geoR or 2) try trying to fit using the functions in the nlme() package Regards P.J. On Tue, 30 Sep 2003, David Pleydell wrote:

Re: [R] FW: error predicting values from the LME

2003-09-30 Thread Spencer Graves
Could you dumb it down to a toy example with 4 observations for a model like y ~ 1 | inter, 2 observations for each of 2 levels of inter? If that works, then you can play with the example that works and the example that doesn't; this is one of the strategies mentioned in Poly (1971) How

Re: [R] debugging R

2003-09-30 Thread Douglas Bates
Francisco J Molina [EMAIL PROTECTED] writes: To debug R from within Xemacs I have tried two different things: C-u M-x R RET - d SPC gdb RET to start an inferior R process with arguments @option{-d gdb} This is the way described in the official documentation of R. and also M-x R

Re: [R] checking generic/method consistency

2003-09-30 Thread Martin Maechler
Arne == Arne Henningsen [EMAIL PROTECTED] on Tue, 30 Sep 2003 14:18:11 +0200 writes: Arne Hi *, Arne thanks for all your answers and discussions. And Arne additionally special thanks to Henrik Bengtsson for Arne writing the first draft of the R Coding Arne Conventions.

[R] Stepwise procedures

2003-09-30 Thread MZodet
Is there a function in R which performs stepwise estimation in ways similar to SAS/STATA (i.e., allows the analyst to specify the significance levels for removal/addition of terms). I've been asked to evaluate two final models: one resulting from a backwards selection in R (stepAIC) and one

Re: [R] checking generic/method consistency

2003-09-30 Thread Arne Henningsen
On Tuesday 30 September 2003 15:23, Martin Maechler wrote: Arne == Arne Henningsen [EMAIL PROTECTED] on Tue, 30 Sep 2003 14:18:11 +0200 writes: Arne Hi *, Arne thanks for all your answers and discussions. And Arne additionally special thanks to Henrik Bengtsson for

[R] Generation of labels for 2 factors

2003-09-30 Thread Xavier Fernández i Marín
Hello, I have a data set with 10 countries, 10 sectors of public policy and 15 years like country sector year a aa 1980 a aa 1981 ... a bb 1980 a bb 1981 ... b aa 1980 b aa 1981 ... b bb 1980 b bb 1981

Re: [R] Stepwise procedures

2003-09-30 Thread Frank E Harrell Jr
On Tue, 30 Sep 2003 09:24:20 -0400 [EMAIL PROTECTED] wrote: Is there a function in R which performs stepwise estimation in ways similar to SAS/STATA (i.e., allows the analyst to specify the significance levels for removal/addition of terms). I've been asked to evaluate two final

[R] ts - unit conversion

2003-09-30 Thread michaell taylor
I've been using R for a while, but now find myself needing to understand time-series objects for a short course I am teaching. I am putting together some daily financial datasets for illustration, but having some trouble in aggregating the data to months or weeks. I am getting a cannot change

[R] Grouped bar plots

2003-09-30 Thread Matjaz Kukar
Dear R experts! I'm wondering if it is possible in R to make a certain type of bar plots, in S-Plus known as Grouped Bar Plots (see the attachment for an example). In S-Plus they can be created with the guiPlot function. Thank you for any answers. Best regards, Matjaz.inline:

Re: [R] Grouped bar plots

2003-09-30 Thread Uwe Ligges
Matjaz Kukar wrote: Dear R experts! I'm wondering if it is possible in R to make a certain type of bar plots, in S-Plus known as Grouped Bar Plots (see the attachment for an example). In S-Plus they can be created with the guiPlot function. Thank you for any answers. ?barplot Uwe Ligges Best

Re: [R] R Production Performance

2003-09-30 Thread Joe Conway
Zitan Broth wrote: Right but R is only preloaded once? Yes. The plr shared library gets loaded and initialized (which in turn loads and initializes libR) only once -- on Postgres's postmaster startup. From that point forward, every new database connection gets a forked copy of the postmaster,

RE: [R] Generation of labels for 2 factors

2003-09-30 Thread Liaw, Andy
Not sure exactly what you want, but try: countrySector - factor(paste(df$country, df$sector, sep=:)) where df is the data frame. This creates a factor for the combination. If you want numeric code, just coerce with as.numeric(). HTH, Andy -Original Message- From: Xavier Fernández

RE: [R] ts - unit conversion

2003-09-30 Thread Heywood, Giles
Financial data, as you point out, is generally irregular, which is in essence what prompted the devlopment of the irregular time-series (its) package, which is posted on CRAN (v0.1.2 posted today, incidentally). In the its class, the time-stamps of the rows of a matrix are represented using the

RE: [R] ts - unit conversion

2003-09-30 Thread Gabor Grothendieck
Change the line in get.hist.quote that pastes the url together so that g=d (i.e. daily) is replaced with g=m (i.e. monthly). That will cause monthly data to be downloaded. With this new get.hist.quote try this: ibm - get.hist.quote(ibm) ibm - ts( ibm[!is.na(rowSums(ibm)),], end=c(2003,8),

[R] Remove comma (,) in data set

2003-09-30 Thread Muhammad Subianto
Dear R-helper, I am new learning R. Now, I have a data set like: 24,2,3,3,1,1,2,3,0,1 45,1,3,10,1,1,3,4,0,1 43,2,3,7,1,1,3,4,0,1 42,3,2,9,1,1,3,3,0,1 36,3,3,8,1,1,3,2,0,1 19,4,4,0,1,1,3,3,0,1 38,2,3,6,1,1,3,2,0,1 21,3,3,1,1,0,3,2,0,1 27,2,3,3,1,1,3,4,0,1 45,1,1,8,1,1,2,2,1,1 ... with 3730 rows I

Re: [R] Remove comma (,) in data set

2003-09-30 Thread Duncan Murdoch
On Tue, 30 Sep 2003 18:51:52 +0200, Muhammad Subianto [EMAIL PROTECTED] wrote : Dear R-helper, I am new learning R. Now, I have a data set like: 24,2,3,3,1,1,2,3,0,1 45,1,3,10,1,1,3,4,0,1 ... with 3730 rows I want to remove comma (,) in data set. The result like: 24 2 3 3 1 1 2 3 0 1 45 1 3

[R] subsetting a matrix

2003-09-30 Thread Rajarshi Guha
Hi, I'm trying to take a set of rows and columns out of a matrix. I hve been using the index aray approach. My overll matrix is X and is 179 x 65. I want to take out 4 columns and 161 rows. Thus I made a 161 x 2 array I and filled it up with the row,col indices. However doing, X[ I ] gives me

[R] lme vs. aov

2003-09-30 Thread array chip
Hi, I have a question about using lme and aov for the following dataset. If I understand correctly, using aov with Error term in the formula is equivalent to using lme with default settings, i.e. both assume compound symmetry correlation structure. And I have found that equivalency in the past.

Re: [R] subsetting a matrix

2003-09-30 Thread Thomas W Blackwell
Rajarshi - Why not simply subscript your matrix X to return the rows and columns you want to keep ? For example, new - X[16:176, c(3,5,7,9)] assuming those are the rows and columns you want. See help(Extract). - tom blackwell - u michigan medical school - ann arbor - On Tue, 30 Sep

[R] generic function object

2003-09-30 Thread Mathieu Drapeau
Hi, I would like to know if there is a way to see the code of a generic function? Like when I only type the function name, it prints the source of that function. But I cannot see the code of a generic function, it prints: ... function (object, ...) standardGeneric(normalize) environment:

Re: [R] Preserving names from p to f in nlm (was: can't get names from vector in nlm calls)

2003-09-30 Thread Spencer Graves
Perhaps someone else can help you figure out how to preserve the names from argument p to function f in calls to nlm, but I'm now beyond the limits of my knowledge. I'm sorry, but I can't help further with this. Spencer Graves Jeff D. Hamann wrote: I started looking

Re: [R] generic function object

2003-09-30 Thread James MacDonald
getMethods(normalize) should do the trick. Best, Jim James W. MacDonald Affymetrix and cDNA Microarray Core University of Michigan Cancer Center 1500 E. Medical Center Drive 7410 CCGC Ann Arbor MI 48109 734-647-5623 Mathieu Drapeau [EMAIL PROTECTED] 09/30/03 02:14PM Hi, I would like to

RE: [R] subsetting a matrix

2003-09-30 Thread Ted Harding
On 30-Sep-03 Rajarshi Guha wrote: Hi, I'm trying to take a set of rows and columns out of a matrix. I hve been using the index aray approach. My overll matrix is X and is 179 x 65. I want to take out 4 columns and 161 rows. Thus I made a 161 x 2 array I and filled it up with the row,col

Re: [R] cluster mgcv update

2003-09-30 Thread Dirk Eddelbuettel
On Tue, Sep 30, 2003 at 02:04:23PM +0200, Martin Wegmann wrote: You need to add atlas2-base-dev: $ apt-get install atlas2-base-dev I installed atlas2-base-dev and g77 but know I get the error messages pasted below. Both (cluster and mgcv) requires lfrtbegin, but that does not seem to

[R] fitdistr, mle's and gamma distribution

2003-09-30 Thread Lourens Olivier Walters
Dear R Users, I am trying to obtain a best-fit analytic distribution for a dataset with 11535459 entries. The data range in value from 1 to 3. I use: fitdistr(data, gamma) to obtain mle's for the parameters. I get the following error: Error in optim(start, mylogfn, x = x, hessian =

Re: [R] subsetting a matrix

2003-09-30 Thread Jason Turner
(Ted Harding) wrote: On 30-Sep-03 Rajarshi Guha wrote: Hi, I'm trying to take a set of rows and columns out of a matrix. I hve been using the index aray approach. My overll matrix is X and is 179 x 65. I want to take out 4 columns and 161 rows. ... This is documented in An Introduction to R,

[R] Problem with R1.8beta

2003-09-30 Thread Phguardiol
Hi, I trying to use r.180beta on a XPPro PC with Athlon 2.0 I downloaded the zip file and when I unzipped it using winzip last version I had this message: C:\Program Files\R\rw1080beta\doc\manual\refman.pdf An error occured while trying to copy a file: The source file is corrupted. Anything to do

Re: [R] fitdistr, mle's and gamma distribution

2003-09-30 Thread Spencer Graves
In my experience, the most likely cause of this problem is that optim may try to test nonpositive values for shape or scale. I avoid this situation by programming the log(likelihood) in terms of log(shape) and log(scale) as follows: gammaLoglik - + function(x, logShape, logScale,

Re: [R] subsetting a matrix

2003-09-30 Thread Rajarshi Guha
On Tue, 2003-09-30 at 15:51, Jason Turner wrote: (Ted Harding) wrote: On 30-Sep-03 Rajarshi Guha wrote: Hi, I'm trying to take a set of rows and columns out of a matrix. I hve been using the index aray approach. My overll matrix is X and is 179 x 65. I want to take out 4 columns and

[R] in addition R1.8b install

2003-09-30 Thread Phguardiol
this problem also occurred later in the installation process: C:\Program Files\R\rw1080beta\doc\manual\refman.pdf C:\Program Files\R\rw1080beta\library\mva\html\loadings.html C:\Program Files\R\rw1080beta\library\ctest\html\fisher.test.html C:\Program Files\R\rw1080beta\library\base\chtml\base.chm

Re: [R] in addition R1.8b install

2003-09-30 Thread Marc Schwartz
On Tue, 2003-09-30 at 15:00, [EMAIL PROTECTED] wrote: this problem also occurred later in the installation process: C:\Program Files\R\rw1080beta\doc\manual\refman.pdf C:\Program Files\R\rw1080beta\library\mva\html\loadings.html C:\Program Files\R\rw1080beta\library\ctest\html\fisher.test.html

Re: [R] in addition R1.8b install

2003-09-30 Thread Marc Schwartz
On Tue, 2003-09-30 at 15:35, Marc Schwartz wrote: On Tue, 2003-09-30 at 15:00, [EMAIL PROTECTED] wrote: this problem also occurred later in the installation process: C:\Program Files\R\rw1080beta\doc\manual\refman.pdf C:\Program Files\R\rw1080beta\library\mva\html\loadings.html C:\Program

Re: [R] fitdistr, mle's and gamma distribution

2003-09-30 Thread Ben Bolker
Spencer Graves's suggestion of using shape and scale parameters on a log scale is a good one. To do specifically what you want (check values for which the objective function is called and see what happens) you can do the following (untested!), which makes a local copy of dgamma that you

[R] truncated multivariate normal

2003-09-30 Thread Juliana Garcia Cespedes
Please, I would like to know how to generate a truncated multivariate normal distribution k - dimensional,X ~ NT(mu, Sigma), where the elements of X to be non-negative (except the first), and the first dimension is strictly larger than zero. Example: X ~ NT_2(mu, Sigma),

Re: [R] fitdistr, mle's and gamma distribution

2003-09-30 Thread Ben Bolker
PS. 11 MILLION entries?? On Tue, 30 Sep 2003, Ben Bolker wrote: Spencer Graves's suggestion of using shape and scale parameters on a log scale is a good one. To do specifically what you want (check values for which the objective function is called and see what happens) you can

[R] Re: Help needed: plotting with no device

2003-09-30 Thread Ross Boylan
Thanks to everyone for your help. I decided to see if the session could be recovered if I connected back from the original, local terminal. The local screen was locked by the KDE screensaver. Either my unlocking it, or the passage of time, seems to have got the process unstuck. It happened

Re: [R] truncated multivariate normal

2003-09-30 Thread Ben Bolker
Do you mean you want to generate deviates from this distribution or that you want the density function (or distribution function)? mvrnorm in package MASS (library(MASS)) will generate multivariate normal deviates. A brute-force approach would just keep picking values until you get ones in

Re: [R] Help needed: plotting with no device

2003-09-30 Thread Ross Boylan
On Mon, 2003-09-29 at 14:25, Luke Tierney wrote: Look at ?Signal and see if that would help. You may need to rename the appropriate .RData beforehand to be safer. luke For those who follow, note it's actually ?Signals you want. It talks about USR1 and USR2 (as did some other people in this

[R] variables

2003-09-30 Thread forkusam
Hi, can someone please help me. I will give a simple example of my problem. p - function() { i - 1 sr - function(){ i-i+3 i- sqrt(i) } cat(i) } This is just an example. My main problem is defining i like a global variable which I can use in the sub- and main functions without any

Re: [R] Problem with R1.8beta

2003-09-30 Thread Duncan Murdoch
On Tue, 30 Sep 2003 15:56:03 EDT, you wrote: Hi, I trying to use r.180beta on a XPPro PC with Athlon 2.0 I downloaded the zip file and when I unzipped it using winzip last version I had this message: C:\Program Files\R\rw1080beta\doc\manual\refman.pdf An error occured while trying to copy a file:

Re: [R] variables

2003-09-30 Thread Jason Turner
forkusam wrote: Hi, can someone please help me. I will give a simple example of my problem. p - function() { i - 1 sr - function(){ i-i+3 i- sqrt(i) } cat(i) } This is just an example. My main problem is defining i like a global variable which I can use in the sub- and main functions

[R] R icon for RedHat 9.

2003-09-30 Thread Carlos J. Gil Bellosta
Dear R-helpers, I just installed R on a RedHat 9 machine and when I was trying add a launcher on the panel to it, I could not find any R icon in png format. Is there any available? Do I have it and I could not find it in my file hierarchy? Thank you for your help. Carlos J. Gil Bellosta

[R] R-icon for RH Linux 9.

2003-09-30 Thread Carlos J. Gil Bellosta
Dear R helpers, It is the first time I have to install R under Linux (Red Hat 9) and as I was setting up the panels, I could not find any png icon to launch the program. Does it exist and I have overlooked it? Can I get it somewhere? Carlos J. Gil Bellosta Sigma Consultores Estadísticos

Re: [R] Remove comma (,) in data set

2003-09-30 Thread Philippe Glaziou
Muhammad Subianto [EMAIL PROTECTED] wrote: I am new learning R. Now, I have a data set like: 24,2,3,3,1,1,2,3,0,1 45,1,3,10,1,1,3,4,0,1 ... with 3730 rows I want to remove comma (,) in data set. The result like: 24 2 3 3 1 1 2 3 0 1 45 1 3 10 1 1 3 4 0 1 ... How can I do it.