Re: [R] Lattice device page options-margins

2005-03-10 Thread Paul Murrell
Hi Bock, Michael wrote: -Original Message- From: Deepayan Sarkar [mailto:[EMAIL PROTECTED] Sent: Wednesday, March 09, 2005 5:09 PM To: r-help@stat.math.ethz.ch Cc: Paul Murrell; Bock, Michael; Sundar Dorai-Raj Subject: Re: [R] Lattice device page options-margins On Wednesday 09 March

RE: [R] Lattice device page options-margins

2005-03-10 Thread Bock, Michael
-Original Message- From: Deepayan Sarkar [mailto:[EMAIL PROTECTED] Sent: Wednesday, March 09, 2005 5:09 PM To: r-help@stat.math.ethz.ch Cc: Paul Murrell; Bock, Michael; Sundar Dorai-Raj Subject: Re: [R] Lattice device page options-margins On Wednesday 09 March 2005 13:54, Paul

Re: [R] plotting

2005-03-10 Thread Marc Schwartz
On Wed, 2005-03-09 at 14:28 -0400, Owen Buchner wrote: I have two questions for you. Firstly I'm having troubles trying to plot more then 1 graph. I'm attempting to make a plot with 9 panels, but i have no clue what type of code to use. If you are using R's base graphics and you want 9 plots

Re: [R] RMySQL installed but not availalable

2005-03-10 Thread Don MacQueen
require(RMySQL) help('MySQL') then see the examples shown by the help page At 9:12 PM -1200 3/9/05, Adriano von Sydow wrote: Hi I run Linux SuSE 9.1, with MSQL 4.0.18 I installed R 2.0.1 and it is working fine I installed RM 0.5-5 package and verified all installed.packages (see below) but when

Re: [R] two-dimensional integration?

2005-03-10 Thread Prof Brian Ripley
Nils, For 2D, see package 'adapt' on CRAN. e.g. adapt(2, c(0,0), c(1,1), functn=function(x) sin(prod(x))*exp(x[1]-x[2])) Package `adapt' will do larger numbers of dimensions, but numerical quadrature is often no more effective than Monte-Carlo methods in more than a few dimensions. For very

Re: [R] two-dimensional integration?

2005-03-10 Thread Dimitris Rizopoulos
This existed even before the Athens Olympics 2004 :) look at package adapt which integrates a function from 2 to 20 dimensions, e.g., library(adapt) f - function(x){ x1 - x[1] x2 - x[2] sin(x1*x2)*exp(x1-x2) } adapt(2, c(0,0), c(1,1), functn=f) I hope it helps. Best, Dimitris

Re: [R] two-dimensional integration?

2005-03-10 Thread Peter Dalgaard
Nils-at-Duke Lid Hjort [EMAIL PROTECTED] writes: I find the one-dimensional integrate very helpful, but often enough I stumble into problems that require two (or more)-dimensional integrals. I suppose there ... (c) Will a thing like the above exist in R before the Tromsoe Olympics in 2014?

[R] Dropping coloumns while redaing dtaa from text file.

2005-03-10 Thread Upasna Sharma
Hi I have a huge text file and .dat file from which I want to read data. I do not need all the columns in the files. I want to extract only some columns from the .txt file or the .dat file, because reading the entire file is becoming very difficult due to memory issues. Is it possible to extract

[R] R: LIST function and LOOPS

2005-03-10 Thread Clark Allan
hi all another simple question. i've written a dummy program so that you get the concept. (the code could be simplfied such that there are no loops. but lets leave the loops in for now.) z1-function(w) { for (i in 1:w) { set.seed(i+6) ss-0 for (j in 1:5) {

Re: [R] contrast matrix for aov

2005-03-10 Thread Prof Brian Ripley
On Wed, 9 Mar 2005, Darren Weber wrote: How do we specify a contrast interaction matrix for an ANOVA model? We have a two-factor, repeated measures design, with Where does `repeated measures' come into this? You appear to have repeated a 2x2 experiment in each of 8 blocks (subjects). Such a

Re: [R] Dropping coloumns while redaing dtaa from text file.

2005-03-10 Thread Peter Dalgaard
Upasna Sharma [EMAIL PROTECTED] writes: Hi I have a huge text file and .dat file from which I want to read data. I do not need all the columns in the files. I want to extract only some columns from the .txt file or the .dat file, because reading the entire file is becoming very difficult

Re: [R] Dropping coloumns while redaing dtaa from text file.

2005-03-10 Thread Prof Brian Ripley
On Thu, 10 Mar 2005, Upasna Sharma wrote: Hi I have a huge text file and .dat file from which I want to read data. I do not need all the columns in the files. I want to extract only some columns from the .txt file or the .dat file, because reading the entire file is becoming very difficult due to

Re: [R] function in order to plot the same graph to postscript and pdf

2005-03-10 Thread Ronny Klein
I guess you would want eval(MYPLOT) there? But the function above works if you say plot.both(myplot = expression(plot(x)), filename = foo) and that is what I said above: set myplot = expression(plot(x)). Z My syntax in the example was a little bit messed, sorry for that. However Gabor

Re: [R] Flattening a list of data frames

2005-03-10 Thread Dimitris Rizopoulos
try this: do.call(rbind, z) Best, Dimitris Dimitris Rizopoulos Ph.D. Student Biostatistical Centre School of Public Health Catholic University of Leuven Address: Kapucijnenvoer 35, Leuven, Belgium Tel: +32/16/336899 Fax: +32/16/337015 Web: http://www.med.kuleuven.ac.be/biostat/

[R] Flattening a list of data frames

2005-03-10 Thread Lorin Hochstein
Hello all, Simple version of my problem: I've got a list of data frames, where each data frame has the same number of columns and the same column names. I'd like to flatten the list into one large data frame. Is there an easy way to do this? Quick example code: a -

Re: [R] How to get standard deviation of rows in a matrix

2005-03-10 Thread Sundar Dorai-Raj
Jagarlamudi, Choudary wrote on 3/9/2005 1:49 PM: Hi all, I am trying to find sd of my rows in a matrix and i get column sd inspite of extracting rows. I tried to do the sqrt(var(x)) but that did'nt work as well, Here is my data genes 15 24 63 40 25 42 46 35 23 53 37 45 30 37 50 55 40 51

Re: [R] contrast matrix for aov

2005-03-10 Thread Peter Dalgaard
Prof Brian Ripley [EMAIL PROTECTED] writes: On Wed, 9 Mar 2005, Darren Weber wrote: How do we specify a contrast interaction matrix for an ANOVA model? We have a two-factor, repeated measures design, with Where does `repeated measures' come into this? You appear to have repeated a

[R] Interval censoring in Survival analysis

2005-03-10 Thread AMOROS NAVARRO, ALEX
Hi, I would like to know if R programme allows doing different kinds of censoring in his last version. What kind of package should I use and is this censoring applicable to parametric and semi-parametric(Cox) models? Thanks, Àlex. __

Re: [R] contrast matrix for aov

2005-03-10 Thread Christophe Pallier
Prof Brian Ripley wrote: On Wed, 9 Mar 2005, Darren Weber wrote: We have a two-factor, repeated measures design, with Where does `repeated measures' come into this? You appear to have repeated a 2x2 experiment in each of 8 blocks (subjects). Such a design is usually analysed with fixed

Re: [R] from long/lat to UTM

2005-03-10 Thread Ted Harding
On 10-Mar-05 Sander Oom wrote: Hi Yyan, The proj4R package by Roger Bivand will allow you to project data in many ways and directions. http://spatial.nhh.no/R/Devel/proj4R-pkg.pdf It uses the proj libraries from: http://www.remotesensing.org/proj/ Not sure where you would derive

Re: [R] How to get standard deviation of rows in a matrix

2005-03-10 Thread Tobias Verbeke
On Wed, 9 Mar 2005 13:49:44 -0600 Jagarlamudi, Choudary [EMAIL PROTECTED] wrote: Hi all, I am trying to find sd of my rows in a matrix and i get column sd inspite of extracting rows. I tried to do the sqrt(var(x)) but that did'nt work as well, Here is my data genes 15 24 63 40

Re: [R] contrast matrix for aov

2005-03-10 Thread Prof Brian Ripley
On Thu, 10 Mar 2005, Christophe Pallier wrote: Prof Brian Ripley wrote: On Wed, 9 Mar 2005, Darren Weber wrote: We have a two-factor, repeated measures design, with Where does `repeated measures' come into this? You appear to have repeated a 2x2 experiment in each of 8 blocks (subjects). Such

Re: [R] contrast matrix for aov

2005-03-10 Thread Prof Brian Ripley
On Thu, 10 Mar 2005, Peter Dalgaard wrote: Prof Brian Ripley [EMAIL PROTECTED] writes: On Wed, 9 Mar 2005, Darren Weber wrote: How do we specify a contrast interaction matrix for an ANOVA model? We have a two-factor, repeated measures design, with Where does `repeated measures' come into this?

RE: [R] How to get standard deviation of rows in a matrix

2005-03-10 Thread Liaw, Andy
You should read ?var and ?sd more carefully. For a data frame or a matrix, var() returns the covariance matrix of the columns, whereas sd() returns the standard deviations of the columns. If you want standard deviations of the rows, you need to transpose the data. Andy From: Jagarlamudi,

RE: [R] Interval censoring in Survival analysis

2005-03-10 Thread Rau, Roland
Hi, -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of AMOROS NAVARRO, ALEX I would like to know if R programme allows doing different kinds of censoring in his last version. What kind of package should I use Yes, try the following:

Re: [R] Lattice device page options-margins

2005-03-10 Thread Chuck Cleland
You can add some space to the left by setting the left.padding component of the layout.widths trellis setting. For example: library(lattice) trellis.device(new=TRUE) trellis.par.set(layout.widths = list(left.padding = 20)) require(stats) ## Tonga Trench Earthquakes Depth -

RE: [R] Structural equation models with R

2005-03-10 Thread André TavaresCorrêa Dias
Dear John, Thanks very much! I don’t know how I didn’t see this before. I checked a thousand of times… There is one more thing that I can not understand. What are the possible reasons for problems on calculation of confidence interval for RMSEA? For the same model I sent before (after

[R] multiple comparisons for lme using multcomp

2005-03-10 Thread Michaël Coeurdassier
Dear R-help list, I would like to perform multiple comparisons for lme. Can you report to me if my way to is correct or not? Please, note that I am not nor a statistician nor a mathematician, so, some understandings are sometimes quite hard for me. According to the previous helps on the

[R] Re: plotting several series on the sames axes

2005-03-10 Thread midnightsun
Many thanks to all those who replied to my 'simple' query, all of which have been posted to the list. __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide!

Re: [R] about kpss.test()

2005-03-10 Thread Achim Zeileis
On Wed, 9 Mar 2005 14:32:09 -0500 (EST) Weiguang Shi wrote: --- Achim Zeileis [EMAIL PROTECTED] wrote: First of all, could you tell me what the KPSS Level in the output of the test means? Thanks. But I meant to also ask about the number after that, e.g., KPSS Level = 0.0027

Re: [R] problem using uniroot with integrate

2005-03-10 Thread Ken Knoblauch
Thank you very much. Yes, that was the problem, partial matching. I saw a warning about that in integrate and some discussion from 1999 in the archives and so added the m0 for integrate but somehow I wasn't bright enough to see that I had the same problem in uniroot. Sorry about no working

Re: [R] multiple comparisons for lme using multcomp

2005-03-10 Thread Peter Dalgaard
Michaël Coeurdassier [EMAIL PROTECTED] writes: summary(csimtest(vect,vcov(lm1),cmatrix=contrMat(table(treatment),type=Tukey),df=59)) Coefficients: Estimate t value Std.Err. p raw p Bonf p adj Al800-Al100 -2.253 -10.4670.213 0.000 0.000 0.000 Al600-Al100

[R] R-help:loading data from text file and viewing it in persp

2005-03-10 Thread Lakshmi Dhevi Baskar
Hi All i would like to get some help on how to use the persp() for 3d surface plotting with the data from txt file. it has format like x1 y1 z1 x2 y2 z2 x3 y3 z3 and so on...sometimes the txt file data may be very large also.. how could i load this txt file and view its 3d

Re: [R] how modify object in parent.env

2005-03-10 Thread Gabor Grothendieck
[I sent this last night but it does not appear to have shown up so this is second attempt. Apologies if this gets posted twice.] The search for the left hand side of the - starts in the parent while search for the right side is in the current environment so: x - 1:3 local( for (i in

[R] Gap statistic

2005-03-10 Thread Nestor Fernandez
Dear All, I need to calculate the optimal number of clusters for a classification based on a large number of observations (tens of thousands). Thibshirani et al. proposed the gap statistic for this purpose. I tried the R-code developed by R. Jörnsten but R hangs with such amount of data (). Is

[R] Gregmisc

2005-03-10 Thread Christian Kamenik
Dear all, I use R 2.0.1 on Windows XP professional. When I want to load the 'Gregmisc' library I get the following error message: Error in library(pkg, character.only = TRUE) : 'gregmisc' is not a valid package -- installed 2.0.0? Can anybody tell me what's wrong with this package? Cheers,

Re: [R] Gregmisc

2005-03-10 Thread Chuck Cleland
It's a bundle of 4 packages (gtools, gdata, gmodels, and gplots), so you need to load one of those packages. See the following: http://cran.us.r-project.org/src/contrib/Descriptions/gregmisc.html Chuck Cleland Christian Kamenik wrote: Dear all, I use R 2.0.1 on Windows XP professional. When I

Re: [R] Gregmisc

2005-03-10 Thread Marc Schwartz
On Thu, 2005-03-10 at 14:44 +0100, Christian Kamenik wrote: Dear all, I use R 2.0.1 on Windows XP professional. When I want to load the 'Gregmisc' library I get the following error message: Error in library(pkg, character.only = TRUE) : 'gregmisc' is not a valid package -- installed

Re: [R] from long/lat to UTM

2005-03-10 Thread Gabor Grothendieck
Ted.Harding at nessie.mcc.ac.uk writes: : : yyan liu wrote: : Hi: :Is there any function in R which can convert the : long/lat to UTM(Universal Transverse Mercator)? :There are quite a few converters on Internet. : However, the interface is designed as input-output : which I can not

Re: [R] Gregmisc

2005-03-10 Thread Uwe Ligges
On Thu, 10 Mar 2005, Christian Kamenik wrote: Dear all, I use R 2.0.1 on Windows XP professional. When I want to load the 'Gregmisc' library I get the following error message: Error in library(pkg, character.only = TRUE) : 'gregmisc' is not a valid package -- installed 2.0.0?

RE: [R] Structural equation models with R

2005-03-10 Thread John Fox
Dear Andre, -Original Message- From: André TavaresCorrêa Dias [mailto:[EMAIL PROTECTED] Sent: Thursday, March 10, 2005 6:14 AM To: John Fox Cc: r-help@stat.math.ethz.ch Subject: RE: [R] Structural equation models with R Dear John, Thanks very much! I don’t know how I didn’t see

Re: [R] install R redhat rpm as a normal user

2005-03-10 Thread Marc Schwartz
On Wed, 2005-03-09 at 16:17 -0800, Eric Hu wrote: Hi, I wonder if anyone has done this before. I have rpm-build installed in my workstation. Thanks. Eric I have not seen any other replies to this, but as far as I know Martyn's RPMS are not relocatable and a test this morning confirms that.

Re: [R] R: LIST function and LOOPS

2005-03-10 Thread Adaikalavan Ramasamy
You will need to capture the value of ss at the end of each 'i' as such z4 -function(w){ output - numeric(w) for (i in 1:w){ set.seed(i+6) # this is redundant line ss-0 for (j in 1:5){ set.seed(j+1+(i-1)*6) r-rnorm(1) ss-ss+r } output[i] -

[R] Vogelsang test? RATS? Time series analysis.

2005-03-10 Thread abunn
Has anybody implemented the Vogelsang test from Vogelsang's 1998 Econometrica paper? Or the tests by Fomby and Vogelsang from their 2002 paper? I have been looking for trends in many time series (1000's) with unknown autocorrelation. In brief, I have fit AR models and classified them as

Re: [R] problem using uniroot with integrate

2005-03-10 Thread Thomas Lumley
On Wed, 9 Mar 2005, Ken Knoblauch wrote: Thank you very much. Yes, that was the problem, partial matching. I saw a warning about that in integrate and some discussion from 1999 in the archives and so added the m0 for integrate but somehow I wasn't bright enough to see that I had the same problem

Re: [R] Interval censoring in Survival analysis

2005-03-10 Thread Thomas Lumley
On Thu, 10 Mar 2005, AMOROS NAVARRO, ALEX wrote: Hi, I would like to know if R programme allows doing different kinds of censoring in his last version. What kind of package should I use and is this censoring applicable to parametric and semi-parametric(Cox) models? To some extent. The survival

[R] Re: [Rd] dealing with package bundles

2005-03-10 Thread Seth Falcon
Liaw, Andy [EMAIL PROTECTED] writes: Given the confusions that some users have regarding package bundles, I'd like to suggest some changes to how package bundles are handled. I wonder if a guideline of naming packages bundles in a manner that identifies them as bundles would help? An example

Re: [R] from long/lat to UTM

2005-03-10 Thread Ted Harding
On 10-Mar-05 Gabor Grothendieck wrote: Ted.Harding at nessie.mcc.ac.uk writes: : [...] : As Sander says, Not sure where you would derive the time zone!. : : Unless you can refer a (long,lat) position to a look-up : table, you can't predict what the zone will be to less : the 1 hour (except of

RE: [R] How could I catch the R data-output stream and presented by other software func

2005-03-10 Thread Warnes, Gregory R
There are quite a few packages which allow R to be accessed from other software systems. Examples include RDCOM (included with Windows R), RSOAP (http://rsoap.sf.net), Rpy (Python, http://rpy.sf.net), RSPerl (http://www.omegahat.org/RSPerl/). Duncan Temple Lang, the developer of RSPerl, has a

[R] Online short course: Advanced Resampling Methods

2005-03-10 Thread Peter C. Bruce
Online short course: Advanced Resampling Methods Phillip Good will offer his online short course “Advanced Resampling Methods” March 18 ­ April 8. In the bootstrap area, it will cover BCA (bias corrected and accelerated) confidence intervals, bootstrap-t intervals, the parametric bootstrap,

Re: [R] from long/lat to UTM

2005-03-10 Thread Greg Snow
Ted Harding [EMAIL PROTECTED] 03/10/05 09:13AM On 10-Mar-05 Gabor Grothendieck wrote: Ted.Harding at nessie.mcc.ac.uk writes: : [...] : As Sander says, Not sure where you would derive the time zone!. : [snip] The only completely general mechanism I can think of would consist of

Re: [R] from long/lat to UTM

2005-03-10 Thread Ted Harding
On 10-Mar-05 Greg Snow wrote: [...] I found a shapefile of time zones at: http://openmap.bbn.com/data/shape/timezone/ A google search on the keywords timezone shapefile finds a few other possibilities. b) a function which, for each name in the list, returns the TZ name The

Re: [R] contrast matrix for aov

2005-03-10 Thread Darren Weber
Prof Brian Ripley wrote: On Wed, 9 Mar 2005, Darren Weber wrote: How do we specify a contrast interaction matrix for an ANOVA model? We have a two-factor, repeated measures design, with Where does `repeated measures' come into this? You appear to have repeated a 2x2 experiment in each of 8

[R] about R CMD check

2005-03-10 Thread xt_wang
hello, everybody, I create a package by myself named var on Linux system. I attach a C code which uses R_ext/Lapack.h in directory src. But when I carry out R CMD check var, it shows wxt0124() Error in .C(wxt1221, PACKAGE = var) : C function name not in DLL for package var Execution halted

Re: [R] contrast matrix for aov

2005-03-10 Thread Darren Weber
As an R newbie (formerly SPSS), I was pleased to find some helpful notes on ANOVA here: http://personality-project.org/r/r.anova.html In my case, I believe the relevant section is: Example 4. Two-Way Within-Subjects ANOVA This is where I noted and copied the error notation. Sorry for any

[R] how to view the syntax of a method which is not a generic method

2005-03-10 Thread Sarah Holte
Hello - I'm trying to modify an option for the lme() or nlme() macros. I want to write my own specification for the variance function and am following homework problem 4, Chapter 5, page 268 of Pinheiro and Bates book on mixed effect. I'm up to point where I've created a new class using an

[R] RODBC autocommit

2005-03-10 Thread Omar Lakkis
How can I turn autocommit off using my RODBC connection to an Informix database? I want to turn autocommit off, insert a thousand or so rows then commit. __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do

[R] Transparent colors OR two series on one histogram

2005-03-10 Thread Nathaniel Street
Hi, I want to be able to plot a single histogram of a measured trait with trait values from two conditions on the same histogram to allow easy comparison. I have previously done this in excel by plotting the two series on a single bargraph having calculated frequencies in bins. You then get

Re: [R] how to view the syntax of a method which is not a generic method

2005-03-10 Thread Seth Falcon
Hi Sarah, Sarah Holte [EMAIL PROTECTED] writes: The book suggests that I use the inialize method for one of the exisiting variance functions for lme(), eg varExp. What I want is the syntax of the initialize.varExp method so that I can edit it to create an initialize method for my newly

Re: [R] Transparent colors OR two series on one histogram

2005-03-10 Thread Thomas Lumley
On Thu, 10 Mar 2005, Nathaniel Street wrote: Hi, I want to be able to plot a single histogram of a measured trait with trait values from two conditions on the same histogram to allow easy comparison. I have previously done this in excel by plotting the two series on a single bargraph having

[R] Question regarding mosaicplot

2005-03-10 Thread Jean Vidal
I tried this : mosaicplot(stoc ~ q9r + segca,data=tmp2,color=T) : works fine. And now, this : mosaicplot(stoc ~ q9r + segca, data=tmp2, color=T, main=Big title) Error in model.frame(formula, rownames, variables, varnames, extras, extranames, : invalid variable type I'm probably stupid

[R] Logistic regression goodness of fit tests

2005-03-10 Thread Trevor Wiens
I was unsure of what suitable goodness-of-fit tests existed in R for logistic regression. After searching the R-help archive I found that using the Design models and resid, could be used to calculate this as follows: d - datadist(mydataframe) options(datadist = 'd') fit - lrm(response ~

Re: [R] Question regarding mosaicplot

2005-03-10 Thread Achim Zeileis
Jean: On Tue, 18 Nov 2003 23:32:05 Jean Vidal wrote: I tried this : mosaicplot(stoc ~ q9r + segca,data=tmp2,color=T) : works fine. And now, this : mosaicplot(stoc ~ q9r + segca, data=tmp2, color=T, main=Big title) Error in model.frame(formula, rownames, variables, varnames, extras,

Re: [R] Logistic regression goodness of fit tests

2005-03-10 Thread Roger Bivand
On Thu, 10 Mar 2005, Trevor Wiens wrote: I was unsure of what suitable goodness-of-fit tests existed in R for logistic regression. After searching the R-help archive I found that using the Design models and resid, could be used to calculate this as follows: d - datadist(mydataframe)

Re: [R] Logistic regression goodness of fit tests

2005-03-10 Thread Frank E Harrell Jr
Trevor Wiens wrote: I was unsure of what suitable goodness-of-fit tests existed in R for logistic regression. After searching the R-help archive I found that using the Design models and resid, could be used to calculate this as follows: d - datadist(mydataframe) options(datadist = 'd') fit -

Re: [R] install R redhat rpm as a normal user

2005-03-10 Thread Eric Hu
Thank you all. Best, Eric On Mar 10, 2005, at 7:21 AM, Marc Schwartz wrote: On Wed, 2005-03-09 at 16:17 -0800, Eric Hu wrote: Hi, I wonder if anyone has done this before. I have rpm-build installed in my workstation. Thanks. Eric I have not seen any other replies to this, but as far as I know

Re: [R] Logistic regression goodness of fit tests

2005-03-10 Thread Trevor Wiens
On Thu, 10 Mar 2005 22:36:09 +0100 (CET) Roger Bivand [EMAIL PROTECTED] wrote: From one geographer to another, and being prepared to bow to better-founded explanations, you seem to have included a variable - the offending variable slr_mean - that is very highly correlated with another.

Re: [R] Logistic regression goodness of fit tests

2005-03-10 Thread Trevor Wiens
On Thu, 10 Mar 2005 16:19:41 -0600 Frank E Harrell Jr [EMAIL PROTECTED] wrote: The goodness of fit test only works on prespecified models. It is not valid when stepwise variable selection is used (unless perhaps you use alpha=0.5). Perhaps I'm blind, but I can't find any reference to an

Re: [R] Is anyone using the MiniR distribution?

2005-03-10 Thread Duncan Murdoch
On Fri, 04 Feb 2005 14:45:35 +, Duncan Murdoch [EMAIL PROTECTED] wrote : On Fri, 4 Feb 2005 14:46:03 +0100, Martin Maechler [EMAIL PROTECTED] wrote : Duncan == Duncan Murdoch [EMAIL PROTECTED] on Fri, 04 Feb 2005 09:50:22 + writes: Duncan The miniR files only include a

[R] Any upcoming R Advanced Programming in Bay area?

2005-03-10 Thread paul king
Hi all, I am taking a new job in San Francisco and wonder if there's an upcoming R advanced course in bay area. Best- Paul __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide!

[R] Setting variable main in plot-function

2005-03-10 Thread T Petersen
Hi, I am plotting the residuals of z like this plot(resid(z)) and I want the title of the graph to be main=Residuals from regression [name of regression] where in this case z is the name of the regression. Is there a way to automaticall put the name of the regression in the title, e.g if the

[R] Ploting a function of two arguments

2005-03-10 Thread Martin C. Martin
Hi, I've written a function: myfun - function(x, y) { // blah blah } and I want to graph it. My plan is to use persp(), and my question is: how do I create the array of values? One possibility is: x - seq(0, 10, by=.1) y - seq(0, 10, by=.1) inputs - somehow create an array of x,y pairs

[R] trouble loading R function code

2005-03-10 Thread David Reinke
I created an R function offline with a text editor and saved it as fhv.R When I type in load(file=fhv.R) I get the message Error: bad restore file magic number (file may be corrupted)-- no data loaded I've checked the file, and it opens up fine in Notepad and other text editors. I also tried this

Re: [R] trouble loading R function code

2005-03-10 Thread Martin C. Martin
try: source(fhv.R) David Reinke wrote: I created an R function offline with a text editor and saved it as fhv.R When I type in load(file=fhv.R) I get the message Error: bad restore file magic number (file may be corrupted)-- no data loaded I've checked the file, and it opens up fine in Notepad

[R] Errors reading data file?

2005-03-10 Thread Chris Bergstresser
Hi all -- I tried loading a data file with the following command: data = read.table(filename.txt, header = TRUE, sep = ,) This appeared to work fine, except it silently skipped 400 records (out of 1200). It turns out, some of the text fields included quotes, and I needed to use 'quote =

[R] Bonferroni simultaneous confidence intervals for multiple regression

2005-03-10 Thread Cara Gormally
Hi, I'm having no luck figuring out how to find Bonferroni simultaneous confidence intervals to obtain a family of estimates in R. Does anyone know how to do this? Thank you! __ R-help@stat.math.ethz.ch mailing list

Re: [R] Setting variable main in plot-function

2005-03-10 Thread Marc Schwartz
On Fri, 2005-03-11 at 00:50 +0100, T Petersen wrote: Hi, I am plotting the residuals of z like this plot(resid(z)) and I want the title of the graph to be main=Residuals from regression [name of regression] where in this case z is the name of the regression. Is there a way to

[R] XML to data frame or list

2005-03-10 Thread bogdan romocea
Dear useRs, I have a simple/RTFM question about XML parsing. Given an XML file, such as (fragment) A100/A B23/B Ctrue/C how do I import it in a data frame or list, so that the values (100, 23, true) can be accessed through the names A, B and C? I installed the XML package and looked over the

[R] How to reply to a posting if I subscribe to the digest version of R-help?

2005-03-10 Thread Keith Wong
Dear listmembers, Is there a recommended method of responding to a thread if I am getting the mailing list in digest form? Will the mailing list program automatically match threads if I copy the subject from the original posting? Thank you. Keith __

[R] sample function

2005-03-10 Thread mirage sell
Hi everyone, I need help. I want to have a uniform kind distribution. When I used sample function I got almost twice many zeros compared to other numbers. What's wrong with my command ? temp -sample(0:12, 2000, replace=T,prob=(rep(1/13,13))) hist(temp) Thanks in advance, Taka,

Re: [R] sample function

2005-03-10 Thread Martin C. Martin
hist is lumping things together. Try: sum(temp == 0) compare to the height of the left most bar. Is this a bug in hist? - Martin mirage sell wrote: Hi everyone, I need help. I want to have a uniform kind distribution. When I used sample function I got almost twice many zeros compared to other

RE: [R] sample function

2005-03-10 Thread Liaw, Andy
It's not the simulated data, but how hist() handled it. If you use truehist() in the MASS package, you don't see the problem. Nor would you see it like this: table(temp)/length(temp) temp 0 1 2 3 4 5 6 7 8 9 10 11 0.0745 0.0745 0.0830

Re: [R] sample function

2005-03-10 Thread David Scott
On Thu, 10 Mar 2005, mirage sell wrote: Hi everyone, I need help. I want to have a uniform kind distribution. When I used sample function I got almost twice many zeros compared to other numbers. What's wrong with my command ? Nothing is wrong with your sampling, it is the display in the

Re: [R] sample function

2005-03-10 Thread David Scott
On Thu, 10 Mar 2005, Martin C. Martin wrote: hist is lumping things together. Try: sum(temp == 0) compare to the height of the left most bar. Is this a bug in hist? No, hist is the wrong thing to use to display this data. Try temp -sample(0:12, 2000, replace=T,prob=(rep(1/13,13)))

Re: [R] sample function

2005-03-10 Thread Marc Schwartz
On Thu, 2005-03-10 at 20:54 -0600, mirage sell wrote: Hi everyone, I need help. I want to have a uniform kind distribution. When I used sample function I got almost twice many zeros compared to other numbers. What's wrong with my command ? temp -sample(0:12, 2000,

Re: [R] Interval censoring in Survival analysis

2005-03-10 Thread A.J. Rossini
On Thu, 10 Mar 2005 08:01:57 -0800 (PST), Thomas Lumley [EMAIL PROTECTED] wrote: On Thu, 10 Mar 2005, AMOROS NAVARRO, ALEX wrote: I am not aware of any package that allows semiparametric regression modelling of interval censored or doubly censored data. I'm aware of 2 such packages, but the

Apologies --- Re: [R] Question regarding mosaicplot

2005-03-10 Thread Jean Vidal
May I present my apologies to all, on this list, and to you Achim ? It's an old message that was sent by accident yesterday when testing the Thunderbird mail program. So, please, ignore this question. It was answered and solved months ago. At 22:16 10/03/2005, Achim Zeileis wrote: Jean: On

[R] help on warning message when using anova

2005-03-10 Thread tpeng
Hi can someone help me with the following warning message I got when I used two way anova: Warning messages: 1: Models with response NULL removed because response differs from model 1 in: anova.lmlist(object, ...) 2: Models with response NULL removed because response differs from model 1 in:

[R] howto: plot location, web search

2005-03-10 Thread Mike R
Hi, new to the list, R and rpy. my first post. using rpy, i find this r.X11(height=3, width=6) allows me to preset the size of a subsequent r.plot(x, y) but i've been unable to find documention on related aspects, such as programatic control of the plot window location on the monitor,

[R] xlab cex bug in Lattice/grid and multipage-output (Windows only?)

2005-03-10 Thread Dieter Menne
The following problem seems to be Windows-specific. Deepayan informed me that it cannot be reproduced on platform i386-pc-linux-gnu. - library(lattice) trellis.device(png, file=sepal%02d.png,color=T, theme=col.whitebg) p=xyplot(Sepal.Length + Sepal.Width ~ Petal.Length + Petal.Width |

RE: [R] howto: plot location, web search

2005-03-10 Thread Mulholland, Tom
If it's specifically R that I am looking for I use the following link http://finzi.psych.upenn.edu/search.html If I were going to use google I would go to the advanced page and insert finzi.psych.upenn.edu into the Domain field so as to restrict the search. As you have found out the letter R

RE: [R] howto: plot location, web search

2005-03-10 Thread Mulholland, Tom
Sorry about the wrong link. The archives are at http://sourceforge.net/mailarchive/forum.php?forum=rpy-list -Original Message- From: Mike R [mailto:[EMAIL PROTECTED] Sent: Friday, 11 March 2005 2:44 PM To: r-help@stat.math.ethz.ch Subject: [R] howto: plot location, web search

RE: [R] about R CMD check

2005-03-10 Thread Liaw, Andy
Please show us at least the prototype of the C function(s) being called, and the .C call in the R code. Just to make sure: The first argument to .C() should be the name of the C function, not the name of the C source file. Andy From: [EMAIL PROTECTED] hello, everybody, I create a