[R] need suggestion about building formual

2005-09-28 Thread Simple
hi, I'm an newbie for R,I want do some fitting in R. I wander if it is possible to write a few of equations but only one formual when fitting Currently,My problem is,in R, is there methods combination a few equations into one formual? For example, y=f1(k); k=f2(t); t=f3(x); although it is

[R] gfortran Makefile for cygwin

2005-09-28 Thread Joel Bremson
Hi all, I'm porting a package that I've worked on for OS X to Cygwin/Windows. This package requires a Makefile. My question is, how can I find out (or what is), the link command? Here is the OS X Makefile: RLIB_LOC=${R_HOME} F90_FILES=\ class_data_frame.f90 \ class_old_dbest.f90 \

Re: [R] scatterplot3d + density() + polygon(color)

2005-09-28 Thread Uwe Ligges
klebyn wrote: Hi R Users, How to use the function polygon() together with the package scatterplot3d? I am trying to color below of the curves defined for the function density(). I tried to use the site: R GRAPH GALLERY as tutorial. I tried to adapt the example of this page:

Re: [R] gfortran Makefile for cygwin

2005-09-28 Thread Uwe Ligges
Joel Bremson wrote: Hi all, I'm porting a package that I've worked on for OS X to Cygwin/Windows. Cygwin is not supported. Please use MinGW's compilers, the tools and import from MkRules, as the manuals indicate. See other packages' Makefile.win files as examples. Messages regarding

Re: [R] Question on lm(): When does R-squared come out as NA?

2005-09-28 Thread Prof Brian Ripley
I've not seen a reply to this, nor ever seen it. Please make a reproducible example available (do see the posting guide). On Sun, 25 Sep 2005, Ajay Narottam Shah wrote: I have a situation with a large dataset (3000+ observations), where I'm doing lags as regressors, where I get: Call:

Re: [R] R CMD build produces tar error under FreeBSD 5.4

2005-09-28 Thread Prof Brian Ripley
On Tue, 27 Sep 2005, Eric van Gyzen wrote: On Sun, Sep 25, 2005 at 12:07:02PM +0100, Prof Brian Ripley wrote: On Sun, 25 Sep 2005, Andrew Robinson wrote: Hi R-helpers, I am trying to build a package under FreeBSD 5.4-RELEASE #0 using R Version 2.1.1. I have constructed a package using

Re: [R] scatterplot3d + density() + polygon(color)

2005-09-28 Thread Martin Maechler
UweL == Uwe Ligges [EMAIL PROTECTED] on Wed, 28 Sep 2005 08:58:16 +0200 writes: UweL klebyn wrote: Hi R Users, How to use the function polygon() together with the package scatterplot3d? I am trying to color below of the curves defined for the

[R] R2WinBUGS: Comparison to WinBUGS

2005-09-28 Thread Hadassa Brunschwig
Hi R-Help! I used R2WinBUGS and WinBUGS directly on the same model just to compare. It seems I am still making a mistake: after running the function bugs() I tried to plot the posteriors of the parameters by using read.bugs() to convert the output to an mcmc object and then plot.mcmc() to plot

Re: [R] [Fwd: R2WinBUGS: Comparison to WinBUGS]

2005-09-28 Thread Sibylle Sturtz
This is due to the following: In bugs(), the default for thinning is n.thin = max(1, floor(n.chains * (n.iter - n.burnin)/1000)) which is 29 for n.iter=12000 and n.burnin=2001 as in your example. Therefore, the number of iterations used for calculation of posterior values is (12000-2001)/29

Re: [R] scatterplot3d + density() + polygon(color)

2005-09-28 Thread Romain Francois
Le 28.09.2005 09:45, Martin Maechler a écrit : UweL == Uwe Ligges [EMAIL PROTECTED] on Wed, 28 Sep 2005 08:58:16 +0200 writes: UweL klebyn wrote: Hi R Users, How to use the function polygon() together with the package scatterplot3d? I am

Re: [R] Dummy quesion about environment

