Re: [R] Sweaving quotes

2010-07-28 Thread Prof Brian Ripley
The alternative is to tell LaTeX what encoding the file is in. For those using UTF-8 locales this means adding the line \usepackage[utf8]{inputenc} Now Murray mentions 'Vista', and so is presumably using cp1252 (the Western-European-language Windows default). That is spported by inputenc, s

Re: [R] Match() on raw objects ?

2010-07-28 Thread bruno Piguet
On Tue, 27 Jul 2010, Prof Brian Ripley wrote : > > Note that raw values in R are really intended to be passed around and not > manipulated: if you do much of the latter, coercing to integer, say, is > likely to be much more efficient. > > Indeed : I tried a minimal benchmark : a <- charToRaw(pas

[R] Displaying Counts of Unused Factors in Contingency Tables with table()

2010-07-28 Thread Na'im R. Tyson
R-philes, I have a question about displaying counts of unused factors using the table() function. I have two vectors with character data in them: local.labels("ah", "ah", "ah~") local.preds("ah", "ah", "ah") If I use the table function as shown below, I get an error because the number of

Re: [R] how to code it??

2010-07-28 Thread Jonathan Christensen
On Wed, Jul 28, 2010 at 2:18 PM, Henrique Dallazuanna wrote: > You've tried: > > diff(c(0, x)) ? > This is clever, but not quite what he's asking for--it converts a sequence of 1's into a 1 followed by zeroes. Jonathan > > On Wed, Jul 28, 2010 at 3:10 PM, Raghu wrote: > > > Hi > > > > I have

[R] package xpose4 in Vista - Update

2010-07-28 Thread Santosh
Dear R experts, There seems to be a problem (please see the error messages below) with the installation of the latest version of xpose4 (version 4.2.1) with the latest R version.(2.11.1)... I didn't face such installation problems when using version R.2.10.1 Would really appreciate your assistanc

[R] Function to return variable name

2010-07-28 Thread Jeremy Miles
I'd like a function that returns the variable name. As in: MyData$Var1 Would return: Var1 There should be a straightforward way to do this, but I can't see it. Thanks, Jeremy -- Jeremy Miles Psychology Research Methods Wiki: www.researchmethodsinpsychology.com __

Re: [R] Pattern recognition

2010-07-28 Thread Pablo Cerdeira
This is a good example of what I'm looking for: [image: dendrogram.jpg] Best On Thu, Jul 29, 2010 at 12:01 AM, Pablo Cerdeira wrote: > > Dear all, > > I'm trying to use some technic to do a pattern recognition over a large > dataset. I really don't have any idea on how to do that using R. > >

[R] Pattern recognition

2010-07-28 Thread Pablo Cerdeira
Dear all, I'm trying to use some technic to do a pattern recognition over a large dataset. I really don't have any idea on how to do that using R. Here is a sample of the data: id,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15 1480010,208,69,180,465,465,241,241,69,584,26,75,578,507,75,284 1480183,208,69,35

Re: [R] Sweaving quotes

2010-07-28 Thread Murray Jorgensen
Thanks, Marc, and also Jay Kearns. After experimentation I think useFancyQuotes = FALSE my be best as TeX style can look a bit funny for R output. Cheers, Murray Marc Schwartz wrote: On Jul 28, 2010, at 7:43 PM, Murray Jorgensen wrote: The significance code line to summary() applied to an

[R] Statistical mailing list

2010-07-28 Thread Ralf B
I am looking for a mailing list for general statistical questions that are not R related. Do you have any suggestions for lists that are busy and helpful and/or lists that you use and recommend? Thanks in advance, Ralf __ R-help@r-project.org mailing li

Re: [R] Sweaving quotes

2010-07-28 Thread Marc Schwartz
On Jul 28, 2010, at 7:43 PM, Murray Jorgensen wrote: > The significance code line to summary() applied to an lm() fitted model > object is > > Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 > > The corresponding line in the LaTeX source produced by Sweave is > > Signif. codes:

Re: [R] Axes=F and plotting dual y axes

2010-07-28 Thread P Ehlers
Johnson, Cedrick W. wrote: That worked. Stupid me forgot that I had the stock ticker 'F' assigned in my workspace. Well.. guess I'll hit myself with a 2x4 now.. Thanks for your help guys.. No, don't do that. Instead, calculate how much time you saved by typing 'F' instead of 'FALSE' and how

[R] Sweaving quotes

2010-07-28 Thread Murray Jorgensen
The significance code line to summary() applied to an lm() fitted model object is Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 The corresponding line in the LaTeX source produced by Sweave is Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 which looks the same

Re: [R] Variance-covariance matrix from GLM

2010-07-28 Thread Bojuan Zhao
Fantastic! it's solved! Thank you very much Bill! Barbara --- On Wed, 7/28/10, bill.venab...@csiro.au wrote: > From: bill.venab...@csiro.au > Subject: RE: [R] Variance-covariance matrix from GLM > To: bojuanz...@yahoo.com, r-help@r-project.org > Date: Wednesday, July 28, 2010, 8:01 PM > ?vco

[R] ggplot2 histograms... a subtle error found

2010-07-28 Thread Mike Williamson
Hello all, I have a peculiar and particular bug that I stumbled across with ggplot2. I cannot seem to replicate it with anything other than my specific data set. Here is the problem: - when I try to plot a histogram, allowing for ggplot2 to decide the binwidths itself, I get the f

Re: [R] Variance-covariance matrix from GLM

2010-07-28 Thread Bill.Venables
?vcov ### now in the stats package You would use V <- vcov(my.glm) -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of Bojuan Zhao Sent: Thursday, 29 July 2010 9:52 AM To: r-help@r-project.org Subject: [R] Variance-covariance matr

[R] Variance-covariance matrix from GLM

2010-07-28 Thread Bojuan Zhao
Hello, Is there a way to obtain the variance-covariance matrix of the estimated parameters from GLM? my.glm<-glm(mat ~X,family = binomial, data =myDATA) out1<-predict(my.glm,se.fit = TRUE) std<-out1$se.fit se.fit is for getting the standard errors of the estimated parameters (\betas). Is the

Re: [R] pass list of args to function call

2010-07-28 Thread Henrique Dallazuanna
Call the function with quotes, so: model.fun <- "glm" rslt <- do.call(model.fun,model.opt) On Wed, Jul 28, 2010 at 7:30 PM, darckeen wrote: > > Thanks for your helpful response. > > I have a related question, i would like to further generalize the function > by having the call arg of do.call f

Re: [R] read.delim()

2010-07-28 Thread Doran, Harold
Thank you, Phil. Unfortunately, there are quotes used properly elsewhere. - Original Message - From: Phil Spector To: Doran, Harold Cc: r-help@r-project.org Sent: Wed Jul 28 18:29:32 2010 Subject: Re: [R] read.delim() Harold - If there aren't any true quoted fields in the file, you

Re: [R] Sweave and scan()

2010-07-28 Thread Murray Jorgensen
Hi Duncan, I at first thought that your suggestion was how the .tex file should be and that this would involve tediously editing the latex source each time I updated the .Rnw. But as you probably intended I find that I can slip the Schunks and Sinputs into the .Rnw and all still works. Thank

[R] Reading timestamp column from MySQL

2010-07-28 Thread harsh yadav
Hi, I am reading a SQL (MySQL) table in R data frame. When I read in the table that has a timestamp data-type field, R gives it the following format:- 1.236887e+12 So when I want to manipulate a column with timestamp = 1236887146615 It returns me multiple rows, as many timestamps gets converte

Re: [R] Sweave and scan()

2010-07-28 Thread Duncan Murdoch
On 28/07/2010 6:33 PM, Murray Jorgensen wrote: Omigod! The archival links shows that this was the same problem that caused me to give up on Sweave about 6 years ago. I guess I never properly assimilated Brian Ripley's comments at the time. I finished up doing this: \begin{verbatim} > height

Re: [R] Beginner stucked with raster + geoR package.

2010-07-28 Thread Alaios
Can someone help me please with step 5? r <- setValues(r,temp) # Unfortunately this ends with the message: Error in setValues(r, temp) : values must be a vector How to make temp a vector? To: r-help@r-project.org Sent: Wed, July 28, 2010 3:55:54 PM Subject: [

Re: [R] Sweave and scan()

2010-07-28 Thread Murray Jorgensen
Omigod! The archival links shows that this was the same problem that caused me to give up on Sweave about 6 years ago. I guess I never properly assimilated Brian Ripley's comments at the time. I finished up doing this: \begin{verbatim} > height = scan() 1: 64 62 66 65 5: 62 69 72 72 70 10: R

Re: [R] pass list of args to function call

2010-07-28 Thread darckeen
Thanks for your helpful response. I have a related question, i would like to further generalize the function by having the call arg of do.call function being referenced. Like this: model.fun <- glm model.opt <- list(y1~1,data=df,family="binomial") rslt <- do.call(model.fun,model.opt) This wo

Re: [R] read.delim()

2010-07-28 Thread Phil Spector
Harold - If there aren't any true quoted fields in the file, you could pass the quote="" option to read.delim(). - Phil Spector Statistical Computing Facility Department

Re: [R] [Rd] R CMD build wiped my computer

2010-07-28 Thread Ray Brownrigg
On Thu, 29 Jul 2010, Duncan Murdoch wrote: > On 28/07/2010 10:01 AM, Jarrod Hadfield wrote: > > Hi Marc, > > > > Thanks for the info on recovery - most of it can pieced together from > > backups but a quick, cheap and easy method of recovery would have been > > nicer. > > > > My main concern is tha

[R] read.delim()

2010-07-28 Thread Doran, Harold
I am reading in a very large file with names in it and R is truncating the number of rows it reads in. The separator in this file is a pipe '|' and so I use dat <- read.delim('pathToMyFile', header= TRUE, sep='|') It turns out that it is reading up to row 61145 and stopping and I think I see w

[R] Model fit

2010-07-28 Thread rwanuza
Does anyone know how to calculated a AIC, BIC, score when using svyglm for logistic regressions? Thanks Rwanuza -- View this message in context: http://r.789695.n4.nabble.com/Model-fit-tp2305668p2305668.html Sent from the R help mailing list archive at Nabble.com. [[alternative HTML v

[R] Adding lines to two graphs alternately

2010-07-28 Thread Dennis Fisher
Colleagues I do a series of calculations, then add some output as a line in a graphic. This is repeated a large number of times. Then, I move to the next panel and do a similar set of calculations, add lines, Many of the calculations are identical for the two panels, so the calculations

[R] RPostgreSQL 0.1-6 installation trouble

2010-07-28 Thread ericksonsm
I am running Mac OS X version 10.6.3. I am running R 2.11.1 GUI 1.34 64 bit. I have RPostgreSQL 0.1-6 (as RPostgreSQL_0.1-6.tar) from this site: http://mac.softpedia.com/get/Developer-Tools/RPostgresql.shtml I have PostgreSQL version 8.4.4-1 for Mac OS X from this site: http://www.enterprise

Re: [R] how to code it??

2010-07-28 Thread raghu
Thanks Gabor. On Wed, Jul 28, 2010 at 9:24 PM, Gabor Grothendieck [via R] < ml-node+2305561-610338075-309...@n4.nabble.com > wrote: > On Wed, Jul 28, 2010 at 2:10 PM, Raghu <[hidden > email]> > wrote: > > > Hi > > > > I have say a large vect

Re: [R] how to code it??

2010-07-28 Thread Raghu
Thanks Matt, I will try putting the other rule into this. On Wed, Jul 28, 2010 at 9:06 PM, Matt Shotwell wrote: > If I take your meaning correctly, you want something like this. > > x <- c(0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 1, 1, 0, 0, 1, 0, 0, 0, > + 1) > > easy <- function(x) { > + state <-

Re: [R] error in f(x,...)

2010-07-28 Thread Gray Calhoun
Hi Nathalie, It should be "subdivisions" in the line: > for (i in 1:length(x)){res[i] <- integrate(f,x[i],2.967, > subdivision=2000)$value} R is trying to pass the argument "subdivision" to your function f. --Gray On Wed, Jul 28, 2010 at 2:27 PM, Nathalie Gimenes wrote: > Dear all, > > I tr

Re: [R] how to code it??

2010-07-28 Thread raghu
Thanks Henrique. On Wed, Jul 28, 2010 at 9:20 PM, Henrique Dallazuanna [via R] < ml-node+2305552-1908682012-309...@n4.nabble.com > wrote: > You've tried: > > diff(c(0, x)) ? > > On Wed, Jul 28, 2010 at 3:10 PM, Raghu <[hidden > email]> > wro

Re: [R] Add an arrow to a plot

2010-07-28 Thread Trying To learn again
Hi many thank¡¡¡ I promise work hard¡¡¡ Today I ´ve been very busy but I will wake up tomorrow on 5:00 to see what you have suggested I need work harder but working and having a children (1,5 years is very exhaust 2010/7/28 Greg Snow > In addition to the arrows function, look at the my.symbols

Re: [R] [Rd] R CMD build wiped my computer

2010-07-28 Thread Kevin Wright
Isn't the trash a feature of the Gnome GUI, not the shell? Using "rm" has always scared me, so I've recently configured my system thus: Install trash-cli. See http://code.google.com/p/trash-cli/ Alias "del" to trash-put. Set a cron job to clean trash after xx

Re: [R] [Rd] R CMD build wiped my computer

2010-07-28 Thread Marc Schwartz
On Jul 28, 2010, at 3:18 PM, Bos, Roger wrote: > Ubuntu also uses ~ as a backup file syntax, but Ubuntu has a trash can > where deleted files are located, so it would be easy to restore them. I > would be surprised if Fedora didn't also have a trash can. Roger, It does, but that is only if y

Re: [R] pass list of args to function call

2010-07-28 Thread Henrique Dallazuanna
Try this: do.call(glm, list(y ~ 1, data = df, family = 'binomial')) On Wed, Jul 28, 2010 at 5:23 PM, darckeen wrote: > > I was wondering how i can get a list of arguments to evaluate within > function > call. I'd like to be able do something like this: > > > #mod <- glm(y~1,data=df,family="b

[R] pass list of args to function call

2010-07-28 Thread darckeen
I was wondering how i can get a list of arguments to evaluate within function call. I'd like to be able do something like this: #mod <- glm(y~1,data=df,family="binomial") opt <- list(family="binomial") mod <- glm(y~1,data=df,opt) Any idea how this can be done properly? TIA -- View this mes

Re: [R] how to code it??

2010-07-28 Thread Gabor Grothendieck
On Wed, Jul 28, 2010 at 2:10 PM, Raghu wrote: > Hi > > I have say a large vector of 3500 digits. Initially the digits are 0s and > 1s. I need to check for a rule to change some of the 0s to -1s in this > vector. But once I change a 0 to -1 then I need to start applying the rule > to change the nex

Re: [R] [Rd] R CMD build wiped my computer

2010-07-28 Thread Bos, Roger
Ubuntu also uses ~ as a backup file syntax, but Ubuntu has a trash can where deleted files are located, so it would be easy to restore them. I would be surprised if Fedora didn't also have a trash can. -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project

Re: [R] how to code it??

2010-07-28 Thread Henrique Dallazuanna
You've tried: diff(c(0, x)) ? On Wed, Jul 28, 2010 at 3:10 PM, Raghu wrote: > Hi > > I have say a large vector of 3500 digits. Initially the digits are 0s and > 1s. I need to check for a rule to change some of the 0s to -1s in this > vector. But once I change a 0 to -1 then I need to start appl

Re: [R] Documenting different OO-aproaches in R as a package?

2010-07-28 Thread Martin Maechler
> "s" == schuster > on Tue, 27 Jul 2010 22:17:09 +0200 writes: s> Hello, s> I see some people including myself confused by the s> different object-oriented approaches in R (S3, S4, OOP, s> R.oo etc.). s> Would it be ok to collect examples and solutions for the

Re: [R] how to code it??

2010-07-28 Thread Matt Shotwell
If I take your meaning correctly, you want something like this. > x <- c(0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 1, 1, 0, 0, 1, 0, 0, 0, + 1) > easy <- function(x) { + state <- 0 + for (i in 1:length(x)) { + if (x[i] == 0) + x[i] <- state + state <- 0 + if (x

Re: [R] question about SVM in e1071

2010-07-28 Thread Jack Luo
Pau, Sorry for getting back to you for this again. I am getting confused about your interpretation of 3). It is obvious from your code that increasing C results in* smaller *number of SVs, this seems to contradict with your interpretation " * Increasing the value of C (...) forces the creation of

[R] install packages in Windows Vista

2010-07-28 Thread Santosh
Dear R experts... I would really appreciate your suggestions in installing a package in Windows Vista... I am unable to install a package on a windows vista based computer, in spite of running it as an administrator. The package "xpose4.2.1" is available from " https://sourceforge.net/projects/xpo

Re: [R] how to generate a random data from a empirical distribition

2010-07-28 Thread Frank Harrell
This is true by definition. Read about the bootstrap which may give you some good background information. Frank E Harrell Jr Professor and ChairmanSchool of Medicine Department of Biostatistics Vanderbilt University On Wed, 28 Jul 2010, xin wei wrote: hi, F

Re: [R] Odd crash with tcl/tk

2010-07-28 Thread Gabor Grothendieck
On Wed, Jul 28, 2010 at 10:06 AM, AndrewPage wrote: > > Hi, > > Recently, I've been trying to use packages in R that require loading the > Tcl/Tk interface.  However, I get a strange result and a crash that I > haven't been able to find discussion about on these boards (or any others). > > When I

[R] error in f(x,...)

2010-07-28 Thread Nathalie Gimenes
Dear all, I tried once to create one variable called bip such that: bip <- cip + (1/f(cip))*fi(f,cip) And this was working. But now, doing the same thing I did before, the software shows me the following message: Error in f(x, ...) : unused argument(s) (subdivision = 2000) I have the variable c

[R] Help Creating a Stacked Bar Chart with Color Coding

2010-07-28 Thread Majonu
I have a data set of the following form: Johnson 4 Smith4 Smith2 Smith3 Garcia 1 Garcia 4 Rodriguez 2 Adams 2 Adams 3 Adams 4 Turner 4 Turner 3 And I'd like to create a stacked bar chart that has scores on the x-axis an

Re: [R] randomisation for matrix

2010-07-28 Thread Knut Krueger
Gray Calhoun schrieb: Hi Knut, I think you're going to have to be more specific. The code matrix(rnorm(25), 5, 5) I found the answer there is a "generator seed" field in Ucinet, I do not know why its possible to set the random generator starting point, by hand and I assume the generator s

Re: [R] how to generate a random data from a empirical distribition

2010-07-28 Thread xin wei
hi, Frank: how can we make sure the randomly sampled data follow the same distribution as the original dataset? i assume each data point has the same prabability to be selected in a simple random sampling scheme. thanks -- View this message in context: http://r.789695.n4.nabble.com/how-to-gene

[R] Odd crash with tcl/tk

2010-07-28 Thread AndrewPage
Hi, Recently, I've been trying to use packages in R that require loading the Tcl/Tk interface. However, I get a strange result and a crash that I haven't been able to find discussion about on these boards (or any others). When I enter library(tcltk), it reads "Loading Tcl/Tk interface ... ", bu

Re: [R] Random Forest - Strata

2010-07-28 Thread Coll
Max, Thanks. Yes what you said is exactly I am looking for, i.e. the first tree fits using data from sites A&B, then predicts on C (and so on). Does that means if I : 1. pass this list as index into trainControl > tmpSiteList [[1]] [1] 1 2 3 4 5 6 7 [[2]] [1] 1 2 3 8 9 10 [[3]] [1] 4 5

[R] how to code it??

2010-07-28 Thread Raghu
Hi I have say a large vector of 3500 digits. Initially the digits are 0s and 1s. I need to check for a rule to change some of the 0s to -1s in this vector. But once I change a 0 to -1 then I need to start applying the rule to change the next 0 only after I see the next 1 in the vector. Say for ex

[R] Out-of-sample predictions with boosting model

2010-07-28 Thread Travis Berge
Hi UseRs - I am new to R, and could use some help making out-of-sample predictions using a boosting model (the mboost command). The issue is complicated by the fact that I have panel data (time by country), and am estimating the model separately for each country. FYI, this is monthly data and I ha

Re: [R] Introductory statistics and introduction to R

2010-07-28 Thread Gabor Grothendieck
On Tue, Jul 27, 2010 at 10:46 AM, Marsh Feldman wrote: > Hi, > > I have a bright, diligent second-year graduate student who wants to learn > statistics and R and will, in effect, be taking a tutorial from me on these > subjects. (If you've seen some of my questions on this list, please don't > lau

Re: [R] columns mapping

2010-07-28 Thread jd6688
Thank you so much. it worked as expected. you have been great help -- View this message in context: http://r.789695.n4.nabble.com/columns-mapping-tp2305213p2305401.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org m

Re: [R] hatching posibility in Panel.Polygon

2010-07-28 Thread Greg Snow
Tufte discusses hash lines in his book "The Visual Display of Quantitative Information" and does a better job of it than I can. The short version is that the hashing can actually produce optical illusion effects that distort the information. (and often don't copy well either). Printing and co

Re: [R] Introductory statistics and introduction to R

2010-07-28 Thread Matthew Keller
Hi Marsh, I taught an intro to R course and have posted all the materials up on the web: http://psych-swiki.colorado.edu:8080/LearnR. Most learning in R comes from doing, not reading, and that's how I structured my course. All the lectures/HWs can be done individually, and the keys are there to c

[R] R CMD build wiped my computer

2010-07-28 Thread R P Herrold
On Wed, 28 Jul 2010, Jarrod Hadfield wrote: My main concern is that this could happen again and that the "bug" is not limited to R 2.9. I would think that an accidental carriage return at the end of a file name (even a temporary one) would be a reasonably common phenomenon (I'm surprised I ha

Re: [R] Axes=F and plotting dual y axes

2010-07-28 Thread Johnson, Cedrick W.
That worked. Stupid me forgot that I had the stock ticker 'F' assigned in my workspace. Well.. guess I'll hit myself with a 2x4 now.. Thanks for your help guys.. -c On 07/28/2010 12:37 PM, Ben Bolker wrote: Johnson, Cedrick W. cedrickjohnson.com> writes: Howdy. Been running into a bit o

Re: [R] finding the next highest number in an array

2010-07-28 Thread Bert Gunter
... just a note: you don't have to first sort the vector to do this: > x <- sample(1:7) > x [1] 3 5 7 6 2 4 1 > which(x==min(x[x>4])) [1] 2 Bert Gunter Genentech Nonclinical Biostatistics On Wed, Jul 28, 2010 at 3:12 AM, Raghu wrote: > Hi > > I have a sorted array ( in ascending order) and I

Re: [R] Introductory statistics and introduction to R

2010-07-28 Thread G. Jay Kerns
Dear Marsh, On Tue, Jul 27, 2010 at 10:46 AM, Marsh Feldman wrote: > Hi, > > I have a bright, diligent second-year graduate student who wants to learn > statistics and R and will, in effect, be taking a tutorial from me on these > subjects. (If you've seen some of my questions on this list, pleas

Re: [R] Axes=F and plotting dual y axes

2010-07-28 Thread D Kelly O'Day
Cedrick I used this script to produce a simple chart with no axes: series2 = c(1:50) series1 = rep(25:74) plot(series1, series2, main="Woo", col="red", xlab="", ylab="", axes=F, type="l") Works fine on Windows XP using R 2.11.1. -- View this message in context: http://r.789695.n4.n

Re: [R] memory problem for scatterplot using ggplot

2010-07-28 Thread David Winsemius
On Jul 28, 2010, at 9:53 AM, Brandon Hurr wrote: It was my understanding that R wasn't really the best thing for absolutely huge datasets. 17.5 million points would probably fall under the category of "absolutely huge." I'm on a little netbook right now (atom/R32) and it failed, but I'll

Re: [R] problem with building package on CRAN

2010-07-28 Thread Uwe Ligges
Package should be online now (or within 24 hours at least). Apologies for the noise, and please just ask me directly in case you get my messages - no need to flood R-help even more Best wishes, Uwe On 26.07.2010 23:25, Romain Francois wrote: Hello, (ccing Rcpp-devel too because this

Re: [R] Axes=F and plotting dual y axes

2010-07-28 Thread Ben Bolker
Johnson, Cedrick W. cedrickjohnson.com> writes: > > Howdy. Been running into a bit of trouble with plotting. Seems that > axes=F is not "working". Whenever I plot (either a dataframe or xts/zoo > series) and I set axes=F along with xlab/ylab="" I still get the default > axes printed in my cha

Re: [R] columns mapping

2010-07-28 Thread Wu Gong
Hi, please try ?merge DF2$mappedColumn <- DF2$name merge(DF1,DF2,all.x=T,sort=F) - A R learner. -- View this message in context: http://r.789695.n4.nabble.com/columns-mapping-tp2305213p2305250.html Sent from the R help mailing list archive at Nabble.com. __

Re: [R] UseR! 2010 - my impressions

2010-07-28 Thread Uwe Ligges
Ravi, you are so right, we had the opportunity to attend an excellent and extremely well organized conference. Congratulation to Kate and the whole team at the NIST!!! Uwe ps. Thanks for starting the thread - it is too easy to forget about these well deserved thanks when we are back in nor

Re: [R] strange error : isS4(x) in gamm function (mgcv package). Variable in data-frame not recognized???

2010-07-28 Thread Joris Meys
Follow up: I finally succeeded to more or less reproduce the error. The origin lied in the fact that I accidently loaded a function while being in browser mode for debugging that function. So something went very much wrong with the namespaces. Teaches me right... Cheers Joris On Wed, Jul 28, 201

[R] Axes=F and plotting dual y axes

2010-07-28 Thread Johnson, Cedrick W.
Howdy. Been running into a bit of trouble with plotting. Seems that axes=F is not "working". Whenever I plot (either a dataframe or xts/zoo series) and I set axes=F along with xlab/ylab="" I still get the default axes printed in my chart. Consider this: #Create some sample data, both 50 units

[R] columns mapping

2010-07-28 Thread jd6688
DF1 name OTHER ABC O KKK O QQQ O DDD O PPP O DF2 name ABC KKK DDD If the names in df1 resides in df2, then add the mapped name to df1 as a separate column, for instance "mappedColumn" the output should be: DF1 name OTHER mappedColumn ABCOABC KKKO KKK QQQ

Re: [R] error: arguments imply differing number

2010-07-28 Thread Jonathan Christensen
Hi, Thanks for including code and data so that we could reproduce what you're doing. Your problem is that you tell ddply to split the dataset by runNumber and cat1, which results in 4 groups. ddply then applies my.summary() to these four groups. One of these groups (cat1 = 1 and runNumber=1) has

Re: [R] Optimization problem with nonlinear constraint

2010-07-28 Thread Ravi Varadhan
For those interested in esoteric of special functions, here is a nice reference on the Lambert W function: http://www.cs.uwaterloo.ca/research/tr/1993/03/W.pdf Ravi. -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of Hans W Borchers S

Re: [R] Optimization problem with nonlinear constraint

2010-07-28 Thread Ravi Varadhan
Very nice, Hans! I didn't know of the existence of Lambert W function (a.k.a Omega function) before. I didn't know that it occurs in the solution of exponential decay with delay: dy/dy = a * y(t - 1). Apparently it is more than 250 years old! Thanks, Ravi. -Original Message- From: r-he

Re: [R] [Rd] R CMD build wiped my computer

2010-07-28 Thread Duncan Murdoch
On 28/07/2010 10:01 AM, Jarrod Hadfield wrote: Hi Marc, Thanks for the info on recovery - most of it can pieced together from backups but a quick, cheap and easy method of recovery would have been nicer. My main concern is that this could happen again and that the "bug" is not limited to

Re: [R] Optimization problem with nonlinear constraint

2010-07-28 Thread Hans W Borchers
Uli Kleinwechter uni-hohenheim.de> writes: > > Dear Ravi, > > As I've already written to you, the problem indeed is to find a solution > to the transcendental equation y = x * T^(x-1), given y and T and the > optimization problem below only a workaround. I don't think optimization is the ri

Re: [R] How get colnames and rownames in Rcpp method?

2010-07-28 Thread Romain Francois
Hi, This is not the appropriate list for questions about Rcpp. See the Rcpp-devel mailing list, but first please think about a reproducible example. Romain Le 28/07/10 14:11, 나여나 a écrit : Hi all, How get colnames and rownames in Rcpp method? attecthed file : RGui.exe capture

Re: [R] memory problem for scatterplot using ggplot

2010-07-28 Thread Mark Connolly
On 07/28/2010 06:13 AM, Edwin Husni Sutanudjaja wrote: Dear all, I have a memory problem in making a scatter plot of my 17.5 million-pair datasets. My intention to use the "ggplot" package and use the "bin2d". Please find the attached script for more details. Could somebody please give me any c

Re: [R] randomisation for matrix

2010-07-28 Thread Gray Calhoun
Hi Knut, I think you're going to have to be more specific. The code matrix(rnorm(25), 5, 5) generates a random 5 by 5 matrix. If you need specific distributions, search through the helpfiles using help.search or RSiteSearch. --Gray On Wed, Jul 28, 2010 at 7:51 AM, Knut Krueger wrote: > Hi

Re: [R] [Rd] R CMD build wiped my computer

2010-07-28 Thread Jarrod Hadfield
Hi Marc, Thanks for the info on recovery - most of it can pieced together from backups but a quick, cheap and easy method of recovery would have been nicer. My main concern is that this could happen again and that the "bug" is not limited to R 2.9. I would think that an accidental carriag

[R] Beginner stucked with raster + geoR package.

2010-07-28 Thread Alaios
Hello everyone. I am trying to build up understanding in R by trying to develop just some simple scenarios. I would like to explain you what I am trying to do and what I did so far. I would like to put inside a RasterLayer (raster package) a Gaussian field (for given covariance) using grf f

Re: [R] How to point a column of dataframe by a "character"

2010-07-28 Thread Gabor Grothendieck
On Wed, Jul 28, 2010 at 8:59 AM, Joshua Wiley wrote: > Hi Tony, > > I am sure there are other ways, but I would create formula objects and > then pass them to lm().  Here's an example: > > mydata <- data.frame(Y = 1:10, X1 = 11:20, X2 = 21:30) > > my.names <- names(mydata)[-1] > > for(i in my.name

Re: [R] memory problem for scatterplot using ggplot

2010-07-28 Thread Brandon Hurr
It was my understanding that R wasn't really the best thing for absolutely huge datasets. 17.5 million points would probably fall under the category of "absolutely huge." I'm on a little netbook right now (atom/R32) and it failed, but I'll try it on my macbookPro/R64 later and see if it's able to

Re: [R] T and F (was: Optimization problem with nonlinear constraint)

2010-07-28 Thread Ravi Varadhan
Hi Patrick, Actually, I feel just the opposite, i.e. it is not a good idea to use `T' for `TRUE'. I have been snared by this trap many a times in my early days with S-Plus and R. It is a good practice to use unabbreviated values. Best, Ravi. -Original Message- From: r-help-boun...@r-

[R] [R-pkgs] IPSUR-1.0 is on CRAN, plus update to RcmdrPlugin.IPSUR

2010-07-28 Thread G. Jay Kerns
IPSUR-1.0 is making its way through CRAN. It is a snapshot of the development version of the following textbook: Title: Introduction to Probability and Statistics using R, First Edition ISBN: 978-0-557-24979-4 Publisher: me The book is targeted for an undergraduate course in probability and stati

[R] error: arguments imply differing number

2010-07-28 Thread jd6688
mydata <- read.table(textConnection(" Id cat1 location item_values p-values sequence a111 1 3002737 100 0.01 1 a112 1 3017821 102 0.05 2 a113 2 3027730 103 0.02 3 a114 2 3036220 104 0.04 4 a115 1 3053984

Re: [R] fisher's posthock test or fisher's combination test

2010-07-28 Thread Knut Krueger
I was looking for: fisher.comb <- function (pvalues) { df=length(pvalues) ch2=(-2*sum(log(pvalues))) return pchisq(ch2, df=df, lower.tail=FALSE) } http://en.wikipedia.org/wiki/Fisher%27s_method for the second part: the calculation of the p-value from the result of fischer combination test . Ki

Re: [R] Installing the Zoo Library.

2010-07-28 Thread Rnoobie
There is no problems getting zoo to install correctly on my windows 7 box, but I can't get it happening on linux (ubuntu) here is the latest issue again... > install.packages("zoo_1.6-4.tar.gz", repos = NULL ); Warning in install.packages("zoo_1.6-4.tar.gz", repos = NULL) : argument 'li

Re: [R] Tinn-R related problem

2010-07-28 Thread Knut Krueger
ForestStats schrieb: Hello, I too am having this problem. Some two minutes ago all was well then all of a sudden I cannot backspace or delete or use arrows etc.. there is a special Tinn-R forum: http://sourceforge.net/projects/tinn-r/forums Kind Regards Knut __

Re: [R] How get colnames and rownames in Rcpp method?

2010-07-28 Thread 나여나
Hi all, How get colnames and rownames in Rcpp method? attecthed file : RGui.exe capture my work environment : R version : 2.11.1 OS : WinXP Pro sp3 Thanks and best regards. Young-Ju, Park from Korea [1][rKWLzcpt.zNp8gmPEwGJCA00] [...@from=dllmain&rcpt

[R] randomisation for matrix

2010-07-28 Thread Knut Krueger
Hi to all, I am looking for a randomisation procedure for a single matrix, including a possibility to set the number of randomisations and the to set the number of row and columns . Knut __ R-help@r-project.org mailing list https://stat.ethz.ch/mai

Re: [R] Glm

2010-07-28 Thread Gray Calhoun
Ana, It's not really clear what you are trying to do. From your description, you aren't estimating anything. The code a <- c(1, 0.6, 0.8) x <- rnorm(2) y <- crossprod(c(1, x), a) generates 'y'; but nothing is estimated. 'y' is going to have the same variance regardless of how many observatio

Re: [R] Odp: Data frame modification

2010-07-28 Thread Petr PIKAL
Hi why do you insist on loops. R is not C. If you want to use loops use C or similar programming languages. It is almost always better to apply whole object approach. Kind and clever people already programmed it (sometimes in C ). x<-rnorm(20) x[c(10,12,13,17)]<-NA x [1] -1.12423790 0.80641

Re: [R] Help with specifiying random effects in lmer - psychology experiment

2010-07-28 Thread nuala brady
Apologies, not sure why that was so garbled first time! - posting again more simply I want to use lmer for analysis of response time (& error data) from a reaction time experiment, and have a question about specifying the structure of random effects in the model. I am using a repeated meas

Re: [R] Using jpeg() without X

2010-07-28 Thread Karl-Dieter Crisman
Something similar happens with png(), I think, if one doesn't compile against certain X headers (for instance on Mac). Is there any resource which describes exactly which graphics devices would work without *any* special options on any OS? I've only been able to find pdf() so far (the second pos

[R] Fwd: How to point a column of dataframe by a "character"

2010-07-28 Thread David Winsemius
(Forgot to copy the list.) Begin forwarded message: From: David Winsemius Date: July 28, 2010 7:44:38 AM EDT To: Tony Subject: Re: [R] How to point a column of dataframe by a "character" On Jul 28, 2010, at 5:35 AM, Tony wrote: Hello, Here is a dilemma I am having for a long time. But, I

Re: [R] [Rd] R CMD build wiped my computer

2010-07-28 Thread Kevin Wright
I doubt if it will be helpful for Jarrod, but when I hear stories like this, I remember reading about a similar event and a remarkable partial recovery here: http://www.justpasha.org/folk/rm.html Kevin On Wed, Jul 28, 2010 at 8:04 AM, Marc Schwartz wrote:

  1   2   >