Re: [R] Execute script and go interactive

2009-02-28 Thread Prof Brian Ripley
On Fri, 27 Feb 2009, Michael Olschimke wrote: Hi, is it possible to execute a script with R and go into interactive mode with the same session? Yes, the mechanism is 'read the script file specified by R_PROFILE_USER (or its defaults) and then continue interactively'. See ?Startup. e.g.

[R] Rprofiling

2009-02-28 Thread Erin Hodgess
Dear R People: Could someone recommend a good reference for Rprofing, please? Thanks, Erin -- Erin Hodgess Associate Professor Department of Computer and Mathematical Sciences University of Houston - Downtown mailto: erinm.hodg...@gmail.com __

Re: [R] Anova Help?

2009-02-28 Thread Tal Galili
Hi Dar. It sounds like anova is a good way to go. The command to look into in R is: ?aov you might also go into manova (for anova for two dependent variables together) But make sure you understand the background of it when you work the commands. Cheers, Tal On Wed, Feb 25, 2009 at 11:06

Re: [R] Advice on graphics to design circle with density-shaded sectors

2009-02-28 Thread Jim Lemon
John Poulsen wrote: Hello, I am looking for some general advice on which graphics package to use to make a figure demonstrating my experimental design. I want to design a circle with 7 sectors inside. Then I will want to shade the sectors depending on densities of observations in the

[R] Using JRI and Java 1.6 on MAC OS X

2009-02-28 Thread dobomode
Dear R-Help, I am trying to get JRI (the rJava interface allowing Java to connect to R) to work. I was able to run it a week ago when I was doing some testing using Java 1.5. However, I am developing a GUI application using some of the new Java 1.6 features and I just can't get JRI to work with

Re: [R] Rprofiling

2009-02-28 Thread Gabor Grothendieck
If f is your function just do this: Rprof(); f(); Rprof(NULL); summaryProf() Also check out the profr package and for benchmarking check out ?system.time and http://rbenchmark.googlecode.com On Sat, Feb 28, 2009 at 4:15 AM, Erin Hodgess erinm.hodg...@gmail.com wrote: Dear R People: Could

Re: [R] Using JRI and Java 1.6 on MAC OS X

2009-02-28 Thread Prof Brian Ripley
There is no fundamental problem with rJava/JRI and Java 1.6: Java 1.6 is all that will work on 64-bit Linux and AMD Solaris platforms, so some of us have only used Java 1.6 for years. This is really an issue for the R-sig-mac list and/or the rJava maintainer (but he is travelling at present).

[R] Review my upgrade plan from 2.8.0 to 2.8.1

2009-02-28 Thread Daren Tan
Hi, On my laptop, R is installed on windows XP SP2 at D:\Program Files\R\R-2.8.0, and all add-on packages are installed at D:\Program Files\R\R-2.8.0.libs. In addition, I have created two environment enviroment to ease upgrading and installation of packages. Packages installed is a mix of those

[R] Implementation of quasi-bayesian maximum likelihood estimation for normal mixtures

2009-02-28 Thread Helena Richter
Hi, as you can see in the topic, I am trying to fit a normal mixture distribution with the approach suggested by Hamilton (1991). Since I couldn't find any existing packages including the quasi-bayesian mle, I have to write my own function. Unfortunately, I have absolutely no experience in

Re: [R] Large 3d array manipulation

2009-02-28 Thread David Winsemius
On Feb 27, 2009, at 10:37 PM, David Winsemius wrote: On Feb 27, 2009, at 7:49 PM, Duncan Murdoch wrote: On 27/02/2009 6:15 PM, Vemuri, Aparna wrote: I have a large 3 dimensional array of size (243,246,768) The first dimension is Rows, second is columns and the third is Time. So for each

Re: [R] Review my upgrade plan from 2.8.0 to 2.8.1

2009-02-28 Thread Gabor Grothendieck
I am on Vista so some of this might be slightly different on XP but I only keep one 2.8.x version, one 2.7.x version, etc. For me, all 2.8.x versions of R go into: %ProgramFiles%\R\R-2.8.x so installing R-2.8.1 is just a matter of installing it over R-2.8.0 and there is nothing more to do. In

[R] lme4 and Variable level detection

2009-02-28 Thread Jeroen Ooms
I am making a little GUI for lme4, and I was wondering if there is a function that automatically detects on which level every variable exists. Furtheremore I got kind of confused about what a random effects model actually calculates. I have some experience with commercial software packages for

Re: [R] Rprofiling

2009-02-28 Thread Dirk Eddelbuettel
Erin, On 28 February 2009 at 03:15, Erin Hodgess wrote: | Dear R People: | | Could someone recommend a good reference for Rprofing, please? Well, you have the notes of the 'high performance computing with R' tutorial at last summer's useR conference giving you examples using both 'base R'

[R] Hotelling-Pabst-Statistic

2009-02-28 Thread Alexandros Droseltis
Hello all! [I am using R version 2.8.0 (2008-10-20).] Is there a package that provides the critical values for the Hotelling-Pabst-Statistic? (Useful for a Spearman Test). Thanks for any hint. Alexandros __ R-help@r-project.org mailing list