2005-09-28 Thread Ron Ophir
Thank you Peter, for the comprehensive explanation. The reason I asked Does 'search do it?' is that as I can run ls(env=environment(h)) I can run ls(env=environment(package:methods)) or ls(package:methods) which I can see by search. I thought maybe what I see by search is all the environments

Re: [R] [Fwd: R2WinBUGS: Comparison to WinBUGS]

2005-09-28 Thread Hadassa Brunschwig
Thanks for the tip. That was already helpful. But I am still not satisfied with the results. I now really changed n.thin to the same I had in WinBUGS. It looks like the commands should now be the same. However, I still get differences of 0.6 in the means of interesting parameters which should not

[R] installation on mac os x

2005-09-28 Thread richard
hello everybody I'm currently doing an internship where i need to build a pipeline between a database and R, programmed in python. For the interface between python and R i installed the RPY package. but i still can't make a connection with R because on the mac os x system i can't install R

Re: [R] Error in make check-all

2005-09-28 Thread Peter Dalgaard
Fernando Mayer [EMAIL PROTECTED] writes: Dear Prof. Ripley, that was exactly what i did. Should i login on the system as root, instead of login as user and became root via console? What should i do to have X11 usable? I don't actually think that this is the problem. You shouldn't generally

Re: [R] anova on binomial LMER objects

2005-09-28 Thread Robert Bagchi
Hi Patrick thanks for your advice. I have now tried glmmPQL, and it worked fine - I'm getting consistent results between plots and models fitted by glmmPQL. Plus it allows predict() and resid() which is another advantage over lmer at present. quick question though: why does one need to use

Re: [R] anova on binomial LMER objects

2005-09-28 Thread Prof Brian Ripley
On Wed, 28 Sep 2005, Robert Bagchi wrote: Hi Patrick thanks for your advice. I have now tried glmmPQL, and it worked fine - I'm getting consistent results between plots and models fitted by glmmPQL. Plus it allows predict() and resid() which is another advantage over lmer at present. quick

Re: [R] p-level in packages mgcv and gam

2005-09-28 Thread Denis Chabot
I only got one reply to my message: No, this won't work. The problem is the usual one with model selection: the p-value is calculated as if the df had been fixed, when really it was estimated. It is likely to be quite hard to get an honest p-value out of something that does adaptive

Re: [R] Question on lm(): When does R-squared come out as NA?

2005-09-28 Thread Ajay Narottam Shah
On Wed, Sep 28, 2005 at 08:23:59AM +0100, Prof Brian Ripley wrote: I've not seen a reply to this, nor ever seen it. Please make a reproducible example available (do see the posting guide). It was a mistake on my part. Just in case others are able to recognise the situation, what was going on

[R] boxplot and xlim confusion?

2005-09-28 Thread Karin Lagesen
I have some code as shown below. Basically, I would like three boxplots to be set next to each other with no ylabels on the two inner plots, and I want the same x axis range on all three. However, it seems like boxplot does not respect the xlim setting. I've tried the various ways I thought would

[R] Random Forest with R

2005-09-28 Thread Louis Ferre
Hi, what is the name of the package that provides Random Forest with R. Sincerely Louis Ferré http://www.univ-tlse2.fr/grimm/smash/ferre/index.html Equipe GRIMM-2254 Département de Math-Info 5 allées Antonio Machado 31058 Toulouse Cedex Tel: 0561504608 0561503982 [[alternative HTML

[R] multidimensional integration

2005-09-28 Thread Marcel Prokopczuk
dear all, i have the following problem: i want to integrate a two-dimensional function. unfortunately R crashes when i try to use adapt() and i get a nice windows message with some hex-code. do anybody of you knows how to avoid this or knows another more stable function than adapt() thanks a lot

[R] confidence variability bands for kernel estimators

2005-09-28 Thread Michael Gälger
I'm using nonparametric regression (packeges ksmooth and ks). My question: is there any way to compute confidence bands (or variability bands) with R. Confidence bands for functions are intervals [CLO(x);CUP(x)] such that with probability 1-alpha the true curve is covered by the band

