Re: [R] boxplot via plot command

2008-12-02 Thread Ken Knoblauch
Antje niederlein-rstat at yahoo.de writes: I've just discovered that the following code leads to boxplot (surprisingly to me). Can anybody explain to me why? Is this documented somewhere? I've never consider this option before. x - rnorm(300) l - c(rep(label1,100), rep(label2,50),

Re: [R] boxplot via plot command

2008-12-02 Thread Antje
Hi, thank you both for your response. I don't want to do anything like this - I just got some code like this from someone else and was wondering about the result. I would have used another approach to create a boxplot like this... Ciao, Antje [EMAIL PROTECTED] schrieb: hi: i'm not well

Re: [R] Projection Pursuit Classification

2008-12-02 Thread Uwe Ligges
ArunPrasad wrote: Hi, Can anyone help me to understand why I am getting the error message like below? Tree.result - PP.Tree(LDA,iris[train,5],iris[train,1:4]) These are *Warning* messages rather than error messages. These warnings come from coding problems in the functions you are

[R] Rcmd build problem: gcc.exe no input files (R.2.8.0 on Win XP)

2008-12-02 Thread Søren Højsgaard
Building a package using R.2.8.0 on Windows XP gives problems when compiling C code: Trying to run 'rcmd build' on R.2.8.0 on Windows XP gives ... -- Making package gRbase adding build stamp to DESCRIPTION making DLL ... gcc -std=gnu99 -Ic:/Programs/R/current/include

[R] merge character strings

2008-12-02 Thread mentor_
Hi, how can I do the following in R?: vec - c(a, b, c) magic. print(str) abc or even better: (a,b,c)? Cheers -- View this message in context: http://www.nabble.com/merge-character-strings-tp20790116p20790116.html Sent from the R help mailing list archive at Nabble.com.

Re: [R] merge character strings

2008-12-02 Thread Dimitris Rizopoulos
check at paste(), e.g., vec - c(a, b, c) paste(vec, collapse = ) I hope it helps. Best, Dimitris mentor_ wrote: Hi, how can I do the following in R?: vec - c(a, b, c) magic. print(str) abc Cheers -- Dimitris Rizopoulos Assistant Professor Department of Biostatistics Erasmus

Re: [R] Examples of advanced data visualization

2008-12-02 Thread Yihui Xie
Since SVG has a lot of elements, Firefox is still under development to support different elements in SVG; see: http://www.mozilla.org/projects/svg/status.html Currently animations won't work under Firefox 3.0.4, and all the rest can be viewed using Firefox 3.0.4. Regards, Yihui -- Yihui Xie

Re: [R] help with package Rwave

2008-12-02 Thread stephen sefick
I would suggest sowas (Douglas Maraun) or wmtsa (Percival and Walden). I don't know for sure but both of these probably rely on Rwave. I like both of the packages mentioned it just depends on what you want to do- wmtsa will do DWT and CWT + a lot more, and sowas is for testing hypotheses and

[R] Problem with tcl/tk and Rcmdr - urgent help required

2008-12-02 Thread beniamino
Hello everyone, i have been searching the whole day, trying ANY solution offered by the internet and before headbanging against the wall would like to try asking for your help. As many others, i am haveing problems loading the Rcmdr package with Macbook running Tiger 10.4.11. X11 version

Re: [R] Can't load Rgraphviz on windows-based R

2008-12-02 Thread ram basnet
Dear all, I am facing the same problem. I followed the all instructions given in help file and some steps of previous discussion (from google search) but I could not install package Rgraphiviz in windows.  I tried even the instructions given to install Bioconductor package then only able to

[R] command files

2008-12-02 Thread b g
Since I'm a SAS programmer, I'm used to creating command files in an editor for submission later. Is there a way to do this in R? I'd need to retain an ouput listing and a log to check for errors. _ Send e-mail faster without

Re: [R] Comparing output from linear regression to output from quasipoisson to determine the model that fits best.

2008-12-02 Thread Uwe Ligges
John Sorkin wrote: R 2.7 Windows XP I have two model that have been run using exactly the same data, both fit using glm(). One model is a linear regression (gaussian(link = identity)) the other a quasipoisson(link = log). I have log likelihoods from each model. Is there any way I can

[R] Limit on number of Rows when reading R dataset

2008-12-02 Thread Harsh
Hello, I am trying to read a dataset with 100,000 rows and around 365 columns into R, using read.table/read.csv. In Windows XP, with R 32 bit, I am able to read only 15266 rows and not more than that. I tried the same in R running in Ubuntu and it does the same and reads only 15266 rows. Using the

[R] fGarch and is.na()

2008-12-02 Thread Felix Lamp
Hello, after loading package fGarch the function is.na.data.frame() is behaving different: sessionInfo() R version 2.8.0 (2008-10-20) i386-pc-mingw32 locale: LC_COLLATE=English_United States.1252;LC_CTYPE=English_United States.1252;LC_MONETARY=English_United

Re: [R] command files

2008-12-02 Thread Blanchette, Marco
Try source('myFirstScript.R') Where myFirstScript.R as the following line x - rnorm(100) y - rnorm(100) plot(x,y) You could also use a editor like emacs with the ess-mode where one buffer can be your script with a live R session in a second buffer. Good luck On 12/2/08 7:21 AM, b g

Re: [R] Can't load Rgraphviz on windows-based R

2008-12-02 Thread Tobias Verbeke
L.S., I am facing the same problem. I followed the all instructions given in help file and some steps of previous discussion (from google search) but I could not install package Rgraphiviz in windows. I tried even the instructions given to install Bioconductor package then only able to

Re: [R] Regarding posting a package to R-forge (with one of the dependent packages not in CRAN)

2008-12-02 Thread Marcelo Perlin
Hi Uwe, Thanks for the reply. I already tried that a couple of weeks ago, with the exactly the same arguments you posted. So far no reply from the author. I going to try again. Regards, Marcelo. On Tue, Dec 2, 2008 at 9:42 AM, Uwe Ligges [EMAIL PROTECTED]wrote: Marcelo Perlin wrote: Hi

Re: [R] Limit on number of Rows when reading R dataset

2008-12-02 Thread Harsh
Thank you Uwe and Prof. Ripley. The problem was solved. The row in question indeed have garbage data, which possibly was truncating the number of lines read. I apologise for the oversight. Thank you once again. Regards Harsh Singhal Bangalore, India On Tue, Dec 2, 2008 at 2:50 PM, Prof Brian

Re: [R] boxplot via plot command

2008-12-02 Thread Ken Knoblauch
Hi, Antje niederlein-rstat at yahoo.de writes: Hi folks, I've just discovered that the following code leads to boxplot (surprisingly to me). Can anybody explain to me why? Is this documented somewhere? I've never consider this option before. x - rnorm(300) l - c(rep(label1,100),

Re: [R] Limit on number of Rows when reading R dataset

2008-12-02 Thread Prof Brian Ripley
Take a look at your dataset at around that row. Perhaps you have an unmatched quote? The limit on the number of rows of a data frame is far larger than 100,000 (2^31-1, but you will run out of address space on a 32-bit platform before that - see ?Memory-limits). On Tue, 2 Dec 2008, Harsh

Re: [R] Can't load Rgraphviz on windows-based R

2008-12-02 Thread Prof Brian Ripley
Rgraphviz is a Bioconductor package, so please ask on the apporpriate mailing list (see the posting guide). On Tue, 2 Dec 2008, Daren Tan wrote: Hi, I have problem loading Rgraphviz. Following the instructions specified by the README in Rgraphviz_1.20.3.tar.gz didn't help either. o. set

Re: [R] Problem with tcl/tk and Rcmdr - urgent help required

2008-12-02 Thread John Fox
Dear Ben, It seems clear that the problem is in loading the tcltk package. To confirm that, does library(tcltk) fail? I'm not very knowledgeable about Mac OS X. I do have a MacBook, with OS X 10.5, however, and encountered no problems installing and running the Rcmdr package. Have you seen the R

[R] Wiley Interdisciplinary Reviews

2008-12-02 Thread hadley wickham
They are obviously growing desperate - I have now been asked to write 5 articles! Hadley -- http://had.co.nz/ __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide

Re: [R] Using grep() to subset lines of text

2008-12-02 Thread Uwe Ligges
ppaarrkk wrote: I have two vectors, a and b. b is a text file. I want to find in b those elements of a which occur at the beginning of the line in b. I have the following code, but it only returns a value for the first value in a, but I want both. Any ideas please. a = c(2,3) b = NULL b[1]

Re: [R] How to solve following equation?

2008-12-02 Thread Hans W. Borchers
I need to solve a equation like this : a = b/(1+x) + c/(1+x)^2 + d/(1+x)^3 where a,b,c,d are known constant. Is there any R-way to do that? Multiplying this expression with (1+x)^3 leads to a polynomial equation. I would certainly recommend the 'PolynomF' package here: #

[R] Left-truncated regression

2008-12-02 Thread Zita
Hi. I am looking for a function for left-truncated data. I have one data set with 2 variables (Hours~Yrs_Ed). I already left-censored the data at 200 and left-truncated it at the same spot, so that I am able to make 2 estimations (one for censoring and one for truncation). I know how to make the

[R] Can't load Rgraphviz on windows-based R

2008-12-02 Thread Daren Tan
Hi, I have problem loading Rgraphviz. Following the instructions specified by the README in Rgraphviz_1.20.3.tar.gz didn't help either. o. set the following Windows environment variables accordingly (control panel - systems - Advanced - Environment Variables ): (a) create new user

[R] boxplot via plot command

2008-12-02 Thread Antje
Hi folks, I've just discovered that the following code leads to boxplot (surprisingly to me). Can anybody explain to me why? Is this documented somewhere? I've never consider this option before. x - rnorm(300) l - c(rep(label1,100), rep(label2,50), rep(label3,150)) df -

Re: [R] Regarding posting a package to R-forge (with one of the dependent packages not in CRAN)

2008-12-02 Thread Uwe Ligges
Marcelo Perlin wrote: Hi Guys, Recently I wrote a package for dealing with Markov Switching Regressions in R and it is included in the Rmetrics project. https://r-forge.r-project.org/projects/rmetrics/ Everything works fine when I use it in computer. But, the package depends on the use of

[R] Question to regular expressions

2008-12-02 Thread Antje
Hi there, I know, this question is not directly an R-help question but probably someone can give me a hint how to deal with the following problem. I have a vector with file/folder names and want to filter for all entries which have 6 numbers in a row and nothing else. folders - c(folder1,

Re: [R] optimization problem

2008-12-02 Thread Mike Prager
Hans W. Borchers [EMAIL PROTECTED] wrote: Why not use one of the global optimizers in R, for instance 'DEoptim', and then apply optim() to find the last six decimals? I am relatively sure that the Differential Evolution operator has a better chance to come near a global optimum than a loop

Re: [R] Limit on number of Rows when reading R dataset

2008-12-02 Thread Uwe Ligges
Harsh wrote: Hello, I am trying to read a dataset with 100,000 rows and around 365 columns into R, using read.table/read.csv. In Windows XP, with R 32 bit, I am able to read only 15266 rows and not more than that. I tried the same in R running in Ubuntu and it does the same and reads only

[R] Problem with the Rmpi package

2008-12-02 Thread Blanchette, Marco
Dear all, I just started to use the snow package to send multiple jobs on our cluster using MPI and the Rpmi package as the communication method. However, the Rmpi package have been behaving strangely. When I try to detach the Rmpi package I get the following error message: library(Rmpi)

Re: [R] Question to regular expressions

2008-12-02 Thread Gabor Grothendieck
Try this: folders - c(folder1, f2, F234562, 12345678, 234567, 912345, 333) grep(^[0-9]{6}$, folders, value = TRUE) [1] 234567 912345 On Tue, Dec 2, 2008 at 10:32 AM, Antje [EMAIL PROTECTED] wrote: Hi there, I know, this question is not directly an R-help question but probably someone

[R] legend idea for latticeExtra

2008-12-02 Thread baptiste auguie
Dear list, I've written a small utility function to add arbitrary legend(s) to a lattice graph (or a combination of them), much like the legend function of base graphics. I though perhaps it could be useful to someone else, or improved by suggestions. I understand this goes against the

Re: [R] Question to regular expressions

2008-12-02 Thread Antje
Hi Gabor, it works! Thank you very much! But I still don't understand the difference between [0-9] and [:digit:]... Ciao, Antje Gabor Grothendieck schrieb: Try this: folders - c(folder1, f2, F234562, 12345678, 234567, 912345, 333) grep(^[0-9]{6}$, folders, value = TRUE) [1] 234567

Re: [R] Question to regular expressions

2008-12-02 Thread Gabor Grothendieck
Try this: ^[[:digit:]]{6}$ On Tue, Dec 2, 2008 at 10:48 AM, Antje [EMAIL PROTECTED] wrote: Hi Gabor, it works! Thank you very much! But I still don't understand the difference between [0-9] and [:digit:]... Ciao, Antje Gabor Grothendieck schrieb: Try this: folders - c(folder1, f2,

Re: [R] exclude a vector value from another vector

2008-12-02 Thread Greg Snow
Others have shown ways to remove your sample from the population, but this may be doing things the long way. If you just want several samples from the same population without overlap between any of the samples, just take one sample of size equal to the sum of the individual sample sizes,

Re: [R] exclude a vector value from another vector

2008-12-02 Thread Henrique Dallazuanna
Try this: setdiff(pop, sample(pop, 2)) On Mon, Dec 1, 2008 at 5:16 PM, Hamid Hamid [EMAIL PROTECTED] wrote: Dear All, I am trying to build a program which will take repeated samples (w/o replacement) from a population of values. The interesting catch is that I would like the sample values

Re: [R] Question to regular expressions

2008-12-02 Thread hadley wickham
On Tue, Dec 2, 2008 at 9:48 AM, Antje [EMAIL PROTECTED] wrote: Hi Gabor, it works! Thank you very much! But I still don't understand the difference between [0-9] and [:digit:]... You might find this site helpful: http://regexp.resource.googlepages.com/analyzer.html Copy in your attempt and

[R] Bug in transform?

2008-12-02 Thread Vitalie Spinu
Dear useRs, Here is a weird behavior of transform function: mtcars1-matcars transform(mtcars1,t1=3,t2=4) Error in data.frame(`_data`, e[!matched]) : arguments imply differing number of rows: 32, 1 instead, this works: mtcars1$t1-0 transform(mtcars1,t1=3,t2=4) also works if applied

Re: [R] command files

2008-12-02 Thread Barry Rowlingson
2008/12/2 b g [EMAIL PROTECTED]: Since I'm a SAS programmer, I'm used to creating command files in an editor for submission later. Is there a way to do this in R? I'd need to retain an ouput listing and a log to check for errors. You probably want R CMD BATCH from a command-line. For

Re: [R] Question to regular expressions

2008-12-02 Thread Prof Brian Ripley
Antje wrote: Hi Gabor, it works! Thank you very much! But I still don't understand the difference between [0-9] and [:digit:]... If all else fails, read the help, here ?regex. Both [0-9] and [[:digit:]] are character classes of digits, but the first contains only arabic numerals. In some

Re: [R] Bug in transform?

2008-12-02 Thread Gabor Grothendieck
Try: cbind(mtcars, t1 = 3, t2 = 4) On Tue, Dec 2, 2008 at 11:14 AM, Vitalie Spinu [EMAIL PROTECTED] wrote: Dear useRs, Here is a weird behavior of transform function: mtcars1-matcars transform(mtcars1,t1=3,t2=4) Error in data.frame(`_data`, e[!matched]) : arguments imply differing

Re: [R] command files

2008-12-02 Thread Gabor Grothendieck
On Tue, Dec 2, 2008 at 11:18 AM, Barry Rowlingson [EMAIL PROTECTED] wrote: 2008/12/2 b g [EMAIL PROTECTED]: Since I'm a SAS programmer, I'm used to creating command files in an editor for submission later. Is there a way to do this in R? I'd need to retain an ouput listing and a log to

Re: [R] Bug in transform?

2008-12-02 Thread Prof Brian Ripley
As the help page says If some of the values are not vectors of the appropriate length, you deserve whatever you get! So you can use mtcars1 - mtcars mtcars1[c(t1, t2)] - cbind(rep(3,32), rep(4, 32)) or even mtcars1 - transform(mtcars, t1=rep(3, 32), t2=rep(4, 32)) Vitalie Spinu

[R] linear functional relationships with heteroscedastic non-Gaussian errors - any packages around?

2008-12-02 Thread Jarle Brinchmann
[apologies if this appears twice] Hi, I have a situation where I have a set of pairs of X Y variables for each of which I have a (fairly) well-defined PDF. The PDF(x_i) 's and PDF(y_i)'s are unfortunately often rather non-Gaussian although most of the time not multi-modal. For these data

Re: [R] Can't load Rgraphviz on windows-based R

2008-12-02 Thread Rory.WINSTON
You can get LoadLibrary failures if a dependency is not found (e.g. a required dll or shlib used by the dll that you are explicitly loading). So make sure that all necessary dependencies (e.g. graphviz libs) are available in the system path. Rory Rory Winston RBS Global Banking Markets

Re: [R] Bug in transform?

2008-12-02 Thread Vitalie Spinu
On Tue, 02 Dec 2008 17:37:44 +0100, Prof Brian Ripley [EMAIL PROTECTED] wrote: As the help page says If some of the values are not vectors of the appropriate length, you deserve whatever you get! So you can use mtcars1 - mtcars mtcars1[c(t1, t2)] - cbind(rep(3,32), rep(4, 32))

Re: [R] legend idea for latticeExtra

2008-12-02 Thread Deepayan Sarkar
On 12/2/08, baptiste auguie [EMAIL PROTECTED] wrote: Dear list, I've written a small utility function to add arbitrary legend(s) to a lattice graph (or a combination of them), much like the legend function of base graphics. I though perhaps it could be useful to someone else, or improved

Re: [R] linear functional relationships with heteroscedastic non-Gaussian errors - any packages around?

2008-12-02 Thread Prof Brian Ripley
I wonder if you are using this term in its correct technical sense. A linear functional relationship is V = a + bU X = U + e Y = V + f e and f are random errors (often but not necessarily independent) with distributions possibly depending on U and V respectively. and pairs from (X,Y) are

Re: [R] linear functional relationships with heteroscedastic non-Gaussian errors - any packages around?

2008-12-02 Thread Prof Brian Ripley
I wonder if you are using this term in its correct technical sense. A linear functional relationship is V = a + bU X = U + e Y = V + f e and f are random errors (often but not necessarily independent) with distributions possibly depending on U and V respectively. and pairs from (X,Y) are

Re: [R] factanal question

2008-12-02 Thread William Revelle
At 11:04 AM -0500 12/1/08, John Fox wrote: Dear Bill, Thanks for pointing out that this functionality is already in the psych package. Shouldn't factor.residuals() avoid this computation for oblique rotations? John, Good suggestion. I will add that in the next revision. Bill Regards,

[R] Asymmetric CIs

2008-12-02 Thread Sarah E. McCormick
Hi, I was wondering if there was some sort of package or function that calculated asymmetric confidence intervals for small proportions. I thought of both the epicalc and epitools package, but I am hoping to find something where you can just plug in a standard error and point estimate and it

Re: [R] Installation of RCurl Windows binary package from BioC extra repos (was Re: [BioC] Rcurl 0.8-1 update for bioconductor 2.7)

2008-12-02 Thread rdb20+
Roger, I had success installing RCurl per Patrick's instructions. I also got DAVIDQuery to install and run correctly but I had to run R CMD INSTALL DAVIDQuery from a cygwin terminal after 1) I replaced cygwin's make program v3.81 with the same version patched for windows from the mingw project.

Re: [R] Question to regular expressions

2008-12-02 Thread Antje
Thanks a lot again to all of you!!! Antje Antje schrieb: Hi Gabor, it works! Thank you very much! But I still don't understand the difference between [0-9] and [:digit:]... Ciao, Antje Gabor Grothendieck schrieb: Try this: folders - c(folder1, f2, F234562, 12345678, 234567, 912345,

[R] predicted probabilities after lmer

2008-12-02 Thread arpino
Dear R-users, I'm using lmer to fit two-level logistic models and I'm interested in predicted probabilities that I get in this way (using fitted): glm1 = lmer(XY$T1~X1 + X2 + X3 + (1|Cind), family=binomial) #estimation of a two-level logit model fit1=fitted(glm1) # I get the fitted

[R] question on lmer function

2008-12-02 Thread Richard Palmer
suppose something like probability(passing test) is driven by 1. fixed effects -- sex 2. district effects - district funding 3. school effects - neighborhood income, racial composition, % two parent families, ... 4. class effects - teacher quality measurement, 5. individual

Re: [R] Left-truncated regression

2008-12-02 Thread Tom La Bone
Look at the sand package, which is available at http://www.csm.ornl.gov/esh/statoed/ and the NADA package, which is available from CRAN. One or both may have items of interest. Tom Zita wrote: Hi. I am looking for a function for left-truncated data. I have one data set with 2

Re: [R] linear functional relationships with heteroscedastic non-Gaussian errors - any packages around?

2008-12-02 Thread Jarle Brinchmann
Thanks for the reply! On Tue, Dec 2, 2008 at 6:34 PM, Prof Brian Ripley [EMAIL PROTECTED] wrote: I wonder if you are using this term in its correct technical sense. A linear functional relationship is V = a + bU X = U + e Y = V + f e and f are random errors (often but not necessarily

Re: [R] Asymmetric CIs

2008-12-02 Thread Frank E Harrell Jr
Sarah E. McCormick wrote: Hi, I was wondering if there was some sort of package or function that calculated asymmetric confidence intervals for small proportions. I thought of both the epicalc and epitools package, but I am hoping to find something where you can just plug in a standard error

Re: [R] question on lmer function

2008-12-02 Thread Doran, Harold
Richard There is much more information that you need to provide before a thoughtful answer can be provided. Maybe you can describe the structure of your data, your outcome variable, etc. There is a vignette in the lmer package called 'Implementation' that will show you some methods for model

[R] question about the tisPlot function in package tis

2008-12-02 Thread Kyle . Matoba
List, I am using the 'tisPlot' function in Jeff Hallman's excellent tis package and was hoping that someone could spare me from having to dig into the code of his 'tisPlot' function. So far as I can tell, the preferred method of controlling the plotting of the x-axis is using the 'xTickFreq'

[R] match

2008-12-02 Thread threshold
Hi, I would like to check which rows of 'types.prev' matrix pop up in 'types', following R in-built procedure. I tried 'match' function but it works in case of the one dimensional vectors. Will appreciate any suggestions. best, robert types edateK [1,] 20060819 12.5 [2,]

Re: [R] question about the tisPlot function in package tis

2008-12-02 Thread Gabor Grothendieck
The development version of zoo has a tis to zoo conversion function that would allow you to plot your data using plot.zoo and xplot.zoo in which case you can use classic graphics axis function (for plot.zoo) or lattice facilities (with xyplot.zoo):

Re: [R] question about the tisPlot function in package tis

2008-12-02 Thread Gabor Grothendieck
Sorry that should have been: library(zoo) z - as.zoo(dat, class = yearmon) plot(z) xyplot(z) On Tue, Dec 2, 2008 at 2:35 PM, Gabor Grothendieck [EMAIL PROTECTED] wrote: The development version of zoo has a tis to zoo conversion function that would allow you to plot your data using plot.zoo and

Re: [R] Bug in transform?

2008-12-02 Thread Peter Dalgaard
Prof Brian Ripley wrote: As the help page says If some of the values are not vectors of the appropriate length, you deserve whatever you get! Yes (did I write that?). It is a bit annoying with things that almost work, though. [snip] I often need to use this for creating new

Re: [R] question on lmer function

2008-12-02 Thread Richard Palmer
Thanks Harold, I will review the lmer vignette again. On Tue, Dec 2, 2008 at 2:12 PM, Doran, Harold [EMAIL PROTECTED] wrote: Richard There is much more information that you need to provide before a thoughtful answer can be provided. Maybe you can describe the structure of your data, your

Re: [R] match

2008-12-02 Thread Peter Alspach
Tena koe Robert You could first apply paste apply(types, 1, paste, collapse=':') and then match. HTH Peter Alspach -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of threshold Sent: Wednesday, 3 December 2008 8:32 a.m. To:

Re: [R] Compiling R for 64-bit Windows (was Re: increasing memory limit in Windows Server 2008 64-bit)

2008-12-02 Thread David M Smith
Michael Styer wrote: So my next question is, does anyone have any thoughts about how significant a project it would be to compile R for 64-bit windows (using, e.g., the Portland Group compiler)? Based on our experiences at REvolution Computing, it's quite a significant project. We're the

[R] Suppressing tick labels?

2008-12-02 Thread Graham Smith
I am trying to suppress the tick labels on the x-axis of the following: barchart(richness[Wood==V]~Sample[Wood==V]) I have tried col.axis=white I have tried removing the axis all together with axes=FALSE I have tried xaxt=n I have also tried labels=c(label1, label2) to replace the default

Re: [R] linear functional relationships with heteroscedastic non-Gaussian errors - any packages around?

2008-12-02 Thread Spencer Graves
Isn't this a special case of structural equation modeling, handled by the 'sem' package? Spencer Jarle Brinchmann wrote: Thanks for the reply! On Tue, Dec 2, 2008 at 6:34 PM, Prof Brian Ripley [EMAIL PROTECTED] wrote: I wonder if you are using this term in its correct

Re: [R] linear functional relationships with heteroscedastic non-Gaussian errors - any packages around?

2008-12-02 Thread Jarle Brinchmann
Yes I think so if the errors were normally distributed. Unfortunately I'm far from that but the combination of sem its bootstrap is a good way to deal with it in the normal case. I must admit as a non-statistician I'm a not 100% sure what the difference (if there is one) between a linear

[R] (no subject)

2008-12-02 Thread jestinah mahachie
I have been using rpart to get trees. May someone help me with how I can do some logistic regression modelling to derive odds ratios on the branches. Also help on how to derive estimates on prediction error. -- Kind Regards, Jestinah Home Address: Brusselsepoortstraat 85 9000 Gent Belgium

[R] date operations

2008-12-02 Thread Christophe Dutang
Hi all, I'm dealing with dates in R (2.7.2), but some basic operations raise a warning. Incompatible methods (+.Date, Ops.difftime) for + I saw this topic in this mailing list, but I do not understand what to do... cf. https://stat.ethz.ch/pipermail/r-help/2008-June/165842.html Do I have

Re: [R] optimization problem (optim vs. nlminb)

2008-12-02 Thread Mike Prager
In case anyone is still reading this thread, I want to add this: In a current problem (a data-shy five-parameter nonlinear optimization), I found nlminb markedly more reliable than optim with method L-BFGS-B. In reviewing the fit I made, I found that optim only came close to its own minimum in

[R] QQ plots and boxcox

2008-12-02 Thread Erin Hodgess
Dear R People: In the DASL library, there is a story about hot dogs. Here are the data: Beef186 495 Beef181 477 Beef176 425 Beef149 322 Beef184 482 Beef190 587 Beef158 370 Beef139 322 Beef175 479 Beef148 375 Beef

[R] ggplot2 facet_wrap problem

2008-12-02 Thread stephen sefick
Hadley, I don't know if I am doing something wrong or if it is ggplot please see the two graphs at the bottom of the page (code). melt.nut - (structure(list(RiverMile = c(119L, 119L, 119L, 119L, 119L, 119L, 119L, 119L, 119L, 148L, 148L, 148L, 148L, 148L, 148L, 148L, 179L, 179L, 179L, 179L, 179L,

Re: [R] Suppressing tick labels?

2008-12-02 Thread William Dunlap
[R] Suppressing tick labels? Graham Smith myotisone at gmail.com Tue Dec 2 21:26:03 CET 2008 I am trying to suppress the tick labels on the x-axis of the following: barchart(richness[Wood==V]~Sample[Wood==V]) Use the scales= argument. E.g. d-data.frame(richness=log(1:12),

Re: [R] QQ plots and boxcox

2008-12-02 Thread John Fox
Dear Erin, What you characterize as snaky might also be called bimodal and short-tailed. Try, e.g., plot(density(residuals(dog1.aov))) to see the bimodality more clearly. A Box-Cox transformation can correct skewness, but won't help here. I hope this helps, John --

Re: [R] question about the tisPlot function in package tis

2008-12-02 Thread Jeffrey J. Hallman
[EMAIL PROTECTED] writes: List, I am using the 'tisPlot' function in Jeff Hallman's excellent tis package and was hoping that someone could spare me from having to dig into the code of his 'tisPlot' function. So far as I can tell, the preferred method of controlling the plotting of the

[R] simultaneous plots

2008-12-02 Thread David Epstein
Is there a good and concise way of making simultaneous plots that are identical, but directed to different devices? I'm writing an R-script that produces a pdf file. I would really like to check visually whether the pdf file shows what I expect. So I would like the same commands to produce a

Re: [R] fGarch and is.na()

2008-12-02 Thread Yohan Chalabi
Thanks for the report, the problem boils down to the call of methods:::bind_activation(TRUE) in one of the depended package. I can reproduce the problem with methods:::bind_activation(TRUE) dfr - data.frame(matrix(0, nrow = 1 , ncol = 1000)) dfr2 - is.na(dfr) I will forward you remark to

[R] sampling from data.frame

2008-12-02 Thread axionator
Hi all, I have a data frame with clustered rows as follows: Cu1 x1 y1 z1 ... Cu1 x2 y2 z2 ... Cu1 x3 y3 z3 ... # end of first cluster Cu1 Cu2 x4 y4 z4 ... Cu2 x5 y5 z5 Cu2 ... # end of second cluster Cu2 Cu3 ... ... cluster-size is 3 in the example above (rows making up a

[R] ggplot2 45deg axis labels

2008-12-02 Thread stephen sefick
I would like to rotate the axis labels 45 deg. -- Stephen Sefick Let's not spend our time and resources thinking about things that are so little or so large that all they really do for us is puff us up and make us feel like gods. We are mammals, and have not exhausted the annoying little

Re: [R] date operations

2008-12-02 Thread jim holtman
'str' is your friend. Look at the results: alpha2=as.Date(2008-12-21) alpha1=as.Date(2008-12-21)-as.Date(2008-10-26) alpha1 Time difference of 56 days alpha2 [1] 2008-12-21 alpha1+alpha2 Time difference of 14290 days Warning message: Incompatible methods (Ops.difftime, +.Date) for +

Re: [R] date operations

2008-12-02 Thread jim holtman
If you want to do the addition, 'unclass' the variable: alpha2+4 [1] 2008-12-25 alpha2 + unclass(alpha1) [1] 2009-02-15 On Tue, Dec 2, 2008 at 4:10 PM, Christophe Dutang [EMAIL PROTECTED] wrote: Hi all, I'm dealing with dates in R (2.7.2), but some basic operations raise a warning.

Re: [R] ggplot2 facet_wrap problem

2008-12-02 Thread ONKELINX, Thierry
Hi Stephen, I think you will need to clarify what your problem is with the second plot. HTH, Thierry -Oorspronkelijk bericht- Van: [EMAIL PROTECTED] namens stephen sefick Verzonden: di 2-12-2008 22:52 Aan: hadley wickham; R-help Onderwerp: [R] ggplot2 facet_wrap problem Hadley, I

Re: [R] ggplot2 45deg axis labels

2008-12-02 Thread ONKELINX, Thierry
Hi Stephen, Have a look at the ggplot2 book. You will find an example on p. 6 of http://had.co.nz/ggplot2/book/polishing.pdf HTH, Thierry -Oorspronkelijk bericht- Van: [EMAIL PROTECTED] namens stephen sefick Verzonden: wo 3-12-2008 0:31 Aan: R-help Onderwerp: [R] ggplot2 45deg axis

Re: [R] ggplot2 45deg axis labels

2008-12-02 Thread hadley wickham
Hi Stephen, Have a look at the polishing your plot for publication chapter of the ggplot2 book. Regards, Hadley On Tue, Dec 2, 2008 at 5:31 PM, stephen sefick [EMAIL PROTECTED] wrote: I would like to rotate the axis labels 45 deg. -- Stephen Sefick Let's not spend our time and resources

Re: [R] sampling from data.frame

2008-12-02 Thread jim holtman
Not sure exactly what you mean by 'sample' since you did not provide an example of the expected output, or input data that could be used. Here is an example of taking one sample from each cluster: df - data.frame(id=paste(C, rep(1:5, each=3), sep=''), data=1:15) # sample 1 from each cluster

Re: [R] ggplot2 facet_wrap problem

2008-12-02 Thread stephen sefick
If you look at the TSS graph in the faceted example and then look at the plot of just the GPP vs. TSS. They are different graphs all together. The one that is not faceted is correct. On Tue, Dec 2, 2008 at 6:36 PM, ONKELINX, Thierry [EMAIL PROTECTED] wrote: Hi Stephen, I think you will need

[R] help on tapply using sample with differing sample-sizes

2008-12-02 Thread Jourdan Gold
Hello, My question likely got buried so I am reposting it in the hopes that someone has an answer. I have thought more about the question and modified my question. I hope tha my specific question is: I am attempting to create a bootstrap procedure for a finite sample using the theory of Rao

[R] Help with maps

2008-12-02 Thread Avram Aelony
A few questions about maps... (1) How can I find a listing of the internal data sets that map() from the maps library contains? For example, usa, county, state, nz all work. Are there any others? (2) Is there an easier, more generalized way to produce this

Re: [R] Help with maps

2008-12-02 Thread hadley wickham
On Tue, Dec 2, 2008 at 6:21 PM, Avram Aelony [EMAIL PROTECTED] wrote: A few questions about maps... (1) How can I find a listing of the internal data sets that map() from the maps library contains? For example, usa, county, state, nz all work. Are there any others? help(package = maps)

[R] how to handle irregularly spaced data as timeseries

2008-12-02 Thread Kirk Wythers
I have a set of modeled climate data recorded at irregular intervals. The format of the data is such that there are monthly measurements for the years 2000, 2020, 2050, 2080, etc. Therefore I have 12 regular records, a skip of some number of years, then 12 more monthly records, another

Re: [R] how to handle irregularly spaced data as timeseries

2008-12-02 Thread Gabor Grothendieck
Have a look at the zoo package. There are three vignettes (pdf documents) included with the package that give many examples of its use. Also see ?read.zoo, ?plot.zoo and ?xyplot.zoo You will need something like: library(zoo) z - read.zoo(myfile, ...whatever...) plot(z) On Tue, Dec 2, 2008 at

[R] for loop looking for file names

2008-12-02 Thread Steven Kennedy
I have a series of csv files in several folders. All begin with a 7 digit number and end with the letter E (eg. 0726016E.csv). I want to be able to read a file in to R, take some of the data out of it and store it in a matrix, then move on to the next file and do the same thing. I was planning

Re: [R] for loop looking for file names

2008-12-02 Thread jim holtman
To get the file names in the current directory try: list.files(pattern=[[:digit:]]{7}E) On Tue, Dec 2, 2008 at 4:11 PM, Steven Kennedy [EMAIL PROTECTED] wrote: I have a series of csv files in several folders. All begin with a 7 digit number and end with the letter E (eg. 0726016E.csv). I

Re: [R] Help with maps

2008-12-02 Thread Avram Aelony
On Tuesday, December 02, 2008, at 04:40PM, hadley wickham [EMAIL PROTECTED] wrote: On Tue, Dec 2, 2008 at 6:21 PM, Avram Aelony [EMAIL PROTECTED] wrote: A few questions about maps... (1) How can I find a listing of the internal data sets that map() from the maps library contains? For

  1   2   >