[R] problem with an old .RData file

2009-02-28 Thread Marco Barbàra
Dear useRs, Today I've encountered a weird unexpected behaviour. It is not a real problem for me and I can workaround it. But I'm very curious of what could have caused it. Starting an R session in a directory with an old .RData my first command was library(MASS) and i got Errore in

Re: [R] Competing risks adjusted for covariates

2009-02-28 Thread Arthur Allignol
You could try the prodlim package. Best regards, Arthur Allignol On Fri, 27 Feb 2009 19:36:31 - Eleni Rapsomaniki er...@medschl.cam.ac.uk wrote: Dear R-users Has anybody implemented a function/package that will compute an individual's risk of an event in the presence of competing

Re: [R] Adjusting confidence intervals for paired t-tests o f multiple endpoints

2009-02-28 Thread Dieter Menne
Erich Studerus erich.studerus at bli.uzh.ch writes: Dear R-users, In a randomized placebo-controlled within-subject design, subjects recieved a psycho-active drug and placebo. Subjects filled out a questionnaire containing 15 scales on four different time points after drug

Re: [R] problem with an old .RData file

2009-02-28 Thread Prof Brian Ripley
When you load your old workspace, try conflicts(detail = TRUE) The workspace contains a lazyload promise from that very old version, with the name of a system function. You need to delete that promise from the workspace, and the command above should help you find which it is. One guess is

Re: [R] lme4 and Variable level detection

2009-02-28 Thread Douglas Bates
On Sat, Feb 28, 2009 at 9:00 AM, Jeroen Ooms j.c.l.o...@uu.nl wrote: I am making a little GUI for lme4, and I was wondering if there is a function that automatically detects on which level every variable exists. Furtheremore I got kind of confused about what a random effects model actually

[R] 3d graph aspect ratio question‏

2009-02-28 Thread Bo Zhou
3d graph aspect ratio questionþ Hi guys, What's the counterpart of rgl package's aspect3d() in Lattice package? aspect3d() :

[R] Extract statistics from lm()

2009-02-28 Thread Bogaso
Hi, perhaps this question was answered previously however I could not find them. My problem is how how to extract a particular statistic from the result given by lm(). For e.g. ctl - c(4.17,5.58,5.18,6.11,4.50,4.61,5.17,4.53,5.33,5.14) trt - c(4.81,4.17,4.41,3.59,5.87,3.83,6.03,4.89,4.32,4.69)

Re: [R] Extract statistics from lm()

2009-02-28 Thread David Winsemius
http://finzi.psych.upenn.edu/R/Rhelp08/archive/147713.html On Feb 28, 2009, at 2:21 PM, Bogaso wrote: Hi, perhaps this question was answered previously however I could not find them. My problem is how how to extract a particular statistic from the result given by lm(). For e.g. ctl -

Re: [R] Review my upgrade plan from 2.8.0 to 2.8.1

2009-02-28 Thread Uwe Ligges
Daren Tan wrote: Hi, On my laptop, R is installed on windows XP SP2 at D:\Program Files\R\R-2.8.0, and all add-on packages are installed at D:\Program Files\R\R-2.8.0.libs. In addition, I have created two environment enviroment to ease upgrading and installation of packages. Packages

Re: [R] Extract statistics from lm()

2009-02-28 Thread Tobias Verbeke
Hi, Hi, perhaps this question was answered previously however I could not find them. My problem is how how to extract a particular statistic from the result given by lm(). For e.g. ctl - c(4.17,5.58,5.18,6.11,4.50,4.61,5.17,4.53,5.33,5.14) trt -

[R] filled.contour and time axis

2009-02-28 Thread Pieter Beck
Dear all, I am using filled.contour and my x-axis consists of a vector of dates, formatted using as.POSIXct. When I assign this vector to x in filled.contour, the x-axis in the plot is correctly formatted as dates. However, when I use plot.axes=axis(1) in filled.contour the x-axis is not

Re: [R] Extract statistics from lm()

2009-02-28 Thread Albyn Jones
Look at the data structure produced by summary() names(summary(lm.D9)) [1] call terms residuals coefficients [5] aliased sigma dfr.squared [9] adj.r.squared fstatisticcov.unscaled Now look at the data structure for the coefficients in the

Re: [R] object .trPaths not found

2009-02-28 Thread Uwe Ligges
rkevinbur...@charter.net wrote: I am running an R script with Tinn-R (2.2.0.1) and I get the error message Error in source(.trPaths[4], echo = TRUE, max.deparse.length = 150) : object .trPaths not found Any solutions? Maybe, but the may depend on your script, your OS, your R version,

Re: [R] filled.contour and time axis

2009-02-28 Thread David Winsemius
The filled.contour help page offers a somewhat more extensive example for using plot.axes than you have attempted: plot.axes = { axis(1, seq(100, 800, by = 100)) axis(2, seq(100, 600, by = 100)) }, And looking at the help page for axis() shows that the default setting for at is NULL, so it

Re: [R] Adjusting confidence intervals for paired t-tests of multiple endpoints