[R] [R-pkgs] package 'ltm' -- version: 0.3-0

2005-09-28 Thread Dimitris Rizopoulos
Dear R users, I'd like to announce the new version of the package ltm (available from CRAN), for fitting Latent Trait Models (including the Rasch and two-parameter logistic models) under the Item Response Theory approach. Three main extra features have been added: (i) now both ltm() and

Re: [R] installation on mac os x

2005-09-28 Thread Seth Falcon
On 28 Sep 2005, [EMAIL PROTECTED] wrote: hello everybody I'm currently doing an internship where i need to build a pipeline between a database and R, programmed in python. For the interface between python and R i installed the RPY package. but i still can't make a connection with R because

Re: [R] Random Forest with R

2005-09-28 Thread TEMPL Matthias
Have a look at http://cran.r-project.org/search.html or http://cran.at.r-project.org/src/contrib/PACKAGES.html -- randomForest Best, Matthias Hi, what is the name of the package that provides Random Forest with R. Sincerely Louis Ferré

Re: [R] Random Forest with R

2005-09-28 Thread Uwe Ligges
Louis Ferre wrote: Hi, what is the name of the package that provides Random Forest with R. Really, what about looking yourself on CRAN or just googling for it? Most surprisingly the name of the package is randomForest. Uwe Ligges Sincerely Louis Ferré

Re: [R] Install and load packages

2005-09-28 Thread Spencer Graves
In Rgui 2.1.1 patched, the fourth drop-down menu after File, Edit and Misc is Packages. The last item on that menu is Install package(s) from local zip files Alternatively, RSiteSearch(install from a local zip file) produced 172 hits, some of which described other ways to do

Re: [R] Random Forest with R

2005-09-28 Thread Wiener, Matthew
It's randomForest. Searching (simple text find) on the packages web page of CRAN using either random or forest would find you this. Hope this helps, Matt Wiener -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Louis Ferre Sent: Tuesday, September 27,

[R] Access to particular predict value

2005-09-28 Thread Toni Viúdez
Hi everybody: I just generate interpolation maps with differents methods, like IDW and kriging, and now i want to compare the predict values versus real, and i don't who is the commant to do it. I want for example, if I pass from console the coordinates of a place, return me the predict value.

Re: [R] graphics guide?

2005-09-28 Thread Don MacQueen
At the R website, CRAN, in the Manuals section, you can download the document titled An Introduction to R, which contains a substantial section on the basics of R graphics. -Don At 3:27 PM +0200 9/27/05, Karin Lagesen wrote: I am trying to create some graphs with R and it seems to be able to

Re: [R] confidence variability bands for kernel estimators

2005-09-28 Thread Liaw, Andy
The `sm' package has functionality to produce pointwise band. The `locfit' package can produce simultaneous band. Both are support software for books, where you'll find more detail. Andy From: Michael Gälger I'm using nonparametric regression (packeges ksmooth and ks). My question: is

[R] Errors in odbcConnectExcel()

2005-09-28 Thread 田中 聡
Dear R-help I would like to read Excel Spreadsheets using odbcConnectExcel() in RODBC, but data in the first row can not be read. For example, I tried to read Excel file 'Book1.xls' in the current Work Directory with the following data (Range(A1:B5) in Excel), 1 19 2 27 3 61 4 76 5 98 My

Re: [R] Errors in odbcConnectExcel()

2005-09-28 Thread Prof Brian Ripley
This is entirely a function of the database format the ODBC Excel driver expects. It is *not* an error in odbcConnectExcel as your subject accuses. To use a spreadsheet as a database you need to have column names. On Thu, 29 Sep 2005, [ISO-2022-JP] ÅÄÃæ Áï wrote: Dear R-help I would like

Re: [R] p-level in packages mgcv and gam

