Re: [R] R-2.6.0 and RWinEdt

2007-10-07 Thread Jean-Louis Abitbol
Dear Listers, I have the same problem on vista home premium. library(RWinEdt) Does not launch R-Winedt startWinEdt(.gW$InstallRoot, .gW$ApplData) Does nor either launch sessionInfo() R version 2.6.0 Patched (2007-10-03 r43075) i386-pc-mingw32 locale:

Re: [R] list matching

2007-10-07 Thread Adrian Dusa
On Sunday 07 October 2007, Gabor Grothendieck wrote: zoo's merge can do a multiway intersection. We turn each component of aa into the times of a dataless zoo object (assuming the elements of each component are unique) and merge them together using all = FALSE which will only leave those

[R] Controlling the Output from table()

2007-10-07 Thread Öhagen Patrik
Dear All, I would like to have some controll over the output from the function table(). I want to controll the order of the cells and I want to include empty cells (if any). Example: I have ordinal data wich takes the values 1,2,3N but I want the output from table to put the

Re: [R] Controlling the Output from table()

2007-10-07 Thread Prof Brian Ripley
On Sun, 7 Oct 2007, Öhagen Patrik wrote: Dear All, I would like to have some controll over the output from the function table(). I want to controll the order of the cells and I want to include empty cells (if any). Example: I have ordinal data wich takes the values 1,2,3N but I want

Re: [R] Controlling the Output from table()

2007-10-07 Thread Öhagen Patrik
Sorry for the vaugness. I have text data (exam grades, say) for different courses and I want to table them. In my table I want the frequency for the lowest grade first and the other frequencies in increaseing (grade) order. Some frequencies might be zero but I want to include those empty cells

Re: [R] Simulate data based on correlation coefficient

2007-10-07 Thread joris . dewolf
Try this: rnormcor - function(x,rho) rnorm(1,rho*x,sqrt(1-rho^2)) a - rnorm(1000,0,1) b - sapply(a, rnormcor, rho = 0.7) cor(a,b) Joris Gustavo

Re: [R] how to use \Sexpr{} with sweave

2007-10-07 Thread Friedrich Schuster
Hi, this works for me (R 2.5.1, Linux): echo=false= data(airquality) library(ctest) kruskal.test(Ozone ~ Month, data = airquality) @ Wir haben hier also \Sexpr{nrow(airquality)} Datenzeilen in airquality. This works \Sexpr{2+6} I think the error is in your latex setup, not in R or Sweave.

Re: [R] Re-ordering factors

2007-10-07 Thread maj
Thanks to James and Phil and Peter for their helpful suggestions. I think that I should also point out one way *not* to do the job: xtabs(Count ~ Education + Age_Group, data=educ) Age_Group Education64 25-34 35-44 45-54 55-64 CompletedHS 7558 16431 1855 9435 8795

Re: [R] Tricky vectorization problem

