Re: [R] online documentation only?

Hi You shall also check setting of your browser, especially proxy. Sometimes help system interacts with particular setting and produce errors or does not start at all. Regards Petr r-help-boun...@r-project.org napsal dne 30.08.2010 22:16:17: Thanks for the replies. I tried help.start()

Re: [R] lost in vegan package

eva epitta at upatras.gr writes: Hi R Helpers, I'm still new to R and i experience many difficulties..I'm using vegan package (R version 2.11) trying to calculate checkerboard units for each species pair of a matrix. I've prepared the function: pair.checker=function (dataset)

Re: [R] Error: cannot allocate vector of size 198.4 Mb

2010/8/31 나여나 dllm...@hanmail.net Hi, All I have a problem of R memory space. I am getting Error: cannot allocate vector of size 198.4 Mb You may want to check circle 2 of the R inferno (found here: http://www.burns-stat.com/pages/Tutor/R_inferno.pdf). [[alternative

[R] binary package build error:target 'xf-' is not a directory

Hi, I built a package in linux and generated its binary using R CMD build --binary testpack. This generated testpack.tar.gz. I tried to install this package in windows using R CMD INSTALL package.tar.gz and I got the error *installing to library 'H:/R-2.11.1/library' *installing

[R] Export Results to a txt file applying a function to a list

Dear users, ***I have a function f to simulate data from a model (example below used only to show my problem) f-function(n,mean1){ a-matrix(rnorm(n, mean1 , sd = 1),ncol=5) b-matrix(runif(n),ncol=5) data-rbind(a,b) out-data out} *I want to simulate 1000 datasets (here only 5) so I

Re: [R] lost in vegan package

Jari Oksanen wrote: eva epitta at upatras.gr writes: Hi R Helpers, I'm still new to R and i experience many difficulties..I'm using vegan package (R version 2.11) trying to calculate checkerboard units for each species pair of a matrix. I've prepared the function:

Re: [R] getNodeSet - what am I doing wrong?

Thanks! but: library(XML) xmlDoc - xmlTreeParse(http://www.unimod.org/xml/unimod_tables.xml;) getNodeSet(xmlDoc,//x:modifications_row, x) Error in function (classes, fdef, mtable) : unable to find an inherited method for function saveXML, for signature XMLDocument ? Thanks, Joh Duncan

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

Hi, R-Helpers, I would like to ask about multiple graphs in one figure. I tried to execute the following codes. xlim - c(1,100) ylim - c(1,4) plot(NA, xlim=xlim, ylim=ylim) x - c(1:100) for(j in seq(1,10,by=1)) { y - j*x^2+log(j) lines(x, y) } In the above codes, I had to

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

On Tue, Aug 31, 2010 at 03:26:12AM -0700, Wonsang You wrote: In the above codes, I had to arbitrarily set up the coordinate range of the figure in advance before calculating the values y. (seexlim and ylim) In results, the figure did not contain all data since most of data were outside the

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

On 08/31/2010 08:26 PM, Wonsang You wrote: Hi, R-Helpers, I would like to ask about multiple graphs in one figure. I tried to execute the following codes. xlim- c(1,100) ylim- c(1,4) plot(NA, xlim=xlim, ylim=ylim) x- c(1:100) for(j in seq(1,10,by=1)) { y- j*x^2+log(j)

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

Hi, It's easy with ggplot2, library(ggplot2) ## create an empty plot p - ggplot(map=aes(x,y)) ## create a dummy list of data.frames with different ranges d - replicate(4, data.frame(x=sample(1:10,1)+rnorm(10), y=sample(1:10,1)+rnorm(10)),

[R] simultaneous estimation

Hi folks, Not sure what this sort of estimation is called. I have a 2-column time-series x(i,t) [with (i=1,2; t=1,...T)], and I want to do the following 'simultaneous' regressions: x(1,t) = (d - 1)(x(1, t-1) - mu(1)) x(2,t) = (d - 1)(x(2, t-1) - mu(2)) And I want to determine the coefficients

Re: [R] simultaneous estimation

On Tue, Aug 31, 2010 at 6:58 AM, murali.me...@avivainvestors.com wrote: Hi folks, Not sure what this sort of estimation is called. I have a 2-column time-series x(i,t) [with (i=1,2; t=1,...T)], and I want to do the following 'simultaneous' regressions: x(1,t) = (d - 1)(x(1, t-1) - mu(1))

Re: [R] simultaneous estimation

On 31/08/2010 6:58 AM, murali.me...@avivainvestors.com wrote: Hi folks, Not sure what this sort of estimation is called. I have a 2-column time-series x(i,t) [with (i=1,2; t=1,...T)], and I want to do the following 'simultaneous' regressions: x(1,t) = (d - 1)(x(1, t-1) - mu(1)) x(2,t) = (d -

Re: [R] About plot graphs

Hi Greg, Thanks for your advice. Performed following test:- attach(women) library(lattice) x11(height=6, width=18) par(mfrow=c(1,3)) plot(weight,height) attach(singer) The following object(s) are masked from women : height plot(weight,height) Error in xy.coords(x, y,

Re: [R] Regarding naive baysian classifier in R

Vishnampettai akron_aadhithya at yahoo.com writes: Hi, I have a small doubt regarding naive Bayes. I am able to classify the data's properly but i am just stuck up with how to get the probability values for naive bayes. In the case of SVM we have attr function that helps in displaying the

[R] Where can I download/install grDevices

Hi folks, Debian 504 64 bit grDevices is NOT on Debian repo $ apt-cache search r-cran | grep grDevices no printout On R $ sudo R install.packages(grDevices, dependencies=TRUE)) Warning in install.packages(grDevices) : argument 'lib' is missing: using '/usr/local/lib/R/site-library' ---

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

On Tue, Aug 31, 2010 at 04:50:39AM -0700, Stephen Liu wrote: Hi folks, Debian 504 64 bit grDevices is NOT on Debian repo $ apt-cache search r-cran | grep grDevices no printout On R $ sudo R install.packages(grDevices, dependencies=TRUE)) Warning in install.packages(grDevices) :

Re: [R] Finding pairs

Something like this may get you started library(igraph) df- data.frame(A=0:7, B=c(rep(NA,3), 0, 1, 3, 7, 2)) # Assuming this is your data b- as.matrix(df) print(b) # A B #[1,] 0 NA #[2,] 1 NA #[3,] 2 NA #[4,] 3 0 #[5,] 4 1 #[6,] 5 3 #[7,] 6 7 #[8,] 7 2 z- graph.edgelist(na.omit(b))

Re: [R] getNodeSet - what am I doing wrong?

Johannes Graumann wrote: Thanks! but: library(XML) xmlDoc - xmlTreeParse(http://www.unimod.org/xml/unimod_tables.xml;) You need to xmlParse() or xmlTreeParse(url, useInternalNodes = TRUE) (which are equivalent) in order to be able to use getNodeSet(). The error you are getting is because

Re: [R] getNodeSet - what am I doing wrong?

Sorry about that - got dropped from my attempts yesterday (see the first example below, that has the useInternalNodes=TRUE) ... Thanks again, Joh Duncan Temple Lang wrote: Johannes Graumann wrote: Thanks! but: library(XML) xmlDoc -

Re: [R] R crashes when communicating with JAGS

Thanks for your answer, Ben. Indeed I tried to use the newest R version, but just to be sure I tried it again just now, with the same result, i.e. R crashes. Unfortunately I cannot even look at the system configurations as this is a university computer and they (stupidly) don't give us

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

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 message appears while I put the code to run for the second WinBUGS. The error message: Error in file(file, ifelse(append, a, w)) : cannot open

[R] [회신] R-help Digest, Vol 90, Issue 35

9쩔첫 1�€횕�€횣쨌횓 횁짝횁횜쨈챘횉횖짹쨀쩔징쩌 짹횢쨔짬횉횛쨈횕쨈횢. 쩐횛�€쨍쨌횓 쨍횧�€횕�€쨘 횁짚쩌짰짹횢 sukgeun.j...@gmail.com �€쨍쨌횓 쨘쨍쨀쨩횁횜쩌쩌쩔채. I no longer work at NFRDI. Please send your mail to me at sukgeun.j...@gmail.com. __ R-help@r-project.org

Re: [R] R crashes when communicating with JAGS

Florian Weiler wrote: Thanks for your answer, Ben. Indeed I tried to use the newest R version, but just to be sure I tried it again just now, with the same result, i.e. R crashes. Unfortunately I cannot even look at the system configurations as this is a university computer and they (stupidly)

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

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.org lenny/main Packages 100 /var/lib/dpkg/status $ dpkg -S

Re: [R] About plot graphs

On Tue, 2010-08-31 at 04:42 -0700, Stephen Liu wrote: Hi Greg, Thanks for your advice. Performed following test:- attach(women) library(lattice) x11(height=6, width=18) par(mfrow=c(1,3)) plot(weight,height) attach(singer) The following object(s) are masked from women

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

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

[R] Munich R Course - Parallel Computing with R

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,

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

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

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

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

[R] Cohen's f

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.:

[R] pasting together 2 character arrays

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.trt1

[R] ggplot sub- and superscript for axis name

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

[R] 3D tomography data

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

Re: [R] pasting together 2 character arrays

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 ... Array3

Re: [R] pasting together 2 character arrays

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='.')

Re: [R] pasting together 2 character arrays

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) j.m...@uea.ac.ukwrote: If possible I would like to combine two different character arrays in combinations Array1 -

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

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] Cohen's f

Iasonas Lamprianou lamprianou at 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:

Re: [R] simultaneous estimation

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

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

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

On Aug 31, 2010, at 10:35 AM, murali.me...@avivainvestors.com murali.me...@avivainvestors.com 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

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

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

[R] wavelet parameters with sowas

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

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

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] solved Re: online documentation only?

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:

[R] Unable to run dev.new

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

Re: [R] simultaneous estimation

On 31/08/2010 11:00 AM, David Winsemius wrote: On Aug 31, 2010, at 10:35 AM, murali.me...@avivainvestors.com murali.me...@avivainvestors.com 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

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

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 PM,

Re: [R] Unable to run dev.new

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

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

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

[R] anova and lm results differ

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] Unable to run dev.new

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 murdoch.dun...@gmail.com To: Stephen Liu sati...@yahoo.com Cc:

Re: [R] simultaneous estimation

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

[R] specify axis range

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] specify axis range

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 m.lathour...@imperial.ac.uk wrote: Dear all, I want to make a plot where I can specify the range of y axis. I

Re: [R] Unable to run dev.new

On Tue, Aug 31, 2010 at 9:10 AM, Stephen Liu sati...@yahoo.com 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

Re: [R] anova and lm results differ

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

[R] TukeyHSD for binomial distribution

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

Re: [R] Remove Loading required package message

On 27/08/10 10:19, Barry Rowlingson wrote: On Fri, Aug 27, 2010 at 9:09 AM, Sébastien Moretti sebastien.more...@unil.ch wrote: [...] require it quietly: require(sp,quietly=TRUE) Doesn't work for me: require(Matrix, quietly=TRUE) Loading required package: lattice

