[R] XLConnect readWorksheet comma decimal sign

2013-11-25 Thread Knut Krueger
Hi to all, how can I read exel files where the decimal sign is comma instead dot. I get the data as ascii and when converting 3,5 with as.numeric the 3,5 will be converted to NA Kind Regards Knut __ R-help@r-project.org mailing list

[R] convert data frame: two variables into _one_ binary variable

2013-11-25 Thread Liviu Andronic
Dear all, I am trying to convert the following data frame into a format more useful for me: library(HSAUR2, lib.loc=C:/Program Files/R/R-3.0.2/library) Loading required package: tools Attaching package: ‘HSAUR2’ The following object is masked _by_ ‘.GlobalEnv’: womensrole

Re: [R] XLConnect readWorksheet comma decimal sign

2013-11-25 Thread PIKAL Petr
Hi Either change comma to dot in Excel (but sometimes Excel is rather reluctant to accept such changes). Or change commaa to dot in R which probably can be easily done by gsub command Or read data with option dec=,. I do not know XLConnect but in read.table it is optional parameter and maybe

Re: [R] testing for bimodal and for dip in between modes in R

2013-11-25 Thread Simon Zehnder
Testing for bimodality is rather testing for unimodality. Hartigan and Hartigan (1985) presented the Dip-Test which is implemented in the R package DipTest with a much better approximation of the test distribution. If the test statistic is too high unimodality is rejected. To estimate the dip

Re: [R] XLConnect readWorksheet comma decimal sign

2013-11-25 Thread Knut Krueger
Am 25.11.2013 11:29, schrieb PIKAL Petr: Hi Either change comma to dot in Excel (but sometimes Excel is rather reluctant to accept such changes). Or change commaa to dot in R which probably can be easily done by gsub command Or read data with option dec=,. I do not know XLConnect but in

Re: [R] Should there be an R-beginners list?

2013-11-25 Thread Arman Eshaghi
I do not agree with a separate beginner's list. But I do stand with moving to stackoverflow, mainly because of the easier google search than current mailing list. It could make it more accessible. On Mon, Nov 25, 2013 at 4:07 AM, John Sorkin jsor...@grecc.umaryland.eduwrote: Mailing list vs.

[R] Independent variable dependent on offset in GLMM

2013-11-25 Thread Jonas Josefsson
Hi! I’m running glmer (lme4) models with biodiversity data and I’m having trouble with understanding/finding information on how the offset() option is implemented. Explicitly, I’m wondering if the offset is only implemented on the dependent variable (as I think it is), or does it also affect

[R] lmer specification for random effects: contradictory reults