2009-02-28 Thread Erich Studerus, Psychiatrische Uni-Klinik
Thanks, but as far as I know the multcomp package and the tukeyHSD function can only handle one dependent variable, whereas I have 60 dependent variables. Regards, Erich On Sat, 28 Feb 2009 17:41:11 + (UTC) Dieter Menne dieter.me...@menne-biomed.de wrote: Erich Studerus erich.studerus

Re: [R] Statistic community?

2009-02-28 Thread John Kane
Usenet groups sci.stat.math, sci.stat.consult ? --- On Wed, 2/25/09, Peters Gj (PSYCHOLOGY) gj.pet...@psychology.unimaas.nl wrote: From: Peters Gj (PSYCHOLOGY) gj.pet...@psychology.unimaas.nl Subject: Re: [R] Statistic community? To: r-h...@stat.math.ethz.ch Received: Wednesday, February

[R] rgl persp3d bounding box color problem.

2009-02-28 Thread Bo Zhou
Hi guys, I hit on a problem when I use rgl. Could you try to run the code here in this link and see why the first persp3d gives a red bounding box and the second shows black? http://rafb.net/p/g1i7ur33.html (sorry for not pasting the code here directly but my previous email got filtered by

Re: [R] filled.contour and time axis

2009-02-28 Thread David Winsemius
Pieter reports that a further suggestion of using, instead, axis.POSIXct(1,x) was helpful. -- David Winsemius On Feb 28, 2009, at 2:53 PM, David Winsemius wrote: The filled.contour help page offers a somewhat more extensive example for using plot.axes than you have attempted: plot.axes

[R] creating log-log survival plots that are not inverted

2009-02-28 Thread Bob Green
I am hoping for some advice regarding how to obtain a log-log survival plot that is not in the inverse. On page 128 of Modelling survival data by Therneau Grambsch there is the an example of the type of desired plot, with a log of the survival curve by years. Marc Schwartz has provided me

Re: [R] object .trPaths not found

2009-02-28 Thread Kingsford Jones
RSiteSearch('trpaths') will lead you to potential solutions... On Fri, Feb 27, 2009 at 7:18 PM, rkevinbur...@charter.net wrote: I am running an R script with Tinn-R (2.2.0.1) and I get the error message Error in source(.trPaths[4], echo = TRUE, max.deparse.length = 150) :  object .trPaths

[R] How to join many records against SQL database

2009-02-28 Thread Andrew Ziem
This is a working example of how to merge records with a SQL database given the constraints 1. The database is too large to pull all the records 2. The database permissions don't allow creating a table for temporarily storing identifiers 3. The R database driver doesn't allow creating temporary

Re: [R] object .trPaths not found

2009-02-28 Thread rkevinburton
FYI. I found the solution. My RProfile.site file could not be written to because of permissions. Whan I selected the Configure/Permanent option in Tinn-R it was silently ignoring the fact that the file could not be written to. When I adjusted the permissions, all was well. Thank you. Kevin

[R] How to create temporary table in MySQL

2009-02-28 Thread Andrew Ziem
Creating a temp table isn't completely intuitive with MySQL 5 and R 2.8..1, but it can be done. library(RMySQL) Loading required package: DBI con - dbConnect(dbDriver(MySQL), dbname = foo, user=me,password=secret) x- data.frame(1:10) colnames(x) -c(x) dbWriteTable(con, #x, x,

[R] Temporary tables with Microsoft SQL?

2009-02-28 Thread Andrew Ziem
I can create a temp table with MySQL and R DBI[1], but I don't see how to do the same with Microsoft SQL 2005 and RODBC. R 2.8.1 creates the table, but then it can never see it. I'm looking to avoid replacing the convenience functions like sqlSave(). [1]

[R] Variable scope.

2009-02-28 Thread rkevinburton
I have a question on scope/reference/value type of variables with 'R'. The issue cam up first when I look at the arima code. I see code like: myupARIMA - function(mod, phi, theta) { . . . . mod } Then armafn - function(p, trans) { . . . . Z -

Re: [R] creating log-log survival plots that are not inverted

2009-02-28 Thread David Winsemius
I think what you want may be produced by this code for InvNormal(S) vs log(time): survplot(fit, fun=qnorm, logT=T, conf = none) That is not what you describe, however. I am worried about the plot on the page you cite, because it is not similar to other log(-log(S)) (complementary log-log)

Re: [R] Adjusting confidence intervals for paired t-tests of multiple endpoints

2009-02-28 Thread Johannes Huesing
Erich Studerus, Psychiatrische Uni-Klinik erich.stude...@bli.uzh.ch [Sat, Feb 28, 2009 at 09:23:05PM CET]: Thanks, but as far as I know the multcomp package and the tukeyHSD function can only handle one dependent variable, whereas I have 60 dependent variables. Couldn't you re-phrase your

Re: [R] Temporary tables with Microsoft SQL?

2009-02-28 Thread Prof Brian Ripley
This is an issue for the R-sig-db list (as was your other posting). RODBC is written to use standard SQL, and Microsft often does not use standard SQL, as in table names like #x. The most productive thing to do would be to modify RODBC to handle the Microsoft peculiarities you want and feed