[R] arrangement of crowded labels

2009-05-20 Thread Thomas Zumbrunn
Dear R users I'm looking for algorithms that assist in spreading out crowded labels, e.g. labels of points in a scatter plot, in order to obtain a nicer visual appearance and better legibility. I'm probably just stuck because I didn't find the right key words for a successful search on the R

Re: [R] arrangement of crowded labels

2009-05-20 Thread Richard . Cotton
I'm looking for algorithms that assist in spreading out crowded labels, e.g. labels of points in a scatter plot, in order to obtain a nicer visual appearance and better legibility. I'm probably just stuck because I didn't find the right key words for a successful search on the R

Re: [R] nlrwr package. Error when fitting the optimal Box-Cox transformation with two variables

2009-05-20 Thread F.Tusell
El mar, 19-05-2009 a las 17:17 +0200, Ikerne del Valle escribió: Dear all: I'm trying to fit the optimal Box-Cox transformation related to nls (see the code below) for the demand of money data in Green (3th Edition) but in the last step R gives the next error message.

[R] Running multiple versions of Rserve on a Linux system

2009-05-20 Thread Harsh
Hi useRs, I would like to know how to run two versions of Rserve on a Linux system. I am using the old JRE client of Rserve which requires me to use Rserve version 0.4-3. I am also using the new REngine API to connect to Rserve version 0.6. I would like to run both these Rserve versions on a

Re: [R] Package Inline under windows

2009-05-20 Thread Duncan Murdoch
_ wrote: Hi all, I installed the package inline (windows-version) but can not compile any code, I alway get an error message ERROR(s) during compilation : source code errors or compiler configuration errors! Unfornutanely there is no description where the package finds a c-compiler nor

[R] Functions returning functions