2005-09-28 Thread Denis Chabot
Hi Yves, Le 05-09-28 à 11:05, Yves Magliulo a écrit : hi, i'll try to help you, i send a mail about this subject last week... and i did not have any response... Sorry, I did not see your message last week. I'm using gam from package mgcv. 1) How to interpret the significance of smooth

Re: [R] p-level in packages mgcv and gam

2005-09-28 Thread Thomas Lumley
On Wed, 28 Sep 2005, Denis Chabot wrote: I only got one reply to my message: No, this won't work. The problem is the usual one with model selection: the p-value is calculated as if the df had been fixed, when really it was estimated. It is likely to be quite hard to get an honest p-value

[R] Change console language ?

2005-09-28 Thread Sérgio Nunes
Hi, I'm pretty new to R, I've just installed version 2.1.1 on Windows. I have a simple (it seems) doubt - how do I change the Console default Language ? It's set to Portuguese but I would like to view it in English. Thanks in advance, Sérgio Nunes __

Re: [R] p-level in packages mgcv and gam

2005-09-28 Thread Peter Dalgaard
Thomas Lumley [EMAIL PROTECTED] writes: Bob, on the other hand, chooses the amount of smoothing depending on the data. When a 4 df smooth fits best he ends up with the same model as Alice and the same p-value. When some other df fits best he ends up with a different model and a *smaller*

Re: [R] Change console language ?

2005-09-28 Thread Prof Brian Ripley
On Wed, 28 Sep 2005, Sérgio Nunes wrote: Hi, I'm pretty new to R, I've just installed version 2.1.1 on Windows. I have a simple (it seems) doubt - how do I change the Console default Language ? It's set to Portuguese but I would like to view it in English. Do read the administration manual

[R] gee models summary

2005-09-28 Thread Dean Sonneborn
I'm running some GEE models but when I request the summary(pcb.gee) all I get are rows and rows of intercorelations and they fill up the screen buffer so I can not even scroll back to see what else might be in the summary. How do I get the summary function to NOT print the intercorrelations?

Re: [R] Error in make check-all

2005-09-28 Thread Fernando Mayer
Dear Peter Dalgaard, first, sorry for my delay to this reply, but I only have acces to this machine in the morning (consider also the time it takes to run make and make check)... It seems that the problem was really with some missing -devel packages. I followed your suggestions and installed

Re: [R] gee models summary

2005-09-28 Thread Renaud Lancelot
This inconvenience does not occur with function geese in package geepack. If you're using function gee in package gee, extract the components you want from the fitted object: library(gee) library(MASS) data(OME) fm - gee(cbind(Correct, Trials-Correct) ~ Loud + Age + OME, id = ID, +

[R] responses to my question on non-central t

2005-09-28 Thread Robert Kinley
The original question is attached at the end of this message ... In the world of R the excellent mvtnorm package answers my original question and many more that I didn't know I needed to ask. In the world of Splus ( version 6 , Windows XP ) this link :-

Re: [R] p-level in packages mgcv and gam

2005-09-28 Thread Denis Chabot
But what about another analogy, that of polynomials? You may not be sure what degree polynomial to use, and you have not decided before analysing your data. You fit different polynomials to your data, checking if added degrees increase r2 sufficiently by doing F-tests. I thought it was the

[R] dumping plots

2005-09-28 Thread kthomps7
i've heard that you can dump a plot file, and in doing so set a pixel parameter (cbx?). is anyone familar with this procedure? i would appreciate any insight. thanks, kevin __ R-help@stat.math.ethz.ch mailing list

Re: [R] dumping plots

2005-09-28 Thread Pierce, Ken
There's the xpinch and ypinch parameters in the windows() device call. However, as an addition to this question, how do you specify output, say with the savePlot function, to maximize the plot area on an 8 1/2 x 11 sheet of paper? Ken ~~~ Kenneth B. Pierce Jr. Research

Re: [R] p-level in packages mgcv and gam

2005-09-28 Thread Liaw, Andy
Just change the df in what Thomas described to the degree of polynomial, and everything he said still applies. Any good book on regression that covers polynomial regression ought to point this out. Andy From: Denis Chabot But what about another analogy, that of polynomials? You may not be

