Re: [R] How to detect if R is running on Mac OS X?

Philippe Grosjean wrote: This question is probably trivial, but I don't find the answer. I have code that is different for Windows, Unix/Linux and Mac OSX. The man page of .Platform tells that .Platform$OS.type is the right way to test for it... but it also tels that it returns either

[R] Kernel density output

Hi, I'm using the density() command for a given vector x and I would like to know how to get the estimated value of the density for each element of the vector x instead of values corresponding to points from a grid. Thanks Florent Bresson __

Re: [R] Kernel density output

Florent Bresson schrieb: Hi, I'm using the density() command for a given vector x and I would like to know how to get the estimated value of the density for each element of the vector x instead of values corresponding to points from a grid. Maybe not the best/most efficient way to do

Re: [R] How to detect if R is running on Mac OS X?

I think there are two different questions here: 1) Is R running under darwin, the underlying OS of MacOS X? You can test that by Sys.info or R.version$platform. It seems unlikely that you would need this. We do need to distinguish darwin at C level, as it is a very unusual 'unix', a

Re: [R] Kernel density output

You can use approx() to interpolate from the grid. On Mon, 15 Jan 2007, Florent Bresson wrote: I'm using the density() command for a given vector x and I would like to know how to get the estimated value of the density for each element of the vector x instead of values corresponding to

[R] Problem with pdf, png, jpeg devices and files named CON on Window s

Hi I cannot seem to create any files that have the name CON before the file extension, i.e. all of the following fail: pdf(CON.pdf) Error in pdf(CON.pdf) : unable to start device pdf jpeg('CON.jpeg') Error in jpeg(CON.jpeg) : unable to start device devWindows png('CON.png') Error in

Re: [R] Problem with pdf, png, jpeg devices and files n amed CON on Window s

Looks like it's a windows' problem. You also cannot create a file with the name con and any extension from the command line. Try echo something con.txt __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do

Re: [R] Problem with pdf, png, jpeg devices and files named CON on Window s

It's not R problem but specification of windows file system. see MSDN: http://msdn2.microsoft.com/en-us/library/aa365247.aspx On 1/15/07, Brandt, T. (Tobias) [EMAIL PROTECTED] wrote: Hi I cannot seem to create any files that have the name CON before the file extension, i.e. all of the

[R] Conflict in .Rprofile documentation FAQ vs. Help?

Hi folks, I note that in the general FAQ's we have 7.25 Why did my .Rprofile stop working when I updated R? Did you read the NEWS file? For functions that are not in the base package you need to specify the correct package namespace, since the code will be run before the packages are loaded.

Re: [R] ks.test not working?

cannot compute correct p-values with ties in: ks.test(x, pgev, fit$mle[1], fit$mle[2], fit$mle[3]) You may want to use the ks.boot function in the Matching package which implements a bootstrap ks-test which provides consistent pvalues (achieved significance levels) when there are ties.

[R] Matrex R adapter released

