Re: [R] glmm: random term, overdispersion and comparisons

2015-10-02 Thread Thierry Onkelinx
Dear Ellen, You're using the Poisson distribution. There is no error (noise) term in a glmm with Poisson distribution. 1) The random part seems to be quite complicated given the sample size. (1|site) is probably sufficient. Note that your design is not nested but crossed. 2) Overdispersion is lik

Re: [R] tcltk table properties

2015-10-02 Thread peter dalgaard
On 02 Oct 2015, at 03:32 , Dan D wrote: > I have a tkwidget table (say, tbl1) that may be reconfigured at various times > depending on user input. Is there an easy way to later extract table > properties? Something like... > > nrow<-tkgetproperties(tbl1, rows) > More like tkcget(tbl1, rows=N

Re: [R] Counting number of rain

2015-10-02 Thread Duncan Murdoch
On 01/10/2015 11:29 PM, Rolf Turner wrote: > On 02/10/15 15:47, David Winsemius wrote: > > > >> On Oct 1, 2015, at 6:22 PM, Rolf Turner wrote: >>> >>> P.S. I have been unable to find a corresponding vector of the names >>> of the days of the week, although I have a very vague recollection >>> of

[R] Rayleigh Distribution

2015-10-02 Thread Jomy Jose
Is it possible to code in R to get Q-Q plot for Rayleigh distribution Jose [[alternative HTML version deleted]] __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the

Re: [R] problem of Mahalanobis distance matching using MatchIT

2015-10-02 Thread Giorgio Garziano
About the "distance = NA" issue, please see if this comment helps: https://lists.gking.harvard.edu/pipermail/matchit/2011-January/000174.html Furthermore, the Mahalanobis NA distance values are hard-wired in the code, file matchit.R: ## no distance for full mahalanobis matching if(fn1=="di

Re: [R] Help with improveProb function in Hmisc in R

2015-10-02 Thread Frank Harrell
Please note that Kirsten is cross-posting to stats.stackexchange.com creating extra work for everyone. -- Frank E Harrell Jr Professor and Chairman School of Medicine Department of *Biostatistics* *Vanderbi

Re: [R] Rayleigh Distribution

2015-10-02 Thread Giorgio Garziano
Similarly to what can be read on help(qqplot), however using Rayleigh distribution: library(VGAM) p <- ppoints(100) x <- qrayleigh(p) y <- rrayleigh(100) qqplot(x,y) qqline(y, distribution = function(p) qrayleigh(p), prob = c(0.1, 0.6), col = 2) -- GG [[alternative HTML version d

Re: [R] Rayleigh Distribution

2015-10-02 Thread Duncan Murdoch
On 02/10/2015 8:24 AM, Jomy Jose wrote: > Is it possible to code in R to get Q-Q plot for Rayleigh distribution If you have a quantile function for any distribution, you can make a QQ plot using ppoints() and the quantile function. The Rayleigh distribution isn't in base R, but Google says it is

Re: [R] Counting number of rain

2015-10-02 Thread David Winsemius
On Oct 2, 2015, at 2:33 AM, Duncan Murdoch wrote: > On 01/10/2015 11:29 PM, Rolf Turner wrote: >> On 02/10/15 15:47, David Winsemius wrote: >> >> >> >>> On Oct 1, 2015, at 6:22 PM, Rolf Turner wrote: P.S. I have been unable to find a corresponding vector of the names of the day

Re: [R] Rayleigh Distribution

2015-10-02 Thread Jomy Jose
Thank you,how to add correct reference line to this Rayleigh Q-Q plot. Jose [[alternative HTML version deleted]] __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read th

[R] Q-Q plot reference line

2015-10-02 Thread Jomy Jose
How can the reference line be drawn for various distributions like Rayleigh and Log logistic Q-Q plots. Thanks Jomy [[alternative HTML version deleted]] __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/

Re: [R] Q-Q plot reference line

2015-10-02 Thread Jeff Newmiller
You sound sad, lost and wandering about in the wilderness. Fortunately there is help for you in the Posting Guide (mentioned at the bottom of every email on this list), which advises you to keep replying on the same thread of emails as long as the topic has not changed, post using plain text rat

Re: [R] Regressing the residuals on the country dummies

2015-10-02 Thread Linus Holtermann
Hi, You have panel data or cross-sectional data? In the case you use cross-sectional data and "countries" are your observations (no repeated measure of them) and you regress the country-dummies on your residuals of the forgone regression, then there are as many regressors as observations. Conse

Re: [R] Rayleigh Distribution

2015-10-02 Thread Johannes Huesing
Jomy Jose [Fri, Oct 02, 2015 at 05:45:58PM CEST]: Thank you,how to add correct reference line to this Rayleigh Q-Q plot. ?abline -- Johannes Hüsing http://derwisch.wikidot.com Threema-ID: VHVJYH3H __ R-help@r-project.org mailing lis

Re: [R] Q-Q plot reference line

2015-10-02 Thread Rolf Turner
On 03/10/15 06:08, Jeff Newmiller wrote: Finally, one definition of insanity is trying to solve a problem by repeating the same strategy that did not work the last time... in this case, sending essentially the same ineffective email to the list multiple times in rapid succession. I hope you ge

Re: [R] Counting number of rain

2015-10-02 Thread Rolf Turner
On 03/10/15 04:42, David Winsemius wrote: On Oct 2, 2015, at 2:33 AM, Duncan Murdoch wrote: The zoo package replaces as.Date.numeric() with a function that assumes an origin of "1970-01-01". There may be other packages that also make a replacement like this. David appears to have one of t

[R] Adding a decimal point at the second mark in R

2015-10-02 Thread Jon1985
Using R, I have two columns (onset, duration) containing milliseconds in a data frame. I want to add a decimal point at the second mark for all data points in both columns i.e. turning 1541ms into 1.541s AND turning 638ms into 0.638s. Could someone tell me how to do this? Thanks Jon -- View t

[R] Package compilation errors

2015-10-02 Thread Noel Hunt
I have a few packages for which compilation has failed. The errors are simple compile-time errors I could fix if the configured source was left untouched and not removed. I have found out how to stop source being removed but I want the whole set of extracted, configured files left, so I can fix th

[R] John M. Chambers Statistical Software Award

2015-10-02 Thread Patrick Breheny
John M. Chambers Statistical Software Award 2016 The Statistical Computing Section of the American Statistical Association announces the competition for the John M. Chambers Statistical Software Award. In 1998 the Association for Computing Machinery presented its Software System Award to John

Re: [R] Adding a decimal point at the second mark in R

2015-10-02 Thread Rolf Turner
On 03/10/15 05:26, Jon1985 wrote: Using R, I have two columns (onset, duration) containing milliseconds in a data frame. I want to add a decimal point at the second mark for all data points in both columns i.e. turning 1541ms into 1.541s AND turning 638ms into 0.638s. Could someone tell me how

[R] "split" and loop functions

2015-10-02 Thread Greg Damico
Hello, I wonder if you might be able to help me. I'm enrolled in an R programming course through Coursera. I've done well so far--though it's been challenging!--but I'm having trouble understanding exactly how "split" and the loop functions (like "lapply") work. I keep getting an error that say

[R] scatter3d

2015-10-02 Thread Jeff Tostenrude
I am using scatter3d in R Commander to plot a group of regression surfaces. However, I only want to display the surfaces. How do I remove the points? [[alternative HTML version deleted]] __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and

Re: [R] scatter3d

2015-10-02 Thread Fox, John
Dear Jeff, I'm tempted to say that a scatterplot without points is an oxymoron, but that wouldn't be very helpful. Actually, the scatter3d() function used by the Rcmdr is in the car package. The Rcmdr 3D scatterplot dialog doesn't provide for suppressing the points, but if you add radius=rep(0

Re: [R] "split" and loop functions

2015-10-02 Thread William Dunlap
You should really be asking Coursera for help with its course. However you can figure this out by breaking down your commands into small steps and seeing what each step gives. This is an advantage of an interactive system like R. E.g., suppose you start with > s1 <- data.frame(state=c("Rhode I

Re: [R] scatter3d

2015-10-02 Thread Jeff Newmiller
Wouldn't you just stop using scatter3d? --- Jeff NewmillerThe . . Go Live... DCN:Basics: ##.#. ##.#. Live Go... Live: OO#.. Dead:

Re: [R] Package compilation errors

2015-10-02 Thread Jeff Newmiller
No, don't get in the way of install.packages(). Download the source versions of the packages and unpack them yourself and fix the source files and update the DESCRIPTION files. Then follow the instructions in the "Writing R Extensions" document at the command line (I.e. R CMD check, R CMD build,

Re: [R] scatter3d

2015-10-02 Thread Jeff Tostenrude
Thank you for the suggestion. Yes, it does seem to be a bit counter-intuitive, but that is the output I am being asked to produce. I have only been using R for a few weeks, so there is probably a better way to do it. Anyway, your suggestion did work, but only up to 337 data points. I am dealing wi