Re: [R] p-level in packages mgcv and gam

2005-09-28 Thread Thomas Lumley
On Wed, 28 Sep 2005, Denis Chabot wrote: But what about another analogy, that of polynomials? You may not be sure what degree polynomial to use, and you have not decided before analysing your data. You fit different polynomials to your data, checking if added degrees increase r2

[R] xyplots

2005-09-28 Thread Nathan Leon Pace, MD, MStat
Hi All, I have a four panel xyplot. I wish to plot each point as an open or filled circle depending on the value of an indicator variable. I assume I need to use panel.superpose(), but I can't figure out the syntax from lattice documentation. Running R 2.1 under Mac OS X 10.4.2. Any

[R] Fast AUC computation

2005-09-28 Thread Nina Paynter
I am doing a simulation with a relatively large data set (20,000 observations) for which I want to calculate the area under the Receiver Operator Curve (AUC) for many parameter combinations. I am using the ROC library and the following commands to generate each AUC: rocobj=rocdemo.sca(truth =

Re: [R] xyplots

2005-09-28 Thread Deepayan Sarkar
On 9/28/05, Nathan Leon Pace, MD, MStat [EMAIL PROTECTED] wrote: Hi All, I have a four panel xyplot. I wish to plot each point as an open or filled circle depending on the value of an indicator variable. I assume I need to use panel.superpose(), but I can't figure out the syntax from

[R] correct syntax

2005-09-28 Thread Stephen Choularton
Hi I am trying to produce a little table like this: 0 1 0 601 408 1 290 2655 but I cannot get the syntax right. Can anyone help. Stephen -- 27/09/2005 [[alternative HTML version deleted]]

Re: [R] is it possible to form matrix of matrices...and multiple arrays

2005-09-28 Thread venomousanimal
booop booop schrieb: Dear sirs, 1...Kindly tell me is it possible to form a matrix which contains a no of matrices.. for eg.. if a,b,c,d are matrices and e is a matrix which contains a,b,c,d as rows and columns.. 2..Is it possible to form array of array of arrays for eg..

[R] R-code for binormla distribution

2005-09-28 Thread Nabil Channouf
Dear users, does any one have a code (S or R) to compute the binormal distribution (or the upper its quadrant area) other than the pmvnorm. Thanks -- Nabil Channouf etudiant en Ph.D. Bureau 3733 Departement d'Informatique et de Recherche Operationnelle (D.I.R.O.) Universite de Montreal, C.P.

[R] xyplots

2005-09-28 Thread Nathan Leon Pace, MD, MStat
Thanks for suggestions on plotting different symbols in an xyplot by Sarkar and Wiener in answer to my question: I have a four panel xyplot. I wish to plot each point as an open or filled circle depending on the value of an indicator variable. Wiener: You don't even need to tell it to use

Re: [R] correct syntax

2005-09-28 Thread Marc Schwartz (via MN)
On Thu, 2005-09-29 at 06:57 +1000, Stephen Choularton wrote: Hi I am trying to produce a little table like this: 0 1 0 601 408 1 290 2655 but I cannot get the syntax right. Can anyone help. Stephen Do you have the raw

Re: [R] correct syntax

2005-09-28 Thread venomousanimal
Stephen Choularton schrieb: Hi I am trying to produce a little table like this: 0 1 0 601 408 1 290 2655 but I cannot get the syntax right. Can anyone help. Stephen vec1= c(0,0,1) vec2= c(NA,601,290) vec3= c(1,408,2655)

[R] After resizing a tkwidget, how do I get the size of the widget??

2005-09-28 Thread Søren Højsgaard
Dear all, I am trying make a small tcltk thing for drawing graphs with vertices, edges etc. I can draw the graph in a window. - I would like it so that if I resize the window, the graph will also be resized. To do this, I guess I need the size of the window, so that I can calculate the new

[R] Problem with memory footprint of qq plot generated with lattice

2005-09-28 Thread nwew
Dear R helpers, I generate a qq plot using the following function call. qqmath(~val|ind,data=xx ,distribution=function(p) qt(p,df=19) ,ylab=Sample Quatinles ,xlab=Theoretical Quantiles ,aspect=1 ,prepanel = prepanel.qqmathline ,panel=function(x,y)

Re: [R] Plot Data Points in boxplots

2005-09-28 Thread Aric Gregson
--On 9/27/05 8:55 PM -0400 Gabor Grothendieck sent: Try this: boxplot(Sepal.Length ~ Species, iris) with(iris, stripchart(Sepal.Length ~ Species, vertical = TRUE, add = TRUE)) Thanks very much for the hint. I did something very similar with decent results. The only problem I am having now

Re: [R] Plot Data Points in boxplots

2005-09-28 Thread Paul Murrell
Hi Aric Gregson wrote: --On 9/27/05 8:55 PM -0400 Gabor Grothendieck sent: Try this: boxplot(Sepal.Length ~ Species, iris) with(iris, stripchart(Sepal.Length ~ Species, vertical = TRUE, add = TRUE)) Thanks very much for the hint. I did something very similar with decent results. The

Re: [R] Problem with memory footprint of qq plot generated with lattice

2005-09-28 Thread dhinds
nwew [EMAIL PROTECTED] wrote: Dear R helpers, I generate a qq plot using the following function call. ... dim(xx) [1] 680237 2 How about doing something like this: fn - function(n,cut=0.001,m=1000) { p - ppoints(n) p - p[pmin(p, 1-p) cut] q -

Re: [R] Plot Data Points in boxplots

2005-09-28 Thread Jim Lemon
Aric Gregson wrote: Hello, I would like to plot my data in a fashion similar to a boxplot, but plot the true data points without a box, just overlay lines of the summary generated by the boxplot. I have less than 10 observations within each group, so I think showing the actual data would

[R] memory issues with large data set

2005-09-28 Thread Christina Yau
Hi, I am running R 2.0.1.1. on Windows. It is a Dell Dimension with a 3.2 Ghz Processor and 4Gb RAM. When using the ReadAffy() function to read in 97 arrays, I get the below error messages: Error: cannot allocate vector of size 393529 Reached total allocation of 1024Mb: see

Re: [R] memory issues with large data set

2005-09-28 Thread James W. MacDonald
Christina Yau wrote: Hi, I am running R 2.0.1.1. on Windows. It is a Dell Dimension with a 3.2 Ghz Processor and 4Gb RAM. This question concerns a BioC package, so the correct listserv is [EMAIL PROTECTED], not the R-help listserv. In the future, you should direct questions about BioC

Re: [R] is it possible to form matrix of matrices...and multiple arrays

2005-09-28 Thread Gabor Grothendieck
On 9/28/05, booop booop [EMAIL PROTECTED] wrote: Dear sirs, 1...Kindly tell me is it possible to form a matrix which contains a no of matrices.. for eg.. if a,b,c,d are matrices and e is a matrix which contains a,b,c,d as rows and columns.. Try this: mat - list(diag(2),

[R] priceIts

2005-09-28 Thread Remigijus Lapinskas
Dear All, There is an example for the priceIts function (the its package) which does not work for me as expected. ?priceIts x1 - priceIts(instrument = c(^ftse), start = 1998-01-01, + quote = Close) Error in validObject(.Object) : invalid class its object: Missing

Re: [R] priceIts

2005-09-28 Thread Mulholland, Tom
Well I downloaded the data using the link in your message which suggests that the code is right. I don't have its loaded (I assume it's from the irregular time series package) so I can't test the code as you have it. Have you tried checking to see it it is an issue with the internet (your

Re: [R] priceIts

2005-09-28 Thread Gabor Grothendieck
The error message indicates that a problem is occurring in download.file. Try to use download.file to download some web page from the internet to see if that works, download.file(http://www.google.com;, destfile = google.htm) If that does NOT work then examine the download.file help page