Re: [R] problem with jpeg and Rbitmap.dll

2010-08-31 Thread Prof Brian Ripley
On Tue, 31 Aug 2010, Erin Hodgess wrote: Dear R People: I am having trouble saving a plot to jpeg via the graphics device. I put it: plot(1:10) and click on the graphics device. I save it to jpeg at 100% quality, put the file name in, and then receive an "Impossible to load Rbitmapp.dll" me

[R] OdfWeave and Locale

2010-08-31 Thread C.H.
Dear all, I have a question regarding the odfweave. I created an odt file (test_input.odt) using OO.o for WIndows XP. The code is generic. Code for odt --- \Sexpr{dim(iris)[1]} <>= odfItemize(levels(iris$Species)) @ Table 1: <>= data.frame(N =tapply(iris$Petal.Length, iris$Speci

Re: [R] non-linear plot parameters

2010-08-31 Thread Marlin Keith Cox
OK, I am thankful and understand what my mistakes were, but now I look at my real data and have tried self starting functions and creating my own functions, resulting in other errors. This is one data set. I have used linear models, but as these data reflect fish consuming energy during a starvat

Re: [R] Can WinBUGS run two processes through R at the same time?

2010-08-31 Thread Uwe Ligges
R can, but you need to specify different working directories if you want to save the inits. On 01.09.2010 02:20, Dejian Zhao wrote: Does winbugs write anything to the file 'inits1.txt'? If yes,possibly the file is exclusively occupied by the first winbugs process. On 2010-9-1 1:12, 潘家群 wrote:

Re: [R] Can WinBUGS run two processes through R at the same time?

2010-08-31 Thread Uwe Ligges
R can, but you need to specify different working directories if you want to save the inits. Uwe Ligges On 31.08.2010 11:25, 潘家群 wrote: Dear all, I want to run two winBugs process through R (by R2WinBUGS package). The first WinBUGS process can successfully perform under R. The error messag

[R] problem with jpeg and Rbitmap.dll

2010-08-31 Thread Erin Hodgess
Dear R People: I am having trouble saving a plot to jpeg via the graphics device. I put it: plot(1:10) and click on the graphics device. I save it to jpeg at 100% quality, put the file name in, and then receive an "Impossible to load Rbitmapp.dll" message. This is on Windows XP on a new Netbo

[R] transaction object - how to coerce this data

2010-08-31 Thread pdb
Hi, I am wanting to look at frequent item sets using the arules package. I need to transform my data into a "transactions" object. The data I read in from a file has 2 columns, an ID and an item. How do I convert data like this into a transactions object? I've tried class? transactions but it o

Re: [R] "mode" function

2010-08-31 Thread David Winsemius
On Aug 31, 2010, at 10:48 PM, Dejian Zhao wrote: R does have a mode function, but it seems NOT to do the same thing as in matlab. Correct. It does not return the value of a vector with maximal proportions of individual elements. ?mode (It appears that the Matlab formalism fn(Name[vec, :

Re: [R] weird to me interaction between time() and %%, %/%

2010-08-31 Thread Gabor Grothendieck
On Tue, Aug 31, 2010 at 7:04 PM, Alistair Gray wrote: > Dear List, > I'm getting weird and unexpected behaviour using time and %%, or %/%. It's > likely I'm not appreciating the nuances of floating point arithmetic. Or it > could be a bug. > Here are two workarounds to get the year time(junk)

Re: [R] vglm

2010-08-31 Thread Dejian Zhao
try fit=vglm(mydata[,"Loss"]~1,pareto1(location=alpha),trace=TRUE,crit="c") On 2010-9-1 3:20, choonhong ang wrote: Hi All, could anybody help me to understand what is this error means ? mydata=read.table("C:/Documents and Settings/angieb/Desktop/CommercialGL/cl_ilf_claimdata.csv",header=TRUE,

[R] problems with JGR and standard R GUI closing during calculations

2010-08-31 Thread Moore, Meagan
Hi R help, I am having problems running calculations in both the standard R GUI (ver 2.11.0) and the Java GUI for R (1.7-2) in Windows XP Service Pack 3 (2.2 GHz processor with 2 GB of RAM). When I try to import my data and run calculations on it, the GUI closes within a minute, sometimes almo

Re: [R] "mode" function

2010-08-31 Thread Dejian Zhao
R does have a mode function, but it seems NOT to do the same thing as in matlab. > A <- matrix(c(1:3,1,3,2,3,3,2),nrow=3,byrow=F) > A [,1] [,2] [,3] [1,]113 [2,]233 [3,]322 > A[,2] [1] 1 3 2 > mode(A) [1] "numeric" > mode(A[,2]) [1] "numeric" On 2010-9-

Re: [R] "mode" function

2010-08-31 Thread David Winsemius
On Aug 31, 2010, at 5:19 PM, Iain Martyn wrote: Hi, I have what I think is a simple question/issue but I have unable to find the answer for it either in the R-manual or by browsing the web. I would like to know if there is a "mode" function in R, such that from a vector or matrix the

Re: [R] weird to me interaction between time() and %%, %/%

2010-08-31 Thread David Scott
Steve is undoubtedly right about this being a floating point problem. One comment I would make is that I think floor is a bit nicer for this, both conceptually and for readability. Also you don't need a really small value to add, half the sampling interval will work: > junk <- ts(1:125, start=

Re: [R] Detecting Growth Trends

2010-08-31 Thread rtist
If the test rejects the null, then it has determined that the new set of incoming data is no longer purely oscillatory in the mean reverting sense (it is now unit root and exhibits growth). Unless I misinterpreted, the OP wants to find a statistical method to determine such behavior beyond purely

[R] "mode" function

2010-08-31 Thread Iain Martyn
Hi, I have what I think is a simple question/issue but I have unable to find the answer for it either in the R-manual or by browsing the web. I would like to know if there is a "mode" function in R, such that from a vector or matrix the function returns the most common value. In other progr

Re: [R] "pairs" with same xlim and ylim scale

2010-08-31 Thread Dejian Zhao
I think you have successfully passed the "xlim" and "ylim" into the function pairs1. Compare the two graphs produced by the codes you provided, you can find the xlim and ylim in the second graph have been reset to the assigned value. It seems that the program halted in producing the second plot

[R] Collapsing levels of categorical variables

2010-08-31 Thread Saeed Abu Nimeh
In this paper [1] the author mentioned a procedure by M. Greenace that can be used to collapse the levels of a categorical variable by setting up a table with the frequency of each level and the proportion of the target value in each level. Then collapsing the table level by level looking at the ch

Re: [R] 3D tomography data

2010-08-31 Thread Kjetil Halvorsen
I guess you need to be more specific about what you want to do. But a search with RSiteSearch("Tomography") leeds to one package names RTOMO, which you can install by calling install.packages("RTOMO",dep=T) then type library(RTOMO) library(help=RTOMO) kjetil On Tue, Aug 31, 2010 at 9:41 AM,

Re: [R] Can WinBUGS run two processes through R at the same time?

2010-08-31 Thread Dejian Zhao
Does winbugs write anything to the file 'inits1.txt'? If yes,possibly the file is exclusively occupied by the first winbugs process. On 2010-9-1 1:12, 潘家群 wrote: Dear all, I want to run two winBugs processes through R (by R2WinBUGS package). The first WinBUGS process can successfully perfor

Re: [R] lattice help required

2010-08-31 Thread Felix Andrews
The development version of the lattice package has a 'ylab.right' option, so you can do something like this to place your y labels: xyplot(1:10 ~ 1:10 | c('a','b'), layout = c(1,2), ylab = list("label one", y = 0.75), ylab.right = list("label two", y = 0.25, rot = 270)) For the row-specif

Re: [R] weird to me interaction between time() and %%, %/%

2010-08-31 Thread Alistair Gray
Steve, thanks for the suggestion. Alistair On 01/09/10 11:20, Steve Taylor wrote: Yes it is a floating point problem. January 2005 is slightly less than 2005, although the gap is somewhat bigger than .Machine$double.eps time(junk) - 2005 # shows the gap is -2.2737e-13 for my machine Try somethi

Re: [R] Retrieving predictions after fitting a model

2010-08-31 Thread Ben Bolker
alexander pswarayi ualberta.ca> writes: > > Hi there > > I am trying to retrieve the predicted values after fitting a model (dose > response model). How do I do it. What is the code for it? For estimated > parameters, one uses : summary(filename). And for the predicted values, what > do I use?

Re: [R] weird to me interaction between time() and %%, %/%

2010-08-31 Thread Steve Taylor
Yes it is a floating point problem. January 2005 is slightly less than 2005, although the gap is somewhat bigger than .Machine$double.eps time(junk) - 2005 # shows the gap is -2.2737e-13 for my machine Try something this: floor(time(junk) + 1e-10) — — — — — — — — —

Re: [R] Reading multiple text files where some files are empty

2010-08-31 Thread Paul
Thanks for the responses. 'Try' seems to be the function that I needed. I had to adapt the code suggested below because vapply wasn't recognised, but this seemed to work: a<-list.files() x<-lapply(a,function(x){try(read.table(x,colClasses='character'))}) bad <- sapply(x, function(x){inherits(x, "

[R] weird to me interaction between time() and %%, %/%

2010-08-31 Thread Alistair Gray
Dear List, I'm getting weird and unexpected behaviour using time and %%, or %/%. It's likely I'm not appreciating the nuances of floating point arithmetic. Or it could be a bug. I'm running > R.version _ platform x86_64-redhat-linux-gnu kernel 2.6.33.8-149.fc13.x86_64

Re: [R] Question regarding significance of a covariate in a coxme survival

2010-08-31 Thread Norm Matloff
In 2010-08-30, C. Peng wrote: > What statistical measure(s) tend to be answering ALL(?) question of > practical interest? None. All I had said was that significance testing doesn't really answer any questions of practical interest. Unfortunately, that doesn't mean there's something to answer a

[R] Retrieving predictions after fitting a model

2010-08-31 Thread alexander pswarayi
Hi there I am trying to retrieve the predicted values after fitting a model (dose response model). How do I do it. What is the code for it? For estimated parameters, one uses : summary(filename). And for the predicted values, what do I use? Many thanks. Alex [[alternative HTML version de

Re: [R] Brown-Forsythe test of equality of MEANS

2010-08-31 Thread JLucke
Yes. I too remember BMD, BMDP, punch cards, computer printouts, and 24-hour turn-around. However, I am also old enough that when someone asks for an antique method, I will gladly supply it if I can, if only for historical reasons. Bert Gunter 08/31/2010 04:41 PM To jlu...@ria.buffalo.edu

Re: [R] Brown-Forsythe test of equality of MEANS

2010-08-31 Thread Iasonas Lamprianou
Combining humor with useful information is nice! Thank you for the advice Bert. I'll try it and see what happens. Are you aware of any free! (on the internet) alternatives to the book you have recommended? Academics are not, usually, known for their deep pockets! Thank you again for the suggesti

Re: [R] Detecting Growth Trends

2010-08-31 Thread Bert Gunter
Error of the 3rd kind, (right answer to wrong question), I think. So what if the test rejects --- then what? I think the poster is looking for some kind of smoother. ?loess, ?smooth.spline and about 400 others may be useful. -- Bert Gunter Genentech Nonclinical Statistics On Tue, Aug 31, 2010 a

Re: [R] Brown-Forsythe test of equality of MEANS

2010-08-31 Thread Bert Gunter
Learned Folks: Well, I've already advertised my ignorance about these matters, so I have nothing to lose by plunging ahead with further Questionable "advice." >From the references cited, Brown-Forsythe originated in the statistical medieval age -- that is, prior to large scale, cheap computing (t

Re: [R] lattice help required

2010-08-31 Thread Kay Cichini
...after digging for a while now, i came to a solution which may be crude, but it works perfectly well. if someone would advise me how to use textGrob avoiding the "\n" in the title, combined with the linewidth = 4 in gpar(), to put the axis title into the right place i'd be very glad! maybe ther

Re: [R] Brown-Forsythe test of equality of MEANS

2010-08-31 Thread Iasonas Lamprianou
Thank you. It seems that nobody has implemented it in R yet. I'll have a look, thank you very much Dr. Iasonas Lamprianou Assistant Professor (Educational Research and Evaluation) Department of Education Sciences European University-Cyprus P.O. Box 22006 1516 Nicosia Cyprus Tel.: +357

Re: [R] any statement equals to 'goto'?

2010-08-31 Thread karena
these seem something that I am looking for, I will try them, thank you!! -- View this message in context: http://r.789695.n4.nabble.com/any-statement-equals-to-goto-tp2402107p2402264.html Sent from the R help mailing list archive at Nabble.com. __ R-h

Re: [R] any statement equals to 'goto'?

2010-08-31 Thread Allan Engelhardt
Sounds like you want help("tryCatch") # Catch the error and do something help("next") # Go to the next value of the surrounding loop Hope this helps a little. Allan On 31/08/10 19:34, karena wrote: I have the following code: ---

Re: [R] vglm

2010-08-31 Thread Joshua Wiley
Hi, My guess is that vglm() cannot find "Loss" because you have not set the data argument. Something like: fit = vglm(Loss ~ 1, pareto1(location = alpha), trace = TRUE, crit = "c", data = mydata) out to work if "Loss" is a variable in the object "mydata" HTH, Josh On Tue, Aug 31, 2010 at 12

Re: [R] vglm

2010-08-31 Thread David Winsemius
On Aug 31, 2010, at 3:20 PM, choonhong ang wrote: Hi All, could anybody help me to understand what is this error means ? mydata=read.table("C:/Documents and Settings/angieb/Desktop/CommercialGL/ cl_ilf_claimdata.csv",header=TRUE,sep=",") names(mydata) [1] "ILFTable""liabLimit" "AnnA

Re: [R] help - read SAS into R

2010-08-31 Thread David Winsemius
On Aug 31, 2010, at 1:25 PM, choonhong ang wrote: HI All, How to read SAS data directly into R ? If the data is in the transport format, then foreign::read.xport has worked well for me. -- David Winsemius, MD West Hartford, CT __ R-help@r-pr

Re: [R] odfWeave

2010-08-31 Thread John Kane
Replying to my own post, it is a path problem. I downloaded the UnZip and Zip binaries, the 5.51-1 UnZip and the most recent Zip and then added a path to the UnZip directory, and odfWeave seems to be working. --- On T>e, 8/31/10, John Kane wrote: > From: John Kane > Subject: Re: [R] odfWeave

Re: [R] Making plots in big scatterplot matrix large enough to see

2010-08-31 Thread Greg Snow
Look at the pairs2 function in the TeachingDemos package, this lets you produce smaller portions of the total scatterplot matrix at a time (with bigger plots), you could print the smaller portions then assemble the full matrix on a large wall, or just use it to look at potentially interesting pa

Re: [R] summary in Hmisc and Latex

2010-08-31 Thread Marc Schwartz
On Aug 31, 2010, at 12:42 PM, moleps wrote: > Dear all, > > With the latest update of Hmisc I no longer have any problems with latex. > However using the ctable option produces latex code that at least on both the > miktex distribution at work and mactex distribution at home refuses to run > d

Re: [R] About plot graphs

2010-08-31 Thread Greg Snow
Stephen, Gavin replied to the specific questions, so I won't repeat what he said. But I think that you really need to understand the search path and scoping better. Remember that computers are stupid, they do exactly what you tell them to do (see fortune(226)). R is such a amazing program and

Re: [R] Brown-Forsythe test of equality of MEANS

2010-08-31 Thread JLucke
The following reference that contains a short Fortran program for the Brown-Forsythe ANOVA Reed, James F., I. & Stark, D. B. Robust alternatives to traditional analyses of variance: Welch $W^*$, James $J_I^*$, James $J_II^*$, and Brown-Forsythe $BF^*$ Computer Methods and Programs in Biomedicine

[R] vglm

2010-08-31 Thread choonhong ang
Hi All, could anybody help me to understand what is this error means ? mydata=read.table("C:/Documents and Settings/angieb/Desktop/CommercialGL/cl_ilf_claimdata.csv",header=TRUE,sep=",") > names(mydata) [1] "ILFTable""liabLimit" "AnnAggLimit" "DedAmt" "Loss" "TIL" > fit=vglm(Loss~1,pa

Re: [R] Detecting Growth Trends

2010-08-31 Thread rtist
You can try a unit root test which test for stationarity in a series. -- View this message in context: http://r.789695.n4.nabble.com/Detecting-Growth-Trends-tp2402080p2402132.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-pr

[R] summary in Hmisc and Latex

2010-08-31 Thread moleps
Dear all, With the latest update of Hmisc I no longer have any problems with latex. However using the ctable option produces latex code that at least on both the miktex distribution at work and mactex distribution at home refuses to run due to an extra blank line inserted between the multicolum

Re: [R] How to Adaptively Set Up the Coordinate Range of Multiple Graphs in One Figure

2010-08-31 Thread Greg Snow
The best approach is to calculate the limits before doing any plotting and set them appropriately. You could use your loop to compute y, but instead of calling lines, store the values, then use a function like matplot (or ggplot2 or lattice) to do the plotting. If you really need to change the

[R] any statement equals to 'goto'?

2010-08-31 Thread karena
I have the following code: - result <- matrix(NA, nrow=1, ncol=5) for(i in 1:(nsnp-1)) { for(j in (i+1):nsnp){ tempsnp1 <- data.lme[,i] tempsnp2 <- data.lme[,j] fm1 <- lme(trait~sex+age+rmtemp.b+f

Re: [R] anova and lm results differ

2010-08-31 Thread Iasonas Lamprianou
So, should I avoid using the Anova option of the Rcmdr and use aov instead (as I have been doing)? Dr. Iasonas Lamprianou Assistant Professor (Educational Research and Evaluation) Department of Education Sciences European University-Cyprus P.O. Box 22006 1516 Nicosia Cyprus Tel.: +357-22-71317

[R] "pairs" with same xlim and ylim scale

2010-08-31 Thread Shi, Tao
Hi list, I have a function which basically is a wrapper of pairs with some useful panel functions. However, I'm having trouble to pass the "xlim" and "ylim" into the function so the x and y axes are in the same scale and 45 degree lines are exactly diagonal. I've looked at some old posts, th

Re: [R] odfWeave

2010-08-31 Thread John Kane
It looks like a path problem. I have not figured out how to handle it in Win 7 (new system last week) but have a look at http://r.789695.n4.nabble.com/odfWeave-Error-td1595848.html#a1595848 which I think addresses the problem. --- On Tue, 8/31/10, Michael Kubovy wrote: > From: Michael Kubovy

[R] Detecting Growth Trends

2010-08-31 Thread Lorenzo Isella
Dear All, I am given some noisy data which (by naked eye) appears to be oscillating first but finally growing. Is there any statistical set (I mean something different from e.g. a linear fit, which would not be convincing at all in my case) to detect growth (possibly without relying on any data

Re: [R] help - read SAS into R

2010-08-31 Thread Nordlund, Dan (DSHS/RDA)
> -Original Message- > From: r-help-boun...@r-project.org [mailto:r-help-boun...@r- > project.org] On Behalf Of choonhong ang > Sent: Tuesday, August 31, 2010 10:25 AM > To: r-help@r-project.org > Subject: [R] help - read SAS into R > > HI All, > > How to read SAS data directly into R ? >

Re: [R] help - read SAS into R

2010-08-31 Thread Erik Iverson
Looks like google has many suggestions. I use Hmisc package's sas.get function, but that requires that you actually have SAS installed. If you do have access to SAS, I've found that to be the easiest. choonhong ang wrote: HI All, How to read SAS data directly into R ? Thank you [[

Re: [R] specify axis range

2010-08-31 Thread Gavin Simpson
On Tue, 2010-08-31 at 17:18 +0100, Lathouri, Maria wrote: > Dear all, > > I want to make a plot where I can specify the range of y axis. I am using the > interaction.plot command and because the values of y axis are from 2 to 500, > I am using a logarithmic scale. > > >interaction.plot(speclong

Re: [R] Greek letters and variable values in lattice plot annotation

2010-08-31 Thread Thomas Maier
Thanks! That does exactly what I needed. On Tue, Aug 31, 2010 at 12:07 PM, David Winsemius wrote: > > On Aug 31, 2010, at 11:34 AM, Thomas Maier wrote: > > I'd like to annotate a lattice plot, e.g. xyplot, with text containing >> both >> Greek letters and variable values. In the base graphics p

Re: [R] Where can I download/install grDevices

2010-08-31 Thread Gavin Simpson
On Tue, 2010-08-31 at 05:23 -0700, Stephen Liu wrote: > Hi Dirk, > > I already have the package installed. > > $ apt-cache policy r-base-core > r-base-core: > Installed: 2.7.1-1+lenny1 > Candidate: 2.7.1-1+lenny1 > Version table: > *** 2.7.1-1+lenny1 0 > 500 http://ftp.hk.debian.or

[R] Can WinBUGS run two processes through R at the same time?

2010-08-31 Thread 潘家群
Dear all, I want to run two winBugs processes through R (by R2WinBUGS package). The first WinBUGS process can successfully perform under R. The error message appears while I put the code to run for the second WinBUGS process. The error message: Error in file(file, ifelse(append, "a", "w")) :

[R] rpart - interpretation of results of tree on survival data

2010-08-31 Thread Ritwik Sinha
Hi All, I am fitting a tree to censored survival data using the rpart package and wanted to better understand the results. I am trying to interpret the output from the tree. I am interested in understanding what "yval" is for a survival tree. I see in the output of summary, the phrase "estimated

Re: [R] help - read SAS into R

2010-08-31 Thread stephen sefick
maybe package foreign - all though I don't have any experience with this. On Tue, Aug 31, 2010 at 12:25 PM, choonhong ang wrote: > HI All, > > How to read SAS data directly into R ? > > Thank you > >        [[alternative HTML version deleted]] > > __ >

[R] odfWeave

2010-08-31 Thread Michael Kubovy
Dear -r-helpers, I'm trying to get my students to use odfWeave. Windows users are having trouble. I would appreciate advice. Comment 1: Hello I have spent a while today trying to get odfWeave to be successful I installed the unzip and zip program but it has not changed the error message.I h

[R] help - read SAS into R

2010-08-31 Thread choonhong ang
HI All, How to read SAS data directly into R ? Thank you [[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 http://www.R-project.org/posting-g

Re: [R] anova and lm results differ

2010-08-31 Thread Iasonas Lamprianou
Thank you, it seems that I need to do some reading... Dr. Iasonas Lamprianou Assistant Professor (Educational Research and Evaluation) Department of Education Sciences European University-Cyprus P.O. Box 22006 1516 Nicosia Cyprus Tel.: +357-22-713178 Fax: +357-22-590539 Honorary Research Fello

Re: [R] Remove "Loading required package" message

2010-08-31 Thread Barry Rowlingson
On Tue, Aug 31, 2010 at 5:21 PM, Allan Engelhardt wrote: > > > On 27/08/10 10:19, Barry Rowlingson wrote: >> >> On Fri, Aug 27, 2010 at 9:09 AM, Sébastien Moretti >>  wrote: >>   [...] >>  require it quietly: >> >> >>> >>> require(sp,quietly=TRUE) >>> > > Doesn't work for me: > >> require("Matrix

Re: [R] simultaneous estimation

2010-08-31 Thread Murali.Menon
Bert, I expect you are correct, burrito notwithstanding (wasn't Taco Bell, was it? :-) The full model adds differences and lags, and incorporates non-zero covariances in the innovations. I only simplified to get an idea of how to implement in R. For anyone interested, I'm looking at the Balv

Re: [R] Remove "Loading required package" message

2010-08-31 Thread Allan Engelhardt
On 27/08/10 10:19, Barry Rowlingson wrote: On Fri, Aug 27, 2010 at 9:09 AM, Sébastien Moretti wrote: [...] require it quietly: require(sp,quietly=TRUE) Doesn't work for me: > require("Matrix", quietly=TRUE) Loading required package: lattice Attaching package: 'Matrix'

[R] TukeyHSD for binomial distribution

2010-08-31 Thread adriana1986
Hi, I am looking for a way to get contrast statements (similar to the output of the TukeyHSD function) for glm fit data with a binomial distribution: glm.fit=glm(y~(treat)*(geno),family=binomial) tukey.fit=TukeyHSD(aov(glm.fit)) My y is a dataframe of number successful and number failed repro

Re: [R] anova and lm results differ

2010-08-31 Thread Prof Brian Ripley
You have confused 'Anova', 'aov' and 'anova'. R has aov() and anova(), package car (called from Rcmdr) has Anova(). They are not the same things: Anova computes (so-called, by SAS) 'type-II or type-III analysis-of-variance tables', something deliberately not available in base R. See http://

Re: [R] Unable to run dev.new

2010-08-31 Thread Joshua Wiley
On Tue, Aug 31, 2010 at 9:10 AM, Stephen Liu wrote: > Hi Duncan, > > Ubuntu 10.04 > R version 2.10.1 (2009-12-14) > > Debian 504 > R version 2.7.1 (2008-06-23) > > > Both R were installed on repo.  What can I do?  Thanks When I was using Lenny (Debian 505 I think) I found this page really helpful

Re: [R] specify axis range

2010-08-31 Thread stephen sefick
If you would provide a reproducible example, even a dummy example, I could give it a shot. Otherwise look for something like "break". On Tue, Aug 31, 2010 at 11:18 AM, Lathouri, Maria wrote: > Dear all, > > I want to make a plot where I can specify the range of y axis. I am using the > interact

[R] specify axis range

2010-08-31 Thread Lathouri, Maria
Dear all, I want to make a plot where I can specify the range of y axis. I am using the interaction.plot command and because the values of y axis are from 2 to 500, I am using a logarithmic scale. >interaction.plot(speclong$Date,speclong$time2,speclong$outcome,xaxt="n",type="l",pch=20,xlab="",

Re: [R] simultaneous estimation

2010-08-31 Thread Bert Gunter
I would hazard the guess that this would be better estimated as a multivariate time series (e.g. AR1) in which the covariance between the two innovation components was NOT assumed to be 0 (nor were their variances assumed to be the same). The R time series task view lists packages to do this, but

Re: [R] Unable to run dev.new

2010-08-31 Thread Stephen Liu
Hi Duncan, Ubuntu 10.04 R version 2.10.1 (2009-12-14) Debian 504 R version 2.7.1 (2008-06-23) Both R were installed on repo. What can I do? Thanks B.R. Stephen L - Original Message From: Duncan Murdoch To: Stephen Liu Cc: "r-help@r-project.org" Sent: Tue, August 31, 2010 11

[R] anova and lm results differ

2010-08-31 Thread Iasonas Lamprianou
Dear all I have found that the two "equivalent" commands do not produce the same results. 1. (I wrote this command by hand, this is what I would do usually) >summary(aov(eduyrs ~ cntry * edf, data=ESS1)) Df Sum Sq Mean Sq F valuePr(>F) cntry 1

Re: [R] Greek letters and variable values in lattice plot annotation

2010-08-31 Thread David Winsemius
On Aug 31, 2010, at 11:34 AM, Thomas Maier wrote: I'd like to annotate a lattice plot, e.g. xyplot, with text containing both Greek letters and variable values. In the base graphics plot command this can be accomplished with e.g. main = substitute(paste(lambda,"=",var),list(var=var)), where

Re: [R] Unable to run dev.new

2010-08-31 Thread Duncan Murdoch
On 31/08/2010 11:52 AM, Stephen Liu wrote: Hi folks, Debian 504 64 bit. Have spent more than 2 hours on dev.new and can't make it to work. > dev.new(height = 6, width = 12) Error in dev.new(height = 6, width = 12) : unused argument(s) (height = 6, width = 12) r-base-core has been install

Re: [R] ggplot sub- and superscript for axis name

2010-08-31 Thread Ista Zahn
Hi Ben, Since you already know LaTeX, I would go with what you know. Try p<-p + scale_x_continuous(name='$\\gamma_{fi}$') library(tikzDevice) ## may too to install first tikz(file="foo.tex", standAlone=TRUE) print(p) dev.off() system("pdflatex foo.tex") Best, Ista On Tue, Aug 31, 2010 at 1:29

Re: [R] simultaneous estimation

2010-08-31 Thread Duncan Murdoch
On 31/08/2010 11:00 AM, David Winsemius wrote: On Aug 31, 2010, at 10:35 AM, > wrote: > Hi Duncan, > > Thanks for your response. > > Indeed, independent normal errors were what I had in mind. As for > variances, if I assume they are the same, would a 'pooled model' > apply in this case? Is

[R] Unable to run dev.new

2010-08-31 Thread Stephen Liu
Hi folks, Debian 504 64 bit. Have spent more than 2 hours on dev.new and can't make it to work. > dev.new(height = 6, width = 12) Error in dev.new(height = 6, width = 12) : unused argument(s) (height = 6, width = 12) r-base-core has been installed. The command should be correct. I have ru

[R] solved Re: online documentation only?

2010-08-31 Thread Bob McCall
The problem was I needed to "connect" even though I was offline and couldn't get online. In this case, "connect" apparently means to start the loopback device. After doing this, the docs work! Must be a quirk of win 7. Thanks for all the replies. Bob -- View this message in context: http://r

[R] Greek letters and variable values in lattice plot annotation

2010-08-31 Thread Thomas Maier
I'd like to annotate a lattice plot, e.g. xyplot, with text containing both Greek letters and variable values. In the base graphics plot command this can be accomplished with e.g. main = substitute(paste(lambda,"=",var),list(var=var)), where var is a variable that contains the value. When I try to

[R] wavelet parameters with sowas

2010-08-31 Thread angela
Hello everybody. My name is Angela. I'm doing wavelet using the sowas library. The problem I have is that I don't know how to choose the paremeters to describes differnt wavelet analysis. How I select the noctave, nvoice, w0, s0,...? I hope my question is clear enough. Thank you very much befor

Re: [R] Making plots in big scatterplot matrix large enough to see

2010-08-31 Thread William Revelle
Jocelyn, In a partial answer to your question, try setting gap=0 in the calls to pairs. This will make the plots closer together. (You might also find pairs.panels in the psych package useful, -- it implements one of the help examples for pairs to report the histogram on the diagonal and

Re: [R] simultaneous estimation

2010-08-31 Thread David Winsemius
On Aug 31, 2010, at 10:35 AM, > wrote: Hi Duncan, Thanks for your response. Indeed, independent normal errors were what I had in mind. As for variances, if I assume they are the same, would a 'pooled model' apply in this case? Is that equivalent to your suggestion of concatenating x(1

[R] VECTOR_ELT() can only be applied to a 'list', not a 'integer (gstat problem)

2010-08-31 Thread Sheila Greene
Hi there, I have a problem using gstat when attempting to fit a regression-kriging model to my data. I have successfully carried out a stepwise regression on my data and I have the model residuals. However when I try to apply a variogram to the residuals I get the following error message; *>>>

Re: [R] simultaneous estimation

2010-08-31 Thread Murali.Menon
Hi Duncan, Thanks for your response. Indeed, independent normal errors were what I had in mind. As for variances, if I assume they are the same, would a 'pooled model' apply in this case? Is that equivalent to your suggestion of concatenating x(1,t) and x(2,t)? Cheers, Murali -Original M

Re: [R] Cohen's f

2010-08-31 Thread Ben Bolker
Iasonas Lamprianou yahoo.com> writes: > > Dear all, > does anyone know how I can use any package in R to compute ANOVA's omega squared or Cohen's f? > > Dr. Iasonas Lamprianou > I don't know whether any of these do what you need, but have you browsed through the results of: library(sos) f

Re: [R] ggplot sub- and superscript for axis name

2010-08-31 Thread David Winsemius
On Aug 31, 2010, at 9:29 AM, Benoit Boulinguiez wrote: Hi all, For publication purpose, I require to label ggplot figures axes with sub- or superscript text. I tried several ways, but never worked so far, to mix character string, sub- or superscripting on it and even worse, mathematical

Re: [R] pasting together 2 character arrays

2010-08-31 Thread Henrique Dallazuanna
Try this: with(expand.grid(Array1, combn(Array2, 2, paste, collapse = '.')[-3]), paste(Var1, Var2, sep = '.')) On Tue, Aug 31, 2010 at 10:16 AM, Maas James Dr (MED) wrote: > If possible I would like to combine two different character arrays in > combinations > > Array1 <- c("height","weight","ag

Re: [R] pasting together 2 character arrays

2010-08-31 Thread Jan van der Laan
Jim, It is not completely clear how you want to handle the items in Array2, but perhaps something like the following does what you needs (or at least points you in the right direction): paste( rep(Array1, each=2), rep(Array2[1], times=8), rep(Array2[2:3], times=4), sep='.') Regard

Re: [R] pasting together 2 character arrays

2010-08-31 Thread David Winsemius
On Aug 31, 2010, at 9:16 AM, Maas James Dr (MED) wrote: If possible I would like to combine two different character arrays in combinations Array1 <- c("height","weight","age","sex") Array2 <- c("trt0","trt1","trt2") I would like to combine these two character vectors to end up with such

[R] 3D tomography data

2010-08-31 Thread Steve
Hi all, I was recently informed about R, as i need a program that can calculate the nearest neighbour in 3D tomography data with its vector. However, I am new to R and it isnt exactly intuitive. Does anyone know of any 3D tutorials that may help me get started? Cheers, Steve -- View this messag

[R] ggplot sub- and superscript for axis name

2010-08-31 Thread Benoit Boulinguiez
Hi all, For publication purpose, I require to label ggplot figures axes with sub- or superscript text. I tried several ways, but never worked so far, to mix character string, sub- or superscripting on it and even worse, mathematical symbols. Let say I want to write the LateX equivalent of

[R] pasting together 2 character arrays

2010-08-31 Thread Maas James Dr (MED)
If possible I would like to combine two different character arrays in combinations Array1 <- c("height","weight","age","sex") Array2 <- c("trt0","trt1","trt2") I would like to combine these two character vectors to end up with such ... Array3 "height.trt0.trt1" "height.trt0.trt2" "weight.trt0

[R] Cohen's f

2010-08-31 Thread Iasonas Lamprianou
Dear all, does anyone know how I can use any package in R to compute ANOVA's omega squared or Cohen's f? Dr. Iasonas Lamprianou Assistant Professor (Educational Research and Evaluation) Department of Education Sciences European University-Cyprus P.O. Box 22006 1516 Nicosia Cyprus Tel.: +357-2

[R] Short survey concerning the use of software engineering in the field of High Performance Computing

2010-08-31 Thread Markus Schmidberger
Already more than 2000 answers. Take the change and help to get a general overview about the use of software engineering. Best Markus -- Dear Colleagues, this is a short survey (21 questions that take about 10 minutes to an

Re: [R] cannot allocate vector of size 381.5 Mb

2010-08-31 Thread Adrienne Wootten
I would suggest two things here: check on the size of other object you may have stored in memory, and get rid of what you don't need. ? ls ? rm also, consider running garbage collection to help free up memory in R gc() I hope this helps! A On Tue, Aug 31, 2010 at 1:56 AM, rusers.sh wrote: >

[R] Munich R Course - Parallel Computing with R

2010-08-31 Thread Markus Schmidberger
Dear R users, at November 8. & 9., 2010, there will be an R course for 'Parallel Computing with R' in Munich at the super computer center. The course is part of the Munich R Course series and gives an introduction to parallel computing with R. Especially the R packages snow, snowfall, multicor

Re: [R] Speeding up prediction of survival estimates when using `survifit'

2010-08-31 Thread Frank Harrell
Frank E Harrell Jr Professor and ChairmanSchool of Medicine Department of Biostatistics Vanderbilt University On Mon, 30 Aug 2010, Ravi Varadhan wrote: Hi, I fit a Cox PH model to estimate the cause-specific hazards (in a competing risks setting). Then , I

  1   2   >