Re: [R] simultaneous estimation

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

Re: [R] Remove Loading required package message

On Tue, Aug 31, 2010 at 5:21 PM, Allan Engelhardt all...@cybaea.com wrote: On 27/08/10 10:19, Barry Rowlingson wrote: On Fri, Aug 27, 2010 at 9:09 AM, Sébastien Moretti sebastien.more...@unil.ch  wrote:   [...]  require it quietly: require(sp,quietly=TRUE) Doesn't work for me:

Re: [R] anova and lm results differ

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

[R] help - read SAS into R

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

[R] odfWeave

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

Re: [R] help - read SAS into R

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

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

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

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

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)) :

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

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.org

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

Thanks! That does exactly what I needed. On Tue, Aug 31, 2010 at 12:07 PM, David Winsemius dwinsem...@comcast.netwrote: 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

Re: [R] specify axis range

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.

Re: [R] help - read SAS into R

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] help - read SAS into R

-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 ? Thank

[R] Detecting Growth Trends

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

Re: [R] odfWeave

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 kub...@virginia.edu wrote:

[R] pairs with same xlim and ylim scale

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, they

Re: [R] anova and lm results differ

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.:

[R] any statement equals to 'goto'?

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 -

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

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

[R] summary in Hmisc and Latex

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

Re: [R] Detecting Growth Trends

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] vglm

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] ILFTableliabLimit AnnAggLimit DedAmt Loss TIL

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

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

Re: [R] About plot graphs

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

Re: [R] summary in Hmisc and Latex

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 due to

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

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

Re: [R] odfWeave

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 Te, 8/31/10, John Kane jrkrid...@yahoo.ca wrote: From: John Kane

Re: [R] help - read SAS into R

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 __

Re: [R] vglm

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] ILFTableliabLimit AnnAggLimit

Re: [R] vglm

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:20 PM,

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

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] any statement equals to 'goto'?

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. __

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

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.:

Re: [R] lattice help required

...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 there

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

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 (to

Re: [R] Detecting Growth Trends

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

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

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

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

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 gunter.ber...@gene.com 08/31/2010 04:41 PM To

[R] Retrieving predictions after fitting a model

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

  1   2   >