Matrex (http://matrex.sourceforge.net/) is a spreadsheet equivalent application that calculates formulas with matrices parameters, not cells. Matrex has charts and presentations (spreadsheets) and is stable and multithreaded. The 12th of january we released MatrexR, that allows to have Matrex as

Re: [R] batch job GLM calculations

Hi, Someone help-me on this task some weeks ago, and it is working find to me. I use something like: - findmodels - function(modeltype = lm, dataset, pattern) { ls - ls(.GlobalEnv, pattern=pattern) mods - ls[sapply(ls, function(x) inherits(get(x), modeltype))] if

Re: [R] changes in the structure of mer objects?

Hi Joris, I suspect you somehow load an older version lme4 or Matrix than you think you are loading. Or then you have an lmer() function or a class definition {from a saved workspace } in your work space. example(lmer) *must* run correctly for the 'lme4' package to get onto CRAN at all,

Re: [R] Controlling size of boxplot when it is added in a plot

Charilaos == Charilaos Skiadas [EMAIL PROTECTED] on Sat, 13 Jan 2007 23:12:00 -0500 writes: Charilaos I am trying to add a boxplot to the bottom of a Charilaos histogram, right between the histogram bars and Charilaos the x axis. Here is the code I am using at the Charilaos

Re: [R] graph results of logistic regression

The code of the function plot.logi.hist is defined in appendix A of the article you referred to. Cheers, Thierry ir. Thierry Onkelinx Instituut voor natuur- en bosonderzoek / Reseach Institute for Nature and Forest

Re: [R] changes in the structure of mer objects?

Thanks Martin, R --vanilla did the trick. By the way, is there a way to check which version of a lme4 or Matrix I am using in a certain instance of R? We are running multiple versions of R on the same server until we are sure all our operational code is behaving well under a new version of R or

[R] Advanced course R programming and Bioconductor in Cambridge UK 30.3.+1.4.

Dear R users developers, Seth Falcon and Martin Morgan are teaching Advanced R Programming and Bioconductor 30 March - 1 April 2007 at the Wellcome Trust Genome Campus, Hinxton, Cambridge, UK. This two-day course focuses on programming skills required to develop software for statistical

[R] How to format R code in LaTex documents

Hi, I am planning on putting some R script in an appendix of a LaTex document. Can anyone recommend me a way of how to format it? Is there a way to keep all line breaks without having to insert \\ in every single line? Thank you! Benjamin __

Re: [R] How to format R code in LaTex documents

Benjamin Dickgiesser wrote: Hi, I am planning on putting some R script in an appendix of a LaTex document. Can anyone recommend me a way of how to format it? Is there a way to keep all line breaks without having to insert \\ in every single line? Thank you! Benjamin Here's one way and

Re: [R] Problem with pdf, png, jpeg devices and files named CON on Window s

This is a restriction of the ISO C interface to the Windows API. E.g. 'check' says ## Furthermore, Uwe Ligges says that Windows still does not allow ## the following DOS device names (by themselves or with possible ## extensions): ## ## NameFunction ##

[R] Re : How to format R code in LaTex documents

I use Tinn-R to insert code in a latex file I use the copy formatted for TEX. But the verbatin environnement in Latex produce what you want. Justin BEM Elève Ingénieur Statisticien Economiste BP 294 Yaoundé. Tél (00237)9597295.

Re: [R] changes in the structure of mer objects?

HaroldD == Doran, Harold [EMAIL PROTECTED] on Mon, 15 Jan 2007 09:18:12 -0500 writes: HaroldD help(package='lme4') will tell you yes, or with more input but less output packageDescription(lme4)$Version Further note sessionInfo() which gives you all versions of all

Re: [R] How to format R code in LaTex documents

On Jan 15, 2007, at 8:42 AM, Benjamin Dickgiesser wrote: Hi, I am planning on putting some R script in an appendix of a LaTex document. Can anyone recommend me a way of how to format it? Is there a way to keep all line breaks without having to insert \\ in every single line? I think the

[R] add points to a 1-weibull graph

Hi, I made a survival analysis using survreg function. Its OK. But, I need to show in the graph 1-weibull curves. Ploting weibull curve a add the estimated points following these steps: fit - survfit(Surv(Time,Censor)~Group) points(fit[1]) points(fit[2]) points(fit[3]) I have 3 groups. But

Re: [R] R2WinBugs and Compare DIC versus BIC or AIC

Thanks for your answer ! But If I have understood correctly, the procedure is correct and this only means that I have already generated the initial values From R. I define a function to generate the initial values as: inits. - function(){list(sdesp=runif(1),sdhet=runif(1),beta0=runif(1),

[R] options(digits) and print.default()

Hello everyone, I use latex() (Hmisc) for report generation and thus have been affected by the problem with rounding decimals described, for example, in this post: http://thread.gmane.org/gmane.comp.lang.r.general/73287/focus=73287 In short, numbers often are printed with 15 or so decimals even

Re: [R] R2WinBugs and Compare DIC versus BIC or AIC

Dae-Jin Lee wrote: Thanks for your answer ! But If I have understood correctly, the procedure is correct and this only means that I have already generated the initial values From R. I define a function to generate the initial values as: inits. -

[R] model simplification in repeated anova

Hi, I've done a repeated ANOVA looking as follows: model-aov(y~a*b*c+Error(d)) Now I want to get rid of some or all of the interactions. But neither AIC nor anova work. What is the appropriate function to test my model? __ R-help@stat.math.ethz.ch

Re: [R] wafer map drawing

Jim and Carlos, Thanks for the leads. Best Regards, Sam From: jim holtman [mailto:[EMAIL PROTECTED] Sent: Friday, January 12, 2007 7:02 PM To: Walker, Sam Cc: r-help Subject: Re: [R] wafer map drawing There is some information in Programming with

Re: [R] model simplification in repeated anova

summary(npk.aovE) summary(np.k.aovE) ## The block stratum has the same total in both summaries. ## In this example, the two suppressed interactions are last and ## both are not significant. ## In order to use anova(model.1, model.2), it is necessary to ## rewrite the SAME models without the Error

Re: [R] changes in the structure of mer objects?

help(package='lme4') will tell you -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of [EMAIL PROTECTED] Sent: Monday, January 15, 2007 8:37 AM To: Martin Maechler Cc: [EMAIL PROTECTED]; r-help@stat.math.ethz.ch; [EMAIL PROTECTED] Subject: Re:

Re: [R] options(digits) and print.default()

On Mon, 15 Jan 2007, Ulrich Keller wrote: Hello everyone, I use latex() (Hmisc) for report generation and thus have been affected by the problem with rounding decimals described, for example, in this post: http://thread.gmane.org/gmane.comp.lang.r.general/73287/focus=73287 In short,

Re: [R] Conflict in .Rprofile documentation FAQ vs. Help?

On Mon, 15 Jan 2007, MHH Stevens wrote: Hi folks, I note that in the general FAQ's we have 7.25 Why did my .Rprofile stop working when I updated R? Did you read the NEWS file? For functions that are not in the base package you need to specify the correct package namespace, since the code

Re: [R] model simplification in repeated anova

## My previous response included only half the file. ## This is based on ?aov ## From Venables and Ripley (2002) p.165. data(npk, package=MASS) ## as a test, not particularly sensible statistically npk.aovE - aov(yield ~ N*P*K + Error(block), npk) np.k.aovE - aov(yield ~ N*P+K +

Re: [R] listing all functions in R

Hello, First, regarding GeSHi syntax highlighting for R, I have done one for the R Wiki (plus the R function that generates the list of keywords automatically). I will attach it to a second email send privately to you, since the mailing list do not accept attachments. For the problem of keeping

[R] feedback on writing r extensions

Hi, To whom should I send feedback/edits concerning the Writing R Extensions web page? http://cran.r-project.org/doc/manuals/R-exts.html Scott __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read

Re: [R] Conflict in .Rprofile documentation FAQ vs. Help?

[Brian Ripley] No one actually said it was a *working* example [...] Do you mean that, whenever we see something presented as an example within or around the R system, we should not take it as dependable unless it is explicitly said to be working? (and it is enclosed in \dontrun{}) Within

[R] spatial correlaton in lme and gls

Hello, I have multiple observations (replicates) per station...and the stations are distributed in space. I've applied models using lme and gls from the nlme library. I've messed around with compound symmetry to capture the correlation within stations. My question is.I have x and y

Re: [R] How to format R code in LaTex documents

Consider using the fancyvrb package if you need additional customization. On 1/16/07, Benjamin Dickgiesser [EMAIL PROTECTED] wrote: Hi, I am planning on putting some R script in an appendix of a LaTex document. Can anyone recommend me a way of how to format it? Is there a way to keep all

[R] How can I generate line segments between two sets of points on the same graph?

I am drawing two plots on the same graphics device window as follows: yLim-range(range(template),range(slope*subSeries) ) plot( template, type=p, pch=19,cex=1.75,lwd=2.5,lty=1, ann=F, axes=F,frame.plot=T,col=red, ylim=yLim,new=T ) title( main = myDotChart of original template and subSeries vs

[R] How do I generate line segments between two sets of points on the same graph

Just by re-examining my question I found the answer ! segments(3, template[3], 3,slope*subSeries[3], col= orange,lwd=2)I corrected the subSeries Points as they appear on the graph withslope*subSeries[3].thank youLloyd L [[alternative HTML version deleted]]

[R] ARIMA xreg and factors

I am using arima to develop a time series regression model, I am using arima b/c I have autocorrelated errors. Several of my independent variables are categorical and I have coded them as factors . When I run ARIMA I don't get any warning or error message, but I do not seem to get estimates for

Re: [R] ARIMA xreg and factors

Try this: set.seed(1) sim.ar - arima.sim(list(ar = c(0.4, 0.4)), n = 1000) z- gl(5, 1, 1000) zm - model.matrix(~z)[,-1] arima(sim.ar, order=c(2,0,0), xreg = zm) On 1/15/07, sj [EMAIL PROTECTED] wrote: I am using arima to develop a time series regression model, I am using arima b/c I have

[R] nonlinear regression: nls, gnls, gnm, other?

Hi all, I'm trying to fit a nonlinear (logistic-like) regression, and I'd like to get some recommendations for which package to use. The expression I want to fit is something like: y ~ A * exp(X * Beta1) / (1 + exp(-(x + X * Beta2 - xmid)/scal)) Basically, it's a logistic function, but I want

Re: [R] How to format R code in LaTex documents

Another alternative is to simply use LyX (www.lyx.org), which is a WYSIWYG (graphical) front-end for LaTeX. There's an code environment where you can paste everything right in, and it'll show in your PDF output exactly as it appears while editing. Also, not an answer to your question, but you

Re: [R] Within-subject factors in lme

Thanks for your suggestion. I guess the model you are fitting here has only a single random effect term, namely subject. If the effect of A depends on S, one needs to include an additional random effects term for the S:A interaction. It is not quite clear what you are after, but the 'if'

Re: [R] Within-subject factors in lme: complemented

The internet address to complement this is: http://wiki.r-project.org/rwiki/doku.php?id=guides:lmer-tests sorry for the double posting! -Original Message- From: Gygax Lorenz ART Sent: Tuesday, January 16, 2007 7:32 AM To: r-help@stat.math.ethz.ch Subject: RE: [R] Within-subject

Re: [R] R editor vs. Tinn-R

Hi On 12 Jan 2007 at 19:11, Martin Becker wrote: Date sent: Fri, 12 Jan 2007 19:11:42 +0100 From: Martin Becker [EMAIL PROTECTED] To: Farrel Buchinsky [EMAIL PROTECTED] Copies to: r-help@stat.math.ethz.ch Subject:Re:

Re: [R] ARIMA xreg and factors

From the help page xreg: Optionally, a vector or matrix of external regressors, which must have the same number of rows as 'x'. Note, not a factor. (It handles your factor like the integer vector it is internally.) If you want a design matrix you need to make one, something