2007-10-07 Thread Dimitris Rizopoulos
a small improvement is to avoid computing the spectral decomposition of `Sigma' (look at code of mvrnorm()) in each iteration, e.g., #set up a collector for subject means in A a - numeric(Ns) #set up a collector for subject means in B b - numeric(Ns) eS - eigen(Sigma, symmetric = TRUE, EISPACK

Re: [R] how to use \Sexpr{} with sweave

2007-10-07 Thread Duncan Murdoch
[EMAIL PROTECTED] wrote: I'm trying to learn Sweave. So far things are going well with the chunks of code identified by = But I'm having trouble with the in-line text use of \Sexpr. Here is a short example .Rnw file: \documentclass[12pt]{article} \usepackage[margin=1.25in]{geometry}

Re: [R] R-2.6.0 and RWinEdt

2007-10-07 Thread Uwe Ligges
Again, a binary includeing a bugfix for R-2.6.0 of RWinEdt is on its way to CRAN master. Uwe Ligges Jean-Louis Abitbol wrote: Dear Listers, I have the same problem on vista home premium. library(RWinEdt) Does not launch R-Winedt startWinEdt(.gW$InstallRoot, .gW$ApplData) Does nor

Re: [R] Variable characterization by the modalities of other variables

2007-10-07 Thread Jorge Manuel de Almeida Magalhães
sorry to all: I'm looking for a way to display the statistical links between categories of observations and active variables (or characteristic modalities). SAPD call this Description of Modalities - procedure DEMOD which statistically characterises each category of the dependent variable

Re: [R] list matching

2007-10-07 Thread Gabor Grothendieck
Perhaps its a version problem? Here is what I get: aa - list(one=c(o, n, e), + tea=c(t, e, a), + thre=c(t, h, r, e)) library(zoo) as.vector(time(do.call(merge, c(lapply(aa, function(x) zoo(,x)), all = FALSE [1] e R.version.string # Vista [1] R version 2.6.0

Re: [R] factanal: error message

2007-10-07 Thread Michael Dewey
At 17:51 06/10/2007, Steve Friedman wrote: Hi everyone, I'm running a factor analysis on a correlation matrix with 32 rows and columns. I get the following error when I issue the command sequence mich.fac1 - factanal(michcor, factor=1) Error in solve.default(cv) : system is computationally

Re: [R] plot the chi-square distribution for n=1 through 10, on one plot

2007-10-07 Thread Michael Dewey
At 23:48 06/10/2007, fb wrote: Dear help list. I would like to plot the chi-square distribution for n=1 through 10, on one plot. How do I go about doing that? I know how to plot one chi-square (using curve() for example), but plotting 10 on one plot, that is too difficult for me. I greatly

Re: [R] list matching

2007-10-07 Thread Adrian Dusa
On Sunday 07 October 2007, you wrote: Perhaps its a version problem? [...] R.version.string # Vista [1] R version 2.6.0 beta (2007-09-23 r42958) packageDescription(zoo)$Version [1] 1.3-2 I don't know, everything is up to date here, too. I'm using Kubuntu Linux and the latest versions

Re: [R] Tricky vectorization problem

2007-10-07 Thread Mike Lawrence
Had a realization this morning that I think achieves ceiling performance and thought I'd share with the list. I was previously generating sample data per participant and then calculating a mean, but of course this could be simplified by simply getting a single value from the sampling

[R] ex-post prediction

2007-10-07 Thread erinys
Dear help-list, I would like to predict ex-post from a model fitted by arima, but the time series prediction methods in package stats have only the argument n.ahead specifying how many time steps ahead to predict . Is there any other possibility to do that? I really appreciate your help!

Re: [R] list matching

2007-10-07 Thread Gabor Grothendieck
I thought the latest zoo fixes were on CRAN but perhaps they are not. Try this: library(zoo) # next line loads latest version of merge.zoo (fixed in August) source(http://r-forge.r-project.org/plugins/scmsvn/viewcvs.php/*checkout*/pkg/R/merge.zoo.R?rev=361root=zoo;)

Re: [R] list matching

2007-10-07 Thread Adrian Dusa
On Sunday 07 October 2007, Gabor Grothendieck wrote: I thought the latest zoo fixes were on CRAN but perhaps they are not. Try this: library(zoo) # next line loads latest version of merge.zoo (fixed in August) source(http://r-forge.r-project.org/plugins/scmsvn/viewcvs.php/*checkout*/

Re: [R] list matching

2007-10-07 Thread Gabor Grothendieck
On 10/7/07, Adrian Dusa [EMAIL PROTECTED] wrote: On Sunday 07 October 2007, Gabor Grothendieck wrote: I thought the latest zoo fixes were on CRAN but perhaps they are not. Try this: library(zoo) # next line loads latest version of merge.zoo (fixed in August)

[R] Fwd: cluster analysis

2007-10-07 Thread amna khan
-- Forwarded message -- From: amna khan [EMAIL PROTECTED] Date: Oct 6, 2007 11:02 AM Subject: cluster analysis To: Alberto Viglione [EMAIL PROTECTED] Dear Sir I am facing problem in using traceWminim() function for cluster analysis. A sample of data frame that I am using is of

Re: [R] R-2.6.0 and RWinEdt

2007-10-07 Thread Uwe Ligges
Patrick Giraudoux wrote: Dear Listers, I have just installed R-2.6.0 and the RWinEdt package 1.7-6 under Windows XP. wait for 1.7-7 which should appear on CRAN real soon now. Uwe The R-WinEdt menu well appears at launching (the command library(RWinEdt) is in .Rprofile), but WinEdt

[R] Trust p-values or not ?

2007-10-07 Thread rocker turtle
Hi, First of all kudos to the creaters/contributors to R ! This is a great package and I am finding it very useful in my research, will love to contribute any modules and dataset which I develop to the project. While doing multiple regression I arrived at the following peculiar situation. Out of

[R] Question about aov

2007-10-07 Thread Alexander Keedy
Hello R gurus, I am a beginner with R. I am doing an ANCOVA analysis using 'aov,' and need some help understanding how 'aov' works. I have a dataset (taken from http://faculty.vassar.edu/lowry/ch17pt2.html) looking at hypnotic induction. The variable 'X' is a measure of how susceptible

Re: [R] how to use \Sexpr{} with sweave

2007-10-07 Thread cryan
Duncan-- The .tex file that you describe as output of Sweave(test.Rnw) is what I had expected. But I get this .tex file when I run Sweave: \documentclass[12pt]{article} \usepackage[margin=1.25in]{geometry} \usepackage{graphicx} \usepackage{Sweave} \begin{document} \begin{Schunk}

[R] Errors in R-2-6-0?

2007-10-07 Thread francogrex
I was happy today to install the new version of R 2-6-0 But I ran into problems I did not have before: setwd(C:/R/DATA/BRugs) library(BRugs) modelCheck(rcapturemodel.txt) Error in file(con, r) : unable to open connection In addition: Warning message: In file(con, r) : cannot open file

Re: [R] Linux editor for R+LaTeX, but not Emacs

2007-10-07 Thread Patrick Connolly
On Sat, 06-Oct-2007 at 11:01AM +0200, Scionforbai wrote: | Hi, | | What features are you missing in emacs that you wish were there? Are | these ESS features or LaTeX related features? | | is it only me or has anyone else the problem that running an R process | within emacs is way much slower

Re: [R] how to use \Sexpr{} with sweave

2007-10-07 Thread cryan
I think some inadvertent cutting and pasting may have made my last e-mail difficult to interpret. Here is the .tex file I get as output of Sweaving my .Rnw file. Note that \Sexpr{2+6} appears where you get (and I want) simply 8. \documentclass[12pt]{article}

Re: [R] Errors in R-2-6-0?

2007-10-07 Thread francogrex
I re-installed version 2-5-0 and it works well. I believe this problem is due to an incompatibility of BRugs with R-2-6-0 francogrex wrote: I was happy today to install the new version of R 2-6-0 But I ran into problems I did not have before: setwd(C:/R/DATA/BRugs) library(BRugs)

Re: [R] Trust p-values or not ?

2007-10-07 Thread John Sorkin
Ankit, (1) Not necessarily. Linear regression has a number of assumptions. I suggest you get a basics statistics textbook and do some reading. A brief summary of the assumptions include: (a) The relation between outcome and predictor variables lie along a line (or plane for a regression with

[R] Arguments to personalised plot()

2007-10-07 Thread Ted Harding
Hi Folks, I'm curious for an explanation of the following -- it's a matter of trying to understand how R parses it. I've written sundry little helper variants of functions, in particular plot(), to save repetitively typing the same options over and over again. For example: plotb -

Re: [R] Linux editor for R+LaTeX, but not Emacs

2007-10-07 Thread Frede Aakmann Tøgersen
Have you ever considered using C-u C-c C-r instead of C-c C-r? Best regards Frede Fra: [EMAIL PROTECTED] på vegne af [EMAIL PROTECTED] Sendt: lø 06-10-2007 23:10 Til: Dirk Eddelbuettel Cc: r-help@r-project.org Emne: Re: [R] Linux editor for R+LaTeX, but not

Re: [R] Trust p-values or not ?

2007-10-07 Thread Peter Dalgaard
rocker turtle wrote: Hi, First of all kudos to the creaters/contributors to R ! This is a great package and I am finding it very useful in my research, will love to contribute any modules and dataset which I develop to the project. While doing multiple regression I arrived at the following

Re: [R] Linux editor for R+LaTeX, but not Emacs

2007-10-07 Thread Deepayan Sarkar
On 10/7/07, Frede Aakmann Tøgersen [EMAIL PROTECTED] wrote: Have you ever considered using C-u C-c C-r instead of C-c C-r? Yes (or equivalently (setq ess-eval-visibly-p nil)), and found that I didn't like it much. In any case, I'm not looking for ugly workarounds, I just wanted to point out that

[R] Old packages to updated R

2007-10-07 Thread Jon Loehrke
I am 100% positive there is a better fix for this problem and it has been discussed. Thank you in advance for bringing it to my attention. When upgrading to newer versions of R the old packages are not brought along in the process. I solved that problem by creating a list of the packages in

Re: [R] Arguments to personalised plot()

2007-10-07 Thread Ted Harding
On 07-Oct-07 21:01:08, Deepayan Sarkar wrote: On 10/7/07, Ted Harding [EMAIL PROTECTED] wrote: Hi Folks, I'm curious for an explanation of the following -- it's a matter of trying to understand how R parses it. I've written sundry little helper variants of functions, in particular plot(),

Re: [R] Arguments to personalised plot()

2007-10-07 Thread Rolf Turner
On 8/10/2007, at 10:01 AM, Deepayan Sarkar wrote: On 10/7/07, Ted Harding [EMAIL PROTECTED] wrote: Hi Folks, I'm curious for an explanation of the following -- it's a matter of trying to understand how R parses it. I've written sundry little helper variants of functions, in particular

Re: [R] Question about aov

2007-10-07 Thread John Fox
Dear Alex, First, the call to lm() is superfluous; you'd get the same output from summary(aov(Y~X+Method,data=hypnotic.induction)). Second, it's not common to use aov() in this manner, since the function is really intended for balanced ANOVAs. More common would be something like model -

Re: [R] Tricky vectorization problem

2007-10-07 Thread Mike Lawrence
Posting this for posterity and to demonstrate that a speed-up of 2 orders of magnitude is indeed possible! I can now run 1e5 monte carlo experiments in the time the old code could only run 1e3. The final change was to replace my use of cor() with .Internal(cor()), which avoids some checks

Re: [R] Arguments to personalised plot()

2007-10-07 Thread Ted Harding
On 07-Oct-07 21:51:29, Rolf Turner wrote: On 8/10/2007, at 10:01 AM, Deepayan Sarkar wrote: [...] I'm not sure exactly which part you didn't expect. Given this definition of plotb, I would expect plotb(x,cos(x)) to expand to plot(x, pch=+, col=blue, cos(x)) and as far as I can tell,

Re: [R] Break [SEC=UNCLASSIFIED]

2007-10-07 Thread Crombie, Joe
Hi, You can exit out of the 'identify()' routine by either: right-click 'stop' Or click on 'stop stop locator' in the top-left of the plot window. The script should continue with line1, line2 etc. Cheers Joe Joe Crombie Information and Risk Sciences Bureau of Rural Science Canberra

[R] multcomp and lme4

2007-10-07 Thread Andrew Dolman
Dear list members, Can anyone please point to an example of how to use glht(multcomp) with lmer objects? I am trying: summary(glht(lmerObject, linfct = mcp(x = Tukey))) as I would for a glm object, but with no luck. Thank you, Andy. [[alternative HTML version deleted]]

[R] R and FDA trials

2007-10-07 Thread Ricardo Pietrobon
Yesterday I just noticed the new document on R and regulatory aspects for biomedical research posted at http://www.r-project.org/doc/R-FDA.pdf Coming from an institution that performs a large number of clinical trials for FDA and being an advocate of R myself, I have found that the following

Re: [R] R and FDA trials

2007-10-07 Thread Frank E Harrell Jr
Ricardo Pietrobon wrote: Yesterday I just noticed the new document on R and regulatory aspects for biomedical research posted at http://www.r-project.org/doc/R-FDA.pdf Coming from an institution that performs a large number of clinical trials for FDA and being an advocate of R myself, I

[R] How can i perform Fisher Exact Test for such data in R?

2007-10-07 Thread Jiantao Shi
Dear all, I am sorry, but this not exact a R question.It's a statistic one. Suppose I have two vectors PA and PB, which are both consist of no-redundant nature numbers. PA and PB have overlaps.GA and GB are drawn from PA and PB, respectively.So GA and GB may have overlaps. The question is whether

[R] special characters in linux using dev.print

2007-10-07 Thread Alexander.Herr
Hi List, I want to print ± in a lattice graph label. This works in windows, but the linux version has problems - it cannot translate the character. Error is invalid input in mbcsToLatin1. I use the standard encoding and also dev.print(file=filename,dev=pdf, encoding=PDFDoc.enc) I am afraid I