Re: [R] Problem with bargraph.CI in Sciplot package

2009-04-11 Thread Dieter Menne
Metconnection simontbate at hotmail.co.uk writes: Hi there, I wonder if anyone can help me. I'm trying to use bargraph.CI in the Sciplot package when there is a missing combination of the factor levels. Unfortunately the standard errors on the plot do not appear to be correct. ...

Re: [R] Levene's test for homogeneity of variances (befor using ANOVA)

2009-04-11 Thread pleco
First you should download the car package from the r project website. Than load it by typing library(car). The actual code for levene's test is levene.test (y, group) Kes Knave wrote: Dear all I have tried to find this function in R, but don't find it by searching in the help

[R] Install odfWeave on Ubuntu 8.10 Intrepid Ibex

2009-04-11 Thread j.k
Hello togheter, I've installed R with this instruction: http://edtechdev.blogspot.com/2009/03/r-in-ubuntu-intrepid.html After that I tried to install odfWeave with this: install.packages(odfWeave) Then there came this warning: Warnung in install.packages(odfWeave) : Argument 'lib' missing:

Re: [R] Install odfWeave on Ubuntu 8.10 Intrepid Ibex

2009-04-11 Thread j.k
I've already found the solution myself: the package odfWeave depends on package XML to configure the package XML the libxml2-devel package in ubuntu is needed. You get that with typing in the console: sudo apt-get install libxml2-devel j.k wrote: Hello togheter, I've installed R with

[R] leave-one-out in R

2009-04-11 Thread Charles Serele
Hi Everyone, I am new in using R and I was wondering if anybody knows how to do a leave-one-out cross-validation in R. Thanks Charles Découvrez les photos les plus intéressantes du jour. [[alternative HTML version deleted]] __

Re: [R] Plot legend

2009-04-11 Thread Jim Lemon
Sherri Heck wrote: Dear All- I am trying to create a legend where the first 9 terms are represented as pch=16 in one column (which I am able to do); but, I would like the last two terms in that column to be 'lines (which are black and red - just as it looks like in the plot). I have read

Re: [R] Coloring X and Y axis

2009-04-11 Thread Jim Lemon
Gundala Viswanath wrote: Dear all, Is there a way to do it? The following code: pdf(test.pdf) plot(1,1,col=red) dev.off() Only colors the plot into red, but not x and y axis (inclusive the tick marks). Hi Gundala, One easy way is to roll your own axes: plot(1,1,col=red)

Re: [R] Coloring X and Y axis

2009-04-11 Thread Umesh Srinivasan
Hi Gundala, If you want your plotting axes to appear in colour, try par(fg = red) # will colour all foreground stuff in red, including the axes. Cheers, Umesh On Sat, Apr 11, 2009 at 3:51 PM, Jim Lemon j...@bitwrit.com.au wrote: Gundala Viswanath wrote: Dear all, Is there a way to do

[R] transfer function model of time series

2009-04-11 Thread Jim Lemon
Hi Rashimba, Mate, do you really think that we are going to do your examination for you? I haven't seen anything that rich since a student handed in a paper with the reference list photocopied from someone else's paper. Jim __ R-help@r-project.org

Re: [R] Problem with bargraph.CI in Sciplot package

2009-04-11 Thread Manuel Morales
On Fri, 2009-04-10 at 15:07 -0700, Metconnection wrote: Hi there, I wonder if anyone can help me. I'm trying to use bargraph.CI in the Sciplot package when there is a missing combination of the factor levels. Unfortunately the standard errors on the plot do not appear to be correct.

Re: [R] MLE for bimodal distribution

2009-04-11 Thread _nico_
Just wanted to thank everyone for their help, I think I mostly managed to solve my problem. -- View this message in context: http://www.nabble.com/MLE-for-bimodal-distribution-tp22954970p23000785.html Sent from the R help mailing list archive at Nabble.com.

Re: [R] leave-one-out in R