2013-11-25 Thread Benedetta Cesqui
Hi All, I was wondering if someone could help me to solve this issue with lmer. In order to understand the best mixed effects model to fit my data, I compared the following options according to the procedures specified in many papers (i.e. Baayen

Re: [R] Should there be an R-beginners list?

2013-11-25 Thread Barry Rowlingson
If you want a vision of an R-beginners list, it is a boot stamping ITS IN THE DOCUMENTATION into a newbies face - forever. http://www.brainyquote.com/quotes/quotes/g/georgeorwe159438.html slight exaggeration perhaps, but most R-beginners would benefit from reading a bit more documentation and

[R] Structural break test Andrews (2003)

2013-11-25 Thread Geoffrey Smith
Dear Friends, I am looking for an R version of the structural break test in Andrews (2003). The excellent strucchange package does not include this test (yet?). Is this test available in another package? If not, might there already be a function written to do this test? Thank you very much.

Re: [R] Should there be an R-beginners list?

2013-11-25 Thread Hadley Wickham
I do not see how it can be illegal to download and duplicate the posts, since all the content is licensed under CC BY-SA. I might have missed something there: http://stackexchange.com/legal If that is really the case, I think I will have to reconsider if I should use it any more. I'm not a

Re: [R] lmer specification for random effects: contradictory reults

2013-11-25 Thread ONKELINX, Thierry
Dear Benedetta, I think you might want (1+T+Z|subject) as random effects rather than (1+T|subject) + (1 + Z|subject). The latter has two random intercepts per subject: a recipe for disaster. Follow-up posts should only go to the mixed models mailing list which I'm cc'ing. Best regards, ir.

Re: [R] Should there be an R-beginners list?

2013-11-25 Thread PIKAL Petr
Hi I doubt if people start to search answers if they often do not search them in help pages and documentation provided. I must agree with Duncan that if Stackoverflow was far more better than this help list most people would seek advice there then here. Is there any evidence in decreasing

Re: [R] Structural break test Andrews (2003)

2013-11-25 Thread Achim Zeileis
On Mon, 25 Nov 2013, Geoffrey Smith wrote: Dear Friends, I am looking for an R version of the structural break test in Andrews (2003). The excellent strucchange package does not include this test (yet?). Is this test available in another package? If not, might there already be a function

Re: [R] Plotting multiple trends on one graph

2013-11-25 Thread John Kane
I am not sure I have grasped what you want but have a look at this using ggplot2 with Duncan's modified data set (with his factor and striptime commands executed) You will probably need to install ggplot2 : install.packages(ggplot2) library(ggplot2) ggplot(dat, aes(gspd_mps, Station )) +

[R] Aggregating spatial data

2013-11-25 Thread IOANNA
Hello all, I have a data frame in the form: s-expand.grid(x=seq(1,100,by=1),y=seq(1,100,by=1)) w-data.frame(x=s$x,y=s$y,z1=rep(c(1,2,3,4),times=length(s$x/4)),z2=seq(1,le ngth(s$x),by=1)) The w$x and w$y represent the location of points and z1 and z2 attributes corresponding to these

Re: [R] Aggregating spatial data

2013-11-25 Thread David Carlson
Something like this? s - expand.grid(x=seq(1,100,by=1),y=seq(1,100,by=1)) w - data.frame(x=s$x,y=s$y,z1=rep(c(1,2,3,4),times=length(s$x/4)), + z2=seq(1,length(s$x),by=1)) w$EW - cut(w$x, breaks=c(.5, 50.5, 100.5), labels=c(West, East)) w$NS - cut(w$y, breaks=c(.5, 50.5, 100.5),

Re: [R] Should there be an R-beginners list?

2013-11-25 Thread Duncan Murdoch
On 25/11/2013 8:47 AM, Hadley Wickham wrote: I do not see how it can be illegal to download and duplicate the posts, since all the content is licensed under CC BY-SA. I might have missed something there: http://stackexchange.com/legal If that is really the case, I think I will have to

Re: [R] Should there be an R-beginners list?

2013-11-25 Thread Marc Schwartz
On Nov 25, 2013, at 7:56 AM, PIKAL Petr petr.pi...@precheza.cz wrote: Hi I doubt if people start to search answers if they often do not search them in help pages and documentation provided. I must agree with Duncan that if Stackoverflow was far more better than this help list most

Re: [R] Should there be an R-beginners list?

2013-11-25 Thread Hadley Wickham
Here's a similar plot for stackoverflow: http://data.stackexchange.com/stackoverflow/query/150130/r-questions-and-answers-per-year#graph and one broken down by month http://data.stackexchange.com/stackoverflow/query/150129/r-questions-and-answers-per-month#graph Hadley On Mon, Nov 25, 2013 at

Re: [R] Should there be an R-beginners list?

2013-11-25 Thread Barry Rowlingson
Joran (on StackOverflow chat, funnily enough) has just pointed us to this: http://www.win.tue.nl/~bvasiles/papers/cscw14.pdf How Social QA Sites are Changing Knowledge Sharing in Open Source Software Communities which includes a graph of postings to R-help and questions tagged '[r]' on

Re: [R] Should there be an R-beginners list?

2013-11-25 Thread Hadley Wickham
Oops, I misunderstood the database schema, and that only includes _questions_ tagged R, not the corresponding answers. Hadley On Mon, Nov 25, 2013 at 9:55 AM, Hadley Wickham h.wick...@gmail.com wrote: Here's a similar plot for stackoverflow:

Re: [R] Should there be an R-beginners list?

2013-11-25 Thread Paul Gilbert
On 13-11-25 06:00 AM, r-help-requ...@r-project.org wrote: Date: Sun, 24 Nov 2013 13:04:43 -0600 From: Yihui Xiex...@yihui.name To: Bert Guntergunter.ber...@gene.com Cc:r-help@r-project.org r-help@r-project.org Subject: Re: [R] Should there be an R-beginners list? Message-ID:

Re: [R] convert data frame: two variables into _one_ binary variable

2013-11-25 Thread Carl Witthoft
In R, as.logical() and other functions treat anything 0 as TRUE. Thus: Rgames foo-sample(0:5,10,rep=TRUE) Rgames foo [1] 0 5 1 0 1 5 2 5 4 5 Rgames as.logical(foo) [1] FALSE TRUE TRUE FALSE TRUE TRUE TRUE TRUE TRUE TRUE For your case, simply (womensrole$agreewomensrole$disagree)

Re: [R] Should there be an R-beginners list?

2013-11-25 Thread Carl Witthoft
A couple things: First, Beginners' lists never work. Beginners invariably can't read (cf. posting guidelines), so they will post to the main list anyway. Second, I see some people prefer to receive email-lists of the topics, and others prefer to work via a webbrowser interface. I'd have to say

Re: [R] convert data frame: two variables into _one_ binary variable

2013-11-25 Thread David Carlson
I think the OP was looking to expand the data frame so that each row was a single observation so that the first row becomes 6 rows, 4-TRUE and 2-FALSE. Something like this womensrole - HSAUR::womensrole step1 - reshape(womensrole, varying=c(agree, disagree), + v.names=Freq, timevar=Agree,

Re: [R] Speeding up code

2013-11-25 Thread MacQueen, Don
ditto to everything Jeff Newmiller said, but I'll take it a little further. I'm guessing that with df - data.frame(31790,31790) you thought you were creating something with 31790 rows and 31790 columns. You weren't. You were creating a data frame with one row and two columns:

Re: [R] convert data frame: two variables into _one_ binary variable

2013-11-25 Thread Thomas Stewart
I do not think Carl's solution answers your question. Try this: z - textConnection( education gender agree disagree sexe 1 0 Male 420 2 1 Male 200 3 2 Male 400 4 3 Male 630 5 4 Male

Re: [R] Plotting multiple trends on one graph

2013-11-25 Thread Law, Jason
Natalie, I'm assuming this is some kind of passive animal sampling? Instream PIT tags for fish? In that case, you can get what I think you want using ggplot2 and something like this: dat$TagID - as.factor(dat$TagID) dat$Station - as.factor(dat$Station) dat$Station2 - as.numeric(dat$Station)

[R] FW: Plotting multiple trends on one graph

2013-11-25 Thread Duncan Mackay
Sending again as apparently did not make the list Duncan -Original Message- From: Duncan Mackay [mailto:dulca...@bigpond.com] Sent: Monday, 25 November 2013 10:06 To: 'Natalie Houghton McNair' Cc: R Subject: RE: [R] Plotting multiple trends on one graph Hi Natalie Here is an option

[R] summary many regressions

2013-11-25 Thread Gary Dong
Dear R users, I have a large data set which includes data from 300 cities. I want to run a biviriate regression for each city and record the coefficient and the adjusted R square. For example, in the following, I have 10 cities represented by numbers from 1 to 10: x = cumsum(c(0, runif(999, -1,

Re: [R] convert data frame: two variables into _one_ binary variable

2013-11-25 Thread arun
Hi, One way would be: womensrole - read.table(text=education gender agree disagree sexe 1    0  Male    4    2    0 2    1  Male    2    0    0 3    2  Male    4    0    0 4    3  Male    6    3    0 5    4  Male    5    5    0 6    5  Male    13   

[R] Fetching data from MySQL via odbcConnect - Error in as.POSIXlt.character(x, tz, ...) :

2013-11-25 Thread Dirk Zimmermann
Hi! I am trying to retrieve data from a MySQL Database using RODBC with the commands odbcConnect and sqlFetch. There are different data files in the database and in some cases it works without any difficulties. Nevertheless I get an error with some data files. Since I'm not familiar with MySQL I

[R] R: lmer specification for random effects: contradictory reults

2013-11-25 Thread Benedetta Cesqui
Dear Thierry, thank you for the quick reply. I have only one question about the approach you proposed. As you suggested, imagine that the model we end up after the model selection procedure is: mod2.1 - lmer(dT_purs ~ T + Z + (1 +T+Z| subject), data =x, REML= FALSE) According to the common

[R] cut2 not binning interval endpoints correctly

2013-11-25 Thread Maximilian Butler
Hi everyone, I am attempting to bin a vector of numbers between 0 and 1 into intervals of 0.001 but many values at the endpoints of the intervals are getting binned into the wrong interval. For example, the first 3 rows are binned incorrectly here: library(Hmisc)

Re: [R] Hey guys

2013-11-25 Thread arun
Hi, Try: dat1- t(dat[,-1])  colnames(dat1) - dat[,1] covmat - cov(dat1) A.K. I'm running into this error and I'm not sure how fix it Error: is.numeric(x) || is.logical(x) is not TRUE This is my data frame geneExpression_Lab10.txt This is the homework instructions in case you want to see

[R] Durbin Watson Test Bound in R

2013-11-25 Thread MaXiaoyue
Hi, How could I use R to check Durbin Watson Test Bound? Best, Rebecca __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide

[R] specify breaks in divergent palette in RColorBrewer: was divergent colors around zero in levelplot()

2013-11-25 Thread Don McKenzie
Bert or anyone else familiar with RColorBrewer: Has anyone tried to accomplish with RColorBrewer what I asked about in my original post (below)? Here is an example cribbed from the levelplot() help examples x - seq(pi/4, 5 * pi, length.out = 100) y - seq(pi/4, 5 * pi, length.out = 100) r -

Re: [R] specify breaks in divergent palette in RColorBrewer: was divergent colors around zero in levelplot()

2013-11-25 Thread Don McKenzie
Never mind. Solved. “cuts” argument back in levelplot(). Duh. On Nov 25, 2013, at 4:27 PM, Don McKenzie d...@u.washington.edu wrote: Bert or anyone else familiar with RColorBrewer: Has anyone tried to accomplish with RColorBrewer what I asked about in my original post (below)? Here

Re: [R] cut2 not binning interval endpoints correctly

2013-11-25 Thread Jim Holtman
FAQ 7.31 Sent from my iPad On Nov 25, 2013, at 9:01, Maximilian Butler maximilian.but...@gmail.com wrote: Hi everyone, I am attempting to bin a vector of numbers between 0 and 1 into intervals of 0.001 but many values at the endpoints of the intervals are getting binned into the wrong

Re: [R] summary many regressions

2013-11-25 Thread David Winsemius
On Nov 25, 2013, at 3:35 PM, Gary Dong wrote: Dear R users, I have a large data set which includes data from 300 cities. I want to run a biviriate regression for each city and record the coefficient and the adjusted R square. For example, in the following, I have 10 cities represented

Re: [R] Durbin Watson Test Bound in R

2013-11-25 Thread Jeff Newmiller
I would first enter ?RSiteSearch at the R console, and learn how to sift through the 5000+ packages at CRAN for myself. Then I would enter RSiteSearch(Durbin) and then study the many options available. --- Jeff

Re: [R] specify breaks in divergent palette in RColorBrewer: was divergent colors around zero in levelplot()

2013-11-25 Thread Don McKenzie
In case anyone cares (?), here is a function to do what I was asking, which doesn’t use colorBrewer but could with some hacking. I’m sure it’s fragile, but it works with well behaved integers and zero in the middle, which was all I needed. The output is a palette that can be passed to

Re: [R] summary many regressions

2013-11-25 Thread arun
Hi, Try: res - do.call(rbind,lapply(split(data,data$city),function(z) {fit_city - lm(y~x,data=z);data.frame(City=unique(z$city),Coefficient=coef(fit_city)[2],Adjusted_R_square= summary(fit_city)$adj.r.squared)})) A.K. On Monday, November 25, 2013 6:37 PM, Gary Dong pdxgary...@gmail.com

Re: [R] specify breaks in divergent palette in RColorBrewer: was divergent colors around zero in levelplot()

2013-11-25 Thread C. Alina Cansler
Don, This looks helpful: https://stat.ethz.ch/pipermail/r-help/2011-March/272361.html Also, here is some code that I had, and tried to make applicable to your question: div.colors -colorRampPalette(c(blue, white, red )) x-seq(-1,12,1);x palette(div.colors(length(x))) y- rep(1,length(x))

Re: [R] Functional Programming patterns

2013-11-25 Thread mohan . radhakrishnan
Thanks. It is very useful. Mohan From: Hadley Wickham h.wick...@gmail.com To: mohan.radhakrish...@polarisft.com Cc: R-help r-help@r-project.org Date: 11/20/2013 06:46 PM Subject:Re: [R] Functional Programming patterns I have some notes on functional programming in R at

[R] Multinomial regression with repeated measures

2013-11-25 Thread Rameela Chandrasekhar
Hi R-users,   Does anyone know of a package in R that could be used to implement a multinomial regression model accounting for repeated measures? Thanks. [[alternative HTML version deleted]] __ R-help@r-project.org mailing list

Re: [R] Durbin Watson Test Bound in R

2013-11-25 Thread Achim Zeileis
On Mon, 25 Nov 2013, MaXiaoyue wrote: Hi, How could I use R to check Durbin Watson Test Bound? Although the bounds approach to the Durbin-Watson test is still described in many text books, there are actually various ways to get p-values for the Durbin-Watson statistic. dwtest() in lmtest

Re: [R] specify breaks in divergent palette in RColorBrewer: was divergent colors around zero in levelplot()

2013-11-25 Thread Achim Zeileis
On Mon, 25 Nov 2013, C. Alina Cansler wrote: Don, This looks helpful: https://stat.ethz.ch/pipermail/r-help/2011-March/272361.html Also, here is some code that I had, and tried to make applicable to your question: div.colors -colorRampPalette(c(blue, white, red )) x-seq(-1,12,1);x

[R] Bootstrap or subsampling using loop?

2013-11-25 Thread Suparna Mitra
Hello R experts, I am trying to do a job where I need to take random subsample from a data and then frequency count of that. Then the median or the frequency from say, 1000 replicates. Should I try this with subsample in loop or bootstrap? My data format is str(Data) 'data.frame': 155752 obs.