2009-05-20 Thread Paulo Grahl
Dear All: I have a question regarding the behavior of functions. Say I define a function that returns another function : A - function(parameters) { # calculations w/ parameters returning 'y' tmpf - function(x) { # function of 'y' } return(tmpf) } The value of the parameters are

Re: [R] *** POSSIBLE SPAM *** Functions returning functions

2009-05-20 Thread Romain Francois
Paulo Grahl wrote: Dear All: I have a question regarding the behavior of functions. Say I define a function that returns another function : A - function(parameters) { # calculations w/ parameters returning 'y' tmpf - function(x) { # function of 'y' } return(tmpf) } The value of

[R] SEM:Standard error of std.coef estimates?

2009-05-20 Thread BdeGroot
Hi, I am currently working with the sem package in R, to create pathway diagrams. Id like to use the standardized path coeffcients. To get these, I use std.coef. However, using this yields only the standardized coefficients, but does not give me the standard error. Does someone know how to get

Re: [R] Too large a data set to be handled by R?

2009-05-20 Thread jim holtman
If your 1500 X 2 matrix is all numeric, it should take up about 240MB of memory. That should easily fit within the 2GB of your laptop and still leave room for several copies that might arise during the processing. Exactly what are you going to be doing with the data? A lot will depend on the

Re: [R] how to calculate means of matrix elements

2009-05-20 Thread Gabor Grothendieck
In that case use a modification of Jim's solution: a - array(cbind(mat1, mat2, mat3), c(3, 3, 3)) apply(a, 1:2, mean, na.rm = TRUE) On Tue, May 19, 2009 at 8:49 AM, dxc13 dx...@health.state.ny.us wrote: Easy enough.  What if some of the matrix elements contained missing values? Then how

Re: [R] Functions returning functions

2009-05-20 Thread Wacek Kusnierczyk
Paulo Grahl wrote: Dear All: I have a question regarding the behavior of functions. Say I define a function that returns another function : A - function(parameters) { # calculations w/ parameters returning 'y' tmpf - function(x) { # function of 'y' } return(tmpf) } The

[R] Where to find a changelog for the survival package

2009-05-20 Thread Heinz Tuechler
Dear All, since some days I try to use the versions 2.35-4 of the survival package instead of versions 2.31, I had installed until now. Several changes in print.survfit, plot.survfit and seemingly in the structure of ratetabels effect some of my syntax files. Is there somewhere a

Re: [R] Functions returning functions

2009-05-20 Thread Gabor Grothendieck
On Wed, May 20, 2009 at 7:48 AM, Wacek Kusnierczyk waclaw.marcin.kusnierc...@idi.ntnu.no wrote: Paulo Grahl wrote: Dear All: I have a question regarding the behavior of functions. Say I define a function that returns another function : A - function(parameters) {      # calculations w/

Re: [R] Functions returning functions

2009-05-20 Thread Romain Francois
Romain Francois wrote: Paulo Grahl wrote: Dear All: I have a question regarding the behavior of functions. Say I define a function that returns another function : A - function(parameters) { # calculations w/ parameters returning 'y' tmpf - function(x) { # function of 'y' }

[R] turning off specific types of warnings

2009-05-20 Thread Eleni Rapsomaniki
Dear R users, I have a long function that among other things uses the survest function from the Design package. This function generates the warning: In survest.cph (...) S.E. and confidence intervals are approximate except at predictor means. Use cph(...,x=T,y=T) (and don't use

Re: [R] turning off specific types of warnings

2009-05-20 Thread jim holtman
?suppressWarnings On Wed, May 20, 2009 at 8:10 AM, Eleni Rapsomaniki er...@medschl.cam.ac.ukwrote: Dear R users, I have a long function that among other things uses the survest function from the Design package. This function generates the warning: In survest.cph (...) S.E. and confidence

Re: [R] turning off specific types of warnings

2009-05-20 Thread Romain Francois
jim holtman wrote: ?suppressWarnings Note that this removes __all__ warnings, if you only want to remove specific warnings, you can adapt suppressWarnings like this: h - function(w) if( any(grepl(ouch, w))) invokeRestart(muffleWarning) withCallingHandlers( warning(ouch), warning = h )

Re: [R] Functions returning functions

2009-05-20 Thread Paulo Grahl
Thanks a lot ! regards, Paulo Gustavo Grahl, CFA On Wed, May 20, 2009 at 8:31 AM, Romain Francois romain.franc...@dbmail.com wrote: Paulo Grahl wrote: Dear All: I have a question regarding the behavior of functions. Say I define a function that returns another function : A -

[R] how to merge replicate spots

2009-05-20 Thread Barbara Cegielska
Hi, Does anybody know how to merge replicate spots in R Bioconductor? I use home-made spotted arrays and every probe is printed in triplicate. I normalize my data using limma package and I do not know how to get mean from these three replicates. Regards, Barbara Cegielska

Re: [R] Extracting correlation in a nlme model

2009-05-20 Thread Christine A.
Hello, for calculation of intraclass correlation coefficient for a simple Random Intercept model, the function ICC1.lme() in the package psychometric can be applied. Best regards, Christine Adrion Kenneth Roy Cabrera Torres wrote: Hi R users: Is there a function to obtain the

Re: [R] Least-square support vector machines regression!

2009-05-20 Thread Thomas Terhoeven-Urselmans
Dear Max, thanks for your valuable comment. I assume that you used the function for regression - not classification. I use Mac OS X plattform (version 10.5.6). The R version is 2.8.1 (I prefer to update to 2.9.1 not 2.9.0). The kernlab package version is 0.9-8. The x and y-input into

[R] Plot data from table with column and row names

2009-05-20 Thread Matt Bishop
Dear All Sorry for what appears a trivial matter - I'm new to R and am stumbling ahead. I have a table of numerical data (36 rows by 12 columns) such as below: GM1 GM2 GM3 GM4 GM5 ...etc GM12 Run1 1 2 1 2 3 ... Run2 2 1 3 2 1 ... ...

[R] Mixed Models

2009-05-20 Thread Alan O'Loughlin
Hello, I've run a Proc Mixed function on a set of data in SAS. The data was a result of an experiment that measured % viability over time and I wanted to compare a Large sample lets say 50L to a small sample say 5L. And compare the % viability between the 2 sizes to see if I get the same

Re: [R] Least-square support vector machines regression!

2009-05-20 Thread Max Kuhn
Thomas, thanks for your valuable comment. I assume that you used the function for regression - not classification. I have been using it for classification and that is the issue. Looking at ?lssvm, it has regression is currently not supported in the details for the type argument. Max

Re: [R] Where to find a changelog for the survival package

2009-05-20 Thread Richard . Cotton
since some days I try to use the versions 2.35-4 of the survival package instead of versions 2.31, I had installed until now. Several changes in print.survfit, plot.survfit and seemingly in the structure of ratetabels effect some of my syntax files. Is there somewhere a documentation of

[R] efficiency when processing ordered data frames

2009-05-20 Thread Brigid Mooney
Hoping for a little insight into how to make sure I have R running as efficiently as possible. Suppose I have a data frame, A, with n rows and m columns, where col1 is a date time stamp. Also suppose that when this data is imported (from a csv or SQL), that the data is already sorted such that

[R] How to google for R stuff?

2009-05-20 Thread Kynn Jones
Hi! I'm new to R programming, though I've been programming in other languages for years. One thing I find most frustrating about R is how difficult it is to use Google (or any other search tool) to look for answers to my R-related questions. With languages with even slightly more distinctive

Re: [R] Plot data from table with column and row names

2009-05-20 Thread Gabor Grothendieck
Here are a few possibilities: Lines - GM1 GM2 GM3 GM4 GM5 Run1 1 2 1 2 3 Run2 2 1 3 2 1 Run3 2 1 1 1 1 DF - read.table(textConnection(Lines), header = TRUE) long - as.data.frame.table(as.matrix(DF)) head(long) with(long,

Re: [R] Functions returning functions

2009-05-20 Thread Stavros Macrakis
On Wed, May 20, 2009 at 7:21 AM, Paulo Grahl pgr...@gmail.com wrote: A - function(parameters) { # calculations w/ parameters returning 'y' tmpf - function(x) { # function of 'y' } return(tmpf) } The value of the parameters are stored in an environment local to the function.

Re: [R] Too large a data set to be handled by R?

2009-05-20 Thread Stavros Macrakis
On Tue, May 19, 2009 at 11:59 PM, tsunhin wong thjw...@gmail.com wrote: In order to save time, I am planning to generate a data set of size 1500 x 2 with each data point a 9-digit decimal number, in order to save my time. I know R is limited to 2^31-1 and that my data set is not going to

Re: [R] How to google for R stuff?

2009-05-20 Thread Jonathan Baron
See http://finzi.psych.upenn.edu/ for many options for searching for R stuff. If you type R into Google, the R home page usually comes up first. I was amazed when this happened 9 years ago. And now, if I type R lmer into Google (without the quotes), a whole bunch of relevant stuff comes up,

[R] Error with regsubset in leaps package - vcov and all.best option (plus calculating VIFs for subsets)

2009-05-20 Thread lara harrup (IAH-P)
Hi all I am hoping this is just a minor problem, I am trying to implement a best subsets regression procedure on some ecological datasets using the regsubsets function in the leaps package. The dataset contains 43 predictor variables plus the response (logcount) all in a dataframe called

Re: [R] How to google for R stuff?

2009-05-20 Thread Ian Fiske
www.rseek.org is the best solution to this that I have found. Ian kynn wrote: Hi! I'm new to R programming, though I've been programming in other languages for years. One thing I find most frustrating about R is how difficult it is to use Google (or any other search tool) to look for

Re: [R] efficiency when processing ordered data frames

2009-05-20 Thread jim holtman
How much is it currently costing you in time to do the selection process? Is it having a large impact on your program? Is it the part that is really consuming the overall time? What is your concern in this area? Here is the timing that it take so select from 10M values those that are less than a

Re: [R] How to google for R stuff?

2009-05-20 Thread Richard . Cotton
One thing I find most frustrating about R is how difficult it is to use Google (or any other search tool) to look for answers to my R-related questions. With languages with even slightly more distinctive names like Perl, Java, Python, Matlab, OCaml, etc., usually including the name of the

Re: [R] How to google for R stuff?

2009-05-20 Thread Luc Villandre
Kynn Jones wrote: Hi! I'm new to R programming, though I've been programming in other languages for years. One thing I find most frustrating about R is how difficult it is to use Google (or any other search tool) to look for answers to my R-related questions. With languages with even

Re: [R] SEM:Standard error of std.coef estimates?

2009-05-20 Thread John Fox
Dear Bastiaan, The standard errors of the standardized coefficients aren't simple because the standard deviations used to standardize the coefficients are also subject to sampling error. I can think of two ways to get standard errors for the standardized coefficients: by the delta method and by

[R] Trouble installing package 'systemfit'

2009-05-20 Thread Rui Wang
I tried unnder R 2.9.0 and R 2.8.1. Both versions failed to install 'systemfit'. I downloaded file 'systemfit_1[1].0-8.zip' to install the package. I am using Windows 2000 service pack 4. The error message is: Error in gzfile(file, r) : cannot open the connection In addition: Warning

Re: [R] Getting lm() to work with a matrix

2009-05-20 Thread Luc Villandre
MikSmith wrote: Hi I'm fairly new to R and am trying to analyse some large spectral datasets using stepwise regression (fairly standard in this area). I have a field sampled dataset, of which a proportion has been held back for validation. I gather than step() needs to be fed a regression model

Re: [R] How to google for R stuff?

2009-05-20 Thread cryan
For Google searches, I find that throwing in the term cran on every search helps weed out irrelevant pages. For example, instead of r residuals I type r cran residuals --Chris Ryan Original message Date: Wed, 20 May 2009 09:43:14 -0400 From: Luc Villandre

Re: [R] Trouble installing package 'systemfit'

2009-05-20 Thread Ronggui Huang
Have you tried another mirror? It seems that downloaded file is incomplete. 2009/5/20 Rui Wang r...@hawaii.edu: I tried unnder R 2.9.0 and R 2.8.1. Both versions failed to install 'systemfit'. I downloaded file 'systemfit_1[1].0-8.zip' to install the package. I am using Windows 2000 service

Re: [R] Trouble installing package 'systemfit'

2009-05-20 Thread Uwe Ligges
Rui Wang wrote: I tried unnder R 2.9.0 and R 2.8.1. Both versions failed to install 'systemfit'. I downloaded file 'systemfit_1[1].0-8.zip' to install the package. I am using Windows 2000 service pack 4. The filename lost its canonical form, pelase try to rename it to: systemfit_1.0-8.zip

Re: [R] Error with regsubset in leaps package - vcov and all.best option (plus calculating VIFs for subsets)

2009-05-20 Thread Thomas Lumley
On Wed, 20 May 2009, lara harrup (IAH-P) wrote: Hi all I am hoping this is just a minor problem, I am trying to implement a best subsets regression procedure on some ecological datasets using the regsubsets function in the leaps package. The dataset contains 43 predictor variables plus

[R] Random effects aov

2009-05-20 Thread Christine Griffiths
Dear All I have a repeated measures design in which abundance was measured repeatedly over 10 months in three treatments (Tortoise A; Tortoise B and control) established in 6 blocks, i.e. crossed fixed effects. My original design incorporated two tortoises per treatment, however as fieldwork

[R] function returning the name of an object

2009-05-20 Thread Antonio Gasparrini
Hello R users, I'm trying to solve an apparently simple problem. I include a matrix in a model formula, and I want then to extract the related coefficients. Basically, I want a function func such that: mymatrix - cbind(rnorm(10),rnorm(10)) y - b - rnorm(10) model - lm(y ~ b + mymatrix)

Re: [R] Getting lm() to work with a matrix

2009-05-20 Thread Berwin A Turlach
G'day Luc, On Wed, 20 May 2009 09:58:41 -0400 Luc Villandre villa...@dms.umontreal.ca wrote: MikSmith wrote: [...] Indeed, functions like /lm()/ require the object fed to the /data/ argument to be either [...] But the data argument is optional and does not need to be specified. In your

Re: [R] Trouble installing package 'systemfit'

2009-05-20 Thread Rui Wang
Renaming the file worked! Thanks! - Original Message - From: Uwe Ligges lig...@statistik.tu-dortmund.de Date: Wednesday, May 20, 2009 10:07 am Subject: Re: [R] Trouble installing package 'systemfit' To: Rui Wang r...@hawaii.edu Cc: r-help@r-project.org Rui Wang wrote: I tried

[R] problem with APPLY

2009-05-20 Thread De France Henri
Hello,   The apply function seems to behave oddly with my code below   NB : H1 is a data frame. (data in the attached file.) # the first lines are: 1 02/01/2008 0.00  0  0 0.00   0 2 03/01/2008 0.00  0  0 0.00   0 3 04/01/2008 0.00  0  0 0.00   0 4 07/01/2008 0.00  0  0

Re: [R] How to google for R stuff?

2009-05-20 Thread spencerg
There is also the relatively new RSiteSearch package. It's RSiteSearch.function searches only help pages of contributed packages but returns the result in a data.frame (of class RSiteSearch) sorted to put the most interesting package first with help pages sorted within packages. If this

Re: [R] SEM:Standard error of std.coef estimates?

2009-05-20 Thread BdeGroot
Dear John, Thank you for your reply. I would like to try the bootstrapping method. I have the boot package, however, I have not used it before. I do not know how to compute the correct code to calculate the standarized errors. If possible and easily achievable, could you please tell me what to

[R] Class for time of day?

2009-05-20 Thread Stavros Macrakis
What is the recommended class for time of day (independent of calendar date)? And what is the recommended way to get the time of day from a POSIXct object? (Not a string representation, but a computable representation.) I have looked in the man page for DateTimeClasses, in the Time Series

[R] meaning of order-vector in object of class hclust

2009-05-20 Thread Pieter van der Spek
I'm trying to export the dendrogram created in R such that I can rebuild it inside an other program. However, I don't fully understand the relation between the merge, height and order vectors. This is what I do grasp. The merge-matrix shows how the original observations were merged into a tree

[R] stationarity tests

2009-05-20 Thread mauede
How can I make sure the residual signal, after subtracting the trend extracted through some technique, is actually trend-free ? I would greatly appreciate any suggestion about some Stationarity tests. I'd like to make sure I have got the difference between ACF and PACF right. In the following

Re: [R] How to google for R stuff?

2009-05-20 Thread Seeliger . Curt
Kynn writes: So I'm curious to learn what strategies R users have found to get around this annoyance. I use Rseek for most of my R questions: http://www.rseek.org/ cur -- Curt Seeliger, Data Ranger Raytheon Information Services - Contractor to ORD seeliger.c...@epa.gov 541/754-4638

Re: [R] arrangement of crowded labels

2009-05-20 Thread Greg Snow
Look at the spread.labs and the dynIdentify and TkIdentify functions in the TeachingDemos package. -- Gregory (Greg) L. Snow Ph.D. Statistical Data Center Intermountain Healthcare greg.s...@imail.org 801.408.8111 -Original Message- From: r-help-boun...@r-project.org

Re: [R] drc results differ for different versions

2009-05-20 Thread Marc Schwartz
On May 20, 2009, at 11:20 AM, Hans Vermeiren wrote: Hello, We use drc to fit dose-response curves, recently we discovered that there are quite different standard error values returned for the same dataset depending on the drc-version / R-version that was used (not clear which factor is

Re: [R] How to google for R stuff?

2009-05-20 Thread spencerg
Dear Dieter: Thanks for the correction. I failed to test the code as written before I posted it. Spencer Graves Dieter Menne wrote: spencerg spencer.graves at prodsyse.com writes: de - RSiteSearch.function(differential equation) des - RSiteSearch.function(differential

[R] evaluating partitions in CLUE

2009-05-20 Thread A M Lavezzi
Hello, I'm starting to use CLUE to compare clusterings. I have some cluster structures obtained outside R, which are simlply vectors of lenght N (N: number of objects), which contain a number indicating the cluster to which each object belongs. I report below a portion of one of these:

[R] drc results differ for different versions

2009-05-20 Thread Hans Vermeiren
Hello, We use drc to fit dose-response curves, recently we discovered that there are quite different standard error values returned for the same dataset depending on the drc-version / R-version that was used (not clear which factor is important) On R 2.9.0 using drc_1.6-3 we get an IC50 of

Re: [R] RGoogleDocs: can now see documents but cannot get content.

2009-05-20 Thread Duncan Temple Lang
Hi Farrel This particular problem is a trivial issue of an argument out of place due to a change in the function definition during the development. There is a new version of the package (0.2-2) and it also uses a slightly different approach (and function) to pull the values into the form of

Re: [R] SEM:Standard error of std.coef estimates?

2009-05-20 Thread John Fox
Dear Bastiaan, I've written an appendix on bootstrapping regression models, available at http://socserv.mcmaster.ca/jfox/Books/Companion/appendix-bootstrapping.pdf , which describes generally how to proceed. In outline, you'll write a function that takes your data matrix (not covariance matrix)

[R] How to load data from Statistics Canada

2009-05-20 Thread guox
We would like to load data from Statistics Canada (http://www.statcan.gc.ca/) using R, for example, Employment and unemployment rates. It seems to me that the tables are displayed in HTML. I was wondering if you know how to load these tables. Thanks, -james

[R] Comparing spatial distributions - permutation test implementation

2009-05-20 Thread JiHO
Hello everyone, I am looking at the joint spatial distribution of 2 kinds of organisms (estimated on a grid of points) and want to test for significant association or dissociation. My first question is: do you know a nice technique to do that, considering that I have a limited number of

Re: [R] problem with APPLY

2009-05-20 Thread Peter Dalgaard
De France Henri wrote: Hello, The apply function seems to behave oddly with my code below NB : H1 is a data frame. (data in the attached file.) # the first lines are: 1 02/01/2008 0.00 0 0 0.00 0 2 03/01/2008 0.00 0 0 0.00 0 3 04/01/2008 0.00 0 0 0.00 0 4

[R] Supressing the enumeration of output in console

2009-05-20 Thread ZeMajik
Hi! Pretty low content question but I've had major trouble finding an answer for it so I hope it's alright. I'm obviously new to R, and have been trying to get rid of the numerated output I get in the console. What I mean more specifically is that X-4;X comes out as [1] 4 and I'd like to get

Re: [R] Class for time of day?

2009-05-20 Thread Gabor Grothendieck
There is a times class in the chron package. Times are measured in fractions of a day so 1/24 is one hour. library(chron) dt - Sys.time() tt - times(format(dt, %H:%M:%S)) tt [1] 12:27:46 tt + 1/24 [1] 13:27:46 There is an article on dates and times in R News 4/1. On Wed, May 20, 2009 at

Re: [R] sem with categorical data

2009-05-20 Thread John Fox
Dear Frank, Did you take a look at the data you generated? The contingency tables are extremely sparse. Try, e.g., with(dat, table(q1, q2)). (I assume, BTW, that you're using the cut2() function in the Hmisc package.) Beyond that, I'm not sure what the point is of your example. You're generating

Re: [R] Supressing the enumeration of output in console

2009-05-20 Thread Greg Snow
Use the 'cat' function to put text on the console without any extra decoration (you are then responsible for all formatting including line feeds). Hope this helps, -- Gregory (Greg) L. Snow Ph.D. Statistical Data Center Intermountain Healthcare greg.s...@imail.org 801.408.8111 -Original

Re: [R] How to google for R stuff?

2009-05-20 Thread Joshua Stults
I was having the same frustration so I made a Google custom search engine for Open Math Tools: http://www.google.com/coop/cse?cx=015659631714008342946:wolsniqtrxc It gives preferences to results from r-project.org, octave.org, maxima.sourceforge.net, and the mailing list archives from those

Re: [R] Class for time of day?

2009-05-20 Thread jim holtman
If you want the hours from a POSIXct, here is one way of doing it; you can create a function for doing it: x - Sys.time() x [1] 2009-05-20 12:17:13 EDT y - difftime(x, trunc(x, units='days'), units='hours') y Time difference of 12.28697 hours as.numeric(y) [1] 12.28697 It depends on what

Re: [R] How to google for R stuff?

2009-05-20 Thread Dieter Menne
spencerg spencer.graves at prodsyse.com writes: de - RSiteSearch.function(differential equation) des - RSiteSearch.function(differential equations) # With de and des, each finds things missed by the other. de. - de | des # combine into one sumDE - packageSum2(de.) # add details on installed

Re: [R] How to load data from Statistics Canada

2009-05-20 Thread Luc Villandre
g...@ucalgary.ca wrote: We would like to load data from Statistics Canada (http://www.statcan.gc.ca/) using R, for example, Employment and unemployment rates. It seems to me that the tables are displayed in HTML. I was wondering if you know how to load these tables. Thanks, -james

[R] dropping empty panels in xyplot

2009-05-20 Thread Steve_Friedman
I note that someone asked for help a few weeks ago regarding the ability to drop empty panels on an xyplot. I did not see a reply to that inquiry so I am asking for assistance with the same problem To prepare the data for the plotting routine, I did the following to restrict the content of the

Re: [R] problem with APPLY

2009-05-20 Thread William Dunlap
-Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of Peter Dalgaard Sent: Wednesday, May 20, 2009 8:16 AM To: De France Henri Cc: r-help@r-project.org Subject: Re: [R] problem with APPLY De France Henri wrote: Hello,

Re: [R] Where to find a changelog for the survival package

2009-05-20 Thread Heinz Tuechler
Thank you Richie. I had seen this before, but my impression is that it's not up to date. I gave a wrong version number in my previous post. I changed from 2.34-1 to 2.35-4. For example, the plot.survfit function lost it's legend parameters, but I don't see this in the changelog. Thanks again,

Re: [R] How to load data from Statistics Canada

2009-05-20 Thread guox
Yes, XML can help load HTML tables. As we knew, tables in Bureau of Labor Statistics (http://www.bls.gov/) can be displayed in HTML and can also be loaded directly using their urls, for example, ftp://ftp.bls.gov/pub/time.series/wp/ I was wondering if tables in Statist Canada have such locations

Re: [R] How to load data from Statistics Canada

2009-05-20 Thread Peter Dalgaard
g...@ucalgary.ca wrote: We would like to load data from Statistics Canada (http://www.statcan.gc.ca/) using R, for example, Employment and unemployment rates. It seems to me that the tables are displayed in HTML. I was wondering if you know how to load these tables. Thanks, I suspect the

[R] Linear constraints for constrasts

2009-05-20 Thread Serguei Kaniovski
Dear List! How can I define contrasts (design matrix) that can all be included, i.e. which do not require a control category be dropped. My application (see below) does not suggest a sensible control category. I am thinking of constraining the (treatment) contrasts to sum up to zero and

[R] minmun p-value for wilcox.test and correlation

2009-05-20 Thread charles78
I have a stupid question on how to get the real p-values for wilcox.test and correlation. the minmun can be reached is 2.2E-16 using the R version 2.6.2. I do not think it is the R version causing this but other issues. Any help is highly appreciated. Charles -- View this message in context:

Re: [R] Predicting complicated GAMMs on response scale

2009-05-20 Thread William Paterson
Creation of Animal category in p.d solved all problems. Plots fine now. The smallest hurdles are often the hardest to get over. Gavin Simpson wrote: On Mon, 2009-05-18 at 11:48 -0700, William Paterson wrote: Hi, I am using GAMMs to show a relationship of temperature differential over

[R] sem with categorical data

2009-05-20 Thread Frank Lawrence
I am trying to run a confirmatory factor analysis using the SEM package. My data are ordinal. I have read http://socserv.mcmaster.ca/jfox/Misc/sem/SEM-paper.pdf. When I apply the hetcor function, I receive the following error: Error in checkmvArgs(lower = lower, upper = upper, mean = mean, corr

[R] error message re: max(i), but code and output seen O.K.

2009-05-20 Thread Kirsten Miles
I have a researcher who is consistently get the warning message: In max(i) : no non-missing arguments to max; returning -Inf Best as I can tell the code is working properly and the output is as expected. I would like some help in understanding why he is getting this error message and what its

[R] Example for parsing XML file?

2009-05-20 Thread Brigid Mooney
Hi, I am trying to parse XML files and read them into R as a data frame, but have been unable to find examples which I could apply successfully. I'm afraid I don't know much about XML, which makes this all the more difficult. If someone could point me in the right direction to a resource

[R] Non-linear regression with latent variable

2009-05-20 Thread Samiul Hasan
Hi Can anyone please suggest me a package where I can estimate a non-linear regression model? One of the independent variables is latent or unobserved. I have an indicator variable for this unobserved variable; however the relationship is known to be non-linear also. In terms of equations my

Re: [R] minmun p-value for wilcox.test and correlation

2009-05-20 Thread Peter Flom
charles78 hu...@hotmail.com wrote I have a stupid question on how to get the real p-values for wilcox.test and correlation. the minmun can be reached is 2.2E-16 using the R version 2.6.2. I do not think it is the R version causing this but other issues. Any help is highly appreciated. Can I

[R] round function seems to produce maximum 2 decimals

2009-05-20 Thread Glenn E Stauffer
I am trying to use round()to force R to display a specific number of decimals, but it seems to display =2 decimals no matter what I specify in the digits argument. As an alternative I tried signif(), but it also produces unexpected results. See example code and results below. Format() works, but

Re: [R] How to load data from Statistics Canada

2009-05-20 Thread Wacek Kusnierczyk
Peter Dalgaard wrote: g...@ucalgary.ca wrote: We would like to load data from Statistics Canada (http://www.statcan.gc.ca/) using R, for example, Employment and unemployment rates. It seems to me that the tables are displayed in HTML. I was wondering if you know how to load these tables.

Re: [R] error message re: max(i), but code and output seen O.K.

2009-05-20 Thread Sundar Dorai-Raj
This error is thrown if the argument to max is either NULL or length zero: [~] Rscript -e max(NULL) [1] -Inf Warning message: In max(NULL) : no non-missing arguments to max; returning -Inf [~] Rscript -e max(numeric(0)) [1] -Inf Warning message: In max(numeric(0)) : no non-missing arguments to

Re: [R] minmun p-value for wilcox.test and correlation

2009-05-20 Thread Peter Dalgaard
charles78 wrote: I have a stupid question on how to get the real p-values for wilcox.test and correlation. the minmun can be reached is 2.2E-16 using the R version 2.6.2. I do not think it is the R version causing this but other issues. Any help is highly appreciated. Charles It's mainly a

Re: [R] error message re: max(i), but code and output seen O.K.

2009-05-20 Thread Rolf Turner
On 21/05/2009, at 6:23 AM, Kirsten Miles wrote: I have a researcher who is consistently get the warning message: In max(i) : no non-missing arguments to max; returning -Inf Best as I can tell the code is working properly and the output is as expected. I would like some help in understanding

Re: [R] How to google for R stuff?

2009-05-20 Thread Kynn Jones
Thank you all very much for the so many useful ideas and resources. KJ [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide

[R] how to get remote ESS graphics to work?

2009-05-20 Thread Matthew Keller
Hi all, My graduate student is logging onto my macpro and running R through ESS aquamacs (with Mx ssh and then Mx ess-remote). Everything is working fine until we get to graphing. We are trying to give him the ability to look at graphics interactively. The ESS manual is not too helpful: If you

Re: [R] Example for parsing XML file?

2009-05-20 Thread Duncan Temple Lang
Hi Brigid. Here are a few commands that should do what you want: bri = xmlParse(myDataFile.xml) tmp = t(xmlSApply(xmlRoot(bri), xmlAttrs))[, -1] dd = as.data.frame(tmp, stringsAsFactors = FALSE, row.names = 1:nrow(tmp)) And then you can convert the columns to whatever

Re: [R] Example for parsing XML file?

2009-05-20 Thread Wacek Kusnierczyk
Brigid Mooney wrote: Hi, I am trying to parse XML files and read them into R as a data frame, but have been unable to find examples which I could apply successfully. I'm afraid I don't know much about XML, which makes this all the more difficult. If someone could point me in the right

Re: [R] How to google for R stuff?

2009-05-20 Thread Emmanuel Charpentier
Le mercredi 20 mai 2009 à 09:02 -0400, Kynn Jones a écrit : Hi! I'm new to R programming, though I've been programming in other languages for years. One thing I find most frustrating about R is how difficult it is to use Google (or any other search tool) to look for answers to my R-related

[R] combining xYplot with map

2009-05-20 Thread eric.archer
I'm using xYplot to create a bubble plot of values that I'd like to visualize on top of a filled-in map of the coast, but I'm too much of a lattice (which I understand xYplot is built on) and mapping newbie to figure out how to begin to make this happen. Below is some sample code that doesn't

Re: [R] round function seems to produce maximum 2 decimals

2009-05-20 Thread Ted Harding
On 20-May-09 20:10:15, Glenn E Stauffer wrote: I am trying to use round()to force R to display a specific number of decimals, but it seems to display =2 decimals no matter what I specify in the digits argument. As an alternative I tried signif(), but it also produces unexpected results. See

Re: [R] round function seems to produce maximum 2 decimals

2009-05-20 Thread Ben Bolker
Glenn E Stauffer wrote: I am trying to use round()to force R to display a specific number of decimals, but it seems to display =2 decimals no matter what I specify in the digits argument. As an alternative I tried signif(), but it also produces unexpected results. See example code and

Re: [R] error message re: max(i), but code and output seen O.K.

2009-05-20 Thread David Scott
On Wed, 20 May 2009, Kirsten Miles wrote: I have a researcher who is consistently get the warning message: In max(i) : no non-missing arguments to max; returning -Inf You get this message when trying to take the maximum of a vector with no elements. I am guessing he is looping through and

Re: [R] combining xYplot with map

2009-05-20 Thread eric.archer
I'll retract my request for help. I managed to solve the problem by calling map followed by a call to points after rescaling the cex. Its the kind of solution that deserves a head-slap. library(Hmisc) library(maps) sample.data - data.frame(lat = c(12.1667, 14.6333, -6.874, 2.6167, 14.6833,

  1   2   >