2009-04-11 Thread Uwe Ligges
See ?errorest in package ipred, for example. Uwe Ligges Charles Serele wrote: Hi Everyone, I am new in using R and I was wondering if anybody knows how to do a leave-one-out cross-validation in R. Thanks Charles Découvrez les photos les plus intéressantes du jour.

Re: [R] Cross-platforms solution to export R graphs

2009-04-11 Thread Philippe Grosjean
..°})) ) ) ) ) ) ( ( ( ( (Prof. Philippe Grosjean ) ) ) ) ) ( ( ( ( (Numerical Ecology of Aquatic Systems ) ) ) ) ) Mons-Hainaut University, Belgium ( ( ( ( ( .. Liviu

Re: [R] Cross-platforms solution to export R graphs

2009-04-11 Thread Philippe Grosjean
Liviu Andronic wrote: Hello, On Thu, Apr 9, 2009 at 3:04 PM, Philippe Grosjean phgrosj...@sciviews.org wrote: Cross-platforms solution to export R graphs There is playwith, and latticist, which seem cross-platform (binaries available for both MacWin). rattle uses latticist. Yes, right.

Re: [R] DierckxSpline fitting with different sets of y-values in one time

2009-04-11 Thread spencerg
Dear Elisabeth: I got the following error from your call to curfit.free.knot: splinefit-with(DF, curfit.free.knot(Age, Connectivity, k=1, g=1, prior=NULL)) Error in xyw.coords(x, y, w) : Need at least 4 distinct x values. Only 3 distinct values found (after combining observations

[R] (no subject)

2009-04-11 Thread Mary Winter
Hi, Is it possible to increase just the text size of the main title on my graph? If so, how?! Along the same lines..along the x axis at a lot of the labels overlap each other, is it possible to make this text smaller of write it on two lines so it doesn't take up as much space along

Re: [R] (no subject)

2009-04-11 Thread Steve henshaw
Hi, If you have a look at the help for par (type ?par at the prompt) and then have a read of the section on 'cex' that should point you in the right direction. You can pass the cex argument to the plot routines as well rather than changing them in par. Steve Mary Winter wrote: Hi,

Re: [R] (no subject)

2009-04-11 Thread Sarah Goslee
Sure, look at ?par - especially the various cex.* options. Sarah On Sat, Apr 11, 2009 at 10:44 AM, Mary Winter statsstud...@hotmail.com wrote: Hi, Is it possible to increase just the text size of the main title on my graph? If so, how?! Along the same lines..along the x axis at a

Re: [R] (no subject)

2009-04-11 Thread ronggui
?par and please note the ctex.main parameter. BTW, it would be better if you use an informative email subject. Best 2009/4/11 Mary Winter statsstud...@hotmail.com: Hi, Is it possible to increase just the text size of the main title on my graph? If so, how?! Along the same

Re: [R] (no subject)

2009-04-11 Thread Jorge Ivan Velez
Dear Mary, Try plot(1,1,main='Title here',cex.main=4) See ?par, especially the cex.* parameter(s). HTH, Jorge On Sat, Apr 11, 2009 at 10:44 AM, Mary Winter statsstud...@hotmail.comwrote: Hi, Is it possible to increase just the text size of the main title on my graph? If so, how?!

[R] Protocol question: formatting of posts

2009-04-11 Thread Carl Witthoft
Here's a random example of some code and the response R produced: sum(2 ^ (which(as.logical(rev(nn))) - 1)) [1] 74 --- Now, maybe I'm old and crochety, but why don't people put a handy linefeed (or CR) between each command line and each response line? I just find it

Re: [R] png with ggplot on windows xp

2009-04-11 Thread ONKELINX, Thierry
Dear Juliet, ggsave uses 300 dpi as default and the png device 72 dpi. The problem seems to be that everything is scaled properly except the pointsize. At least on WinXP. A work around is to increase the theme base_size from 12 to 50 (12 * 300 / 72 = 50). p - p + theme_gray(50) ggsave(file =

[R] Labeling points on plot on relative warp scores?

2009-04-11 Thread stephalope
Hi there, I am plotting relative warp scores (equivalent to pca scores) and I want to label (color code and shape) the points by group. I can't figure out how to do this beyond simple plotting. plot(RW1, RW2); Do I need to make vectors of each group and then plot them separately onto the same

Re: [R] split a character variable into several character variable by a character

2009-04-11 Thread Darren Norris
just an alternative try gsub ?gsub from Francisco's example: dat-read.table(clipboard, header=T)#Read from your email gsub(-.*,,as.character(dat$popcode))# gives the BCPy01 part of column popcode gsub(.*-,,as.character(dat$popcode)) # gives the 01 part of column popcode then to add these

Re: [R] Protocol question: formatting of posts

2009-04-11 Thread Marc Schwartz
On Apr 11, 2009, at 10:38 AM, Carl Witthoft wrote: Here's a random example of some code and the response R produced: sum(2 ^ (which(as.logical(rev(nn))) - 1)) [1] 74 --- Now, maybe I'm old and crochety, but why don't people put a handy linefeed (or CR) between each

[R] question related to fitting overdispersion count data using lmer quasipoisson

2009-04-11 Thread Sean Zhang
Dear R-helpers: I have a question related to fitting overdispersed count data using lmer. Basically, I simulate an overdispsed data set by adding an observation-level normal random shock into exp(+rnorm()). Then I fit a lmer quasipoisson model. The estimation results are very off (see model

[R] Sean / Re: question related to fitting overdispersion count data using lmer quasipoisson

2009-04-11 Thread Sean Zhang
Hey Buddy, Hope you have been doing well since last contact. If you have the answer to the following question, please let me know. If you have chance to travel up north. let me know. best, -Sean -- Forwarded message -- From: Sean Zhang seane...@gmail.com Date: Sat, Apr 11,

[R] Error in R CMD check 2.8.1

2009-04-11 Thread Michael Friendly
Env: R 2.8.1, Win Xp, Eclipse/StatET In a .Rd file, I have an example containing the lines: # calculate Y M, using polynomial contrasts trends - as.matrix(VocabGrowth) %*% poly(8:11, degree=3) colnames(trends)- c(Linear, Quad, Cubic) [At the risk of a long message, I'll append the complete

Re: [R] Labeling points on plot on relative warp scores?

2009-04-11 Thread Uwe Ligges
Specify corresponding vectors for the argument cex and pch. Best, Uwe Ligges stephalope wrote: Hi there, I am plotting relative warp scores (equivalent to pca scores) and I want to label (color code and shape) the points by group. I can't figure out how to do this beyond simple plotting.

Re: [R] transfer function model of time series

2009-04-11 Thread Rashmiba Zala
Please remove my R posting from the website as soon as possible. No more i need help. - Original Message - From: Rashmiba Zala rz...@uoguelph.ca To: r-help@R-project.org Sent: Saturday, April 11, 2009 8:27:37 AM (GMT-0800) Auto-Detected Subject: Re: transfer function model of time series

Re: [R] Error in R CMD check 2.8.1

2009-04-11 Thread Duncan Murdoch
On 11/04/2009 12:33 PM, Michael Friendly wrote: Env: R 2.8.1, Win Xp, Eclipse/StatET In a .Rd file, I have an example containing the lines: # calculate Y M, using polynomial contrasts trends - as.matrix(VocabGrowth) %*% poly(8:11, degree=3) colnames(trends)- c(Linear, Quad, Cubic) [At the

[R] Help with postscript (huge file size)

2009-04-11 Thread Talita Perciano
Dear users, I'm generating some images in R to put into a document that I'm producing using Latex. This document in Latex is following a predefined model, which does not accept compilation with pdflatex, so I have to compile with latex - dvi - pdf. Because of that, I have to generate the images

[R] Fatal Error (was Error in saveLog(currentLogFileName)

2009-04-11 Thread Prew, Paul
Hello, I hope someone can advise me on how to get R running on my Windows PC again. I uninstalled R2.8.1, and also the previous version 2.8 that were on my hard drive. I reinstalled 2.8.1, but when I started the program, the message below came up, and R shut down. Fatal error: unable to

[R] How to save independent R sessions in BATCH mode?

2009-04-11 Thread mcnda839
Hi Rers, I would like to use the BATCH mode to R automatically: open 2 R sessions (ex. data1.Rdata, data2.Rdata), produce some calculations (…), save 2 sessions (sol1.Rdata, sol2.Rdata), My problem when using R CMD BATCH master.R is that when I use master.R as a master file, with:

Re: [R] Cross-platforms solution to export R graphs

2009-04-11 Thread Liviu Andronic
On Sat, Apr 11, 2009 at 3:20 PM, Philippe Grosjean phgrosj...@sciviews.org wrote: format (PDF) to another one (SVG). In Inkscape, you use File - Open... for the first step, and File - Save as... for the second. Since it is a vector format, your graph should not look pixelised. Yes, this is

Re: [R] How to save independent R sessions in BATCH mode?

2009-04-11 Thread Sarah Goslee
You need to delete all the files initially loaded in data1.Rdata using rm() eg rm(list=ls(pat=.*)) Sarah On Sat, Apr 11, 2009 at 1:35 PM, mcnda...@mncn.csic.es wrote: Hi Rers, I would like to use the BATCH mode to R automatically: open 2 R sessions (ex. data1.Rdata, data2.Rdata), produce

Re: [R] Help with postscript (huge file size)

2009-04-11 Thread Ben Bolker
Talita Perciano wrote: Dear users, I'm generating some images in R to put into a document that I'm producing using Latex. This document in Latex is following a predefined model, which does not accept compilation with pdflatex, so I have to compile with latex - dvi - pdf. Because of

Re: [R] Help with postscript (huge file size)

2009-04-11 Thread Talita Perciano
Thank you for the answer. Just to clear things out, I'm generating plots of rgb images. Talita 2009/4/11 Ben Bolker bol...@ufl.edu Talita Perciano wrote: Dear users, I'm generating some images in R to put into a document that I'm producing using Latex. This document in Latex is

Re: [R] Help with postscript (huge file size)

2009-04-11 Thread Ben Bolker
Do you mean you're importing jpegs or other bitmaps into R and writing them out (possibly with annotation etc.) as PostScript? Can you give a small example of some sort? It would help for giving advice. Talita Perciano wrote: Thank you for the answer. Just to clear things out, I'm

Re: [R] Fatal Error (was Error in saveLog(currentLogFileName)

2009-04-11 Thread Uwe Ligges
Prew, Paul wrote: Hello, I hope someone can advise me on how to get R running on my Windows PC again. I uninstalled R2.8.1, and also the previous version 2.8 that were on my hard drive. I reinstalled 2.8.1, but when I started the program, the message below came up, and R shut down.

Re: [R] Help with postscript (huge file size)

2009-04-11 Thread Talita Perciano
Ok... I'm using the rimage package to manipulate an image. So, the image I have in R is of the type imagematrix, which is a matrix with the pixel values of the R, G anf B bands. What I'm doing is applying some operation (like laplacian filter for example) and plotting the result as an image:

Re: [R] Cross-platforms solution to export R graphs

2009-04-11 Thread Philippe Grosjean
Liviu Andronic wrote: On Sat, Apr 11, 2009 at 3:20 PM, Philippe Grosjean phgrosj...@sciviews.org wrote: format (PDF) to another one (SVG). In Inkscape, you use File - Open... for the first step, and File - Save as... for the second. Since it is a vector format, your graph should not look

Re: [R] Curve fitting,FDA for biological data

2009-04-11 Thread spencerg
Dear Dr Gkikopoulos: 1. Have you looked at bioconductor.org? They have substantive extensions to R specifically for genomic data, which I assume would include chromosome. 2. To identify periodicities at different timescales, I agree with Stephen that spectrum would likely

[R] who happenly read these two paper Mohsen Pourahmadi (biometrika1999, 2000)

2009-04-11 Thread Manli Yan
http://biomet.oxfordjournals.org/cgi/reprint/86/3/677 biometrika1999 http://biomet.oxfordjournals.org/cgi/reprint/94/4/1006 biometrika2000 Hi All: I just want to try some luck. I am currenly working on my project,one part of my project is to reanalysis the kenward cattle data by using the

[R] _Possible_ work-round for normalizePath error (was Re: [Rd] Package (PR#13475))

2009-04-11 Thread S Ellison
Uwe Ligges lig...@statistik.tu-dortmund.de 04/10/09 4:21 PM Well, you need to ask Symantec to fix Norton, hence this is the wrong address. That would not help other R users who were looking in the archives for help with what looked like a mysterious normalizePath error, would it? As I said,

Re: [R] leave-one-out in R

2009-04-11 Thread Max Kuhn
You can use the caret or ipred packages. caret has it setup for specific models and can be used with parallel processing. ipred can do loo for any model since you manually code the model. Max On Apr 11, 2009, at 3:03 AM, Charles Serele chser...@yahoo.ca wrote: Hi Everyone, I am new in

[R] looking for one-liner for strsplit and regex

2009-04-11 Thread Juliet Hannah
Hi, I have a line such as: myline - 0.10 1.5000 0.6000 538 0.369404 and I would like to put the numbers into a vector. Some combination of tabs and spaces occur between the numbers. I tried: try1 - strsplit(myline,[[:blank:]]+) try1 [[1]] [1]

[R] predict.lda posteriors equal to 1

2009-04-11 Thread gcam032
Hi all, I was wondering if anyone was able to explain to me why when I run predict.lda() on some new data I get lots of posterior probabilities = 1. This seems markedly unrealistic. What do the posteriors mean in the context of the predict.lda() function? Thanks Gareth -- View this message

[R] predict.lda posteriors equal to 1

2009-04-11 Thread gcam032
Hi all, I was wondering if anyone was able to explain to me why when I run predict.lda() on some new data I get lots of posterior probabilities = 1. This seems markedly unrealistic. What do the posteriors mean in the context of the predict.lda() function? Thanks Gareth -- View this message

Re: [R] predict.lda posteriors equal to 1

2009-04-11 Thread Bill.Venables
If you try library(MASS) example(lda) predict(z)$posterior You will see there are very many posterior probabilities there, too, which differe from 1 by round-off error at most. So your data are not unique. I don't get why you don't get it. Bill Venables

Re: [R] looking for one-liner for strsplit and regex

2009-04-11 Thread David Winsemius
What you have now is a list with one element which happens to be a character vector. Try: vec - scan(textConnection(myline)) Read 5 items vec [1] 0.10 1.50 0.60 538.00 0.369404 -- David Winsemius On Apr 11, 2009, at 8:23 PM, Juliet Hannah wrote: Hi, I have a line

Re: [R] Fatal Error (was Error in saveLog(currentLogFileName)

2009-04-11 Thread Prew, Paul
Dear Uwe and David, Thank you, deleting the .RData file did the trick. It's fantastic that you are so readily willing to help. Paul From: Uwe Ligges [mailto:lig...@statistik.tu-dortmund.de] Sent: Sat 4/11/2009 3:27 PM To: Prew, Paul Cc: r-help@r-project.org

Re: [R] predict.lda posteriors equal to 1

2009-04-11 Thread Gareth Campbell
Thanks Bill, Yes in my data too I get some different posteriors also. I guess what I was getting at was that when I'm dealing with multi-elemental data with 16 variables a probability of 1 is a bit more certain an assignment than I was expecting. I'm a student, so perhaps I don't have quite the