Re: [R] R studio and git

2015-09-25 Thread David Winsemius
This is not the RStudio support list. That function is offered on a web-hosted support list. -- David On Sep 25, 2015, at 6:19 PM, Glenn Schultz wrote: > Hello all, > > Maybe not the right place. I have a project that is under git version > control in R studio. I just updated to recent OS

[R] flexmix - concomitant model and significance of variables

2015-09-25 Thread Damien Jourdain
Dear All, I am new to this forum and to flexmix. I am using flexmix to make a cluster analysis (Model Based). The data for the clustering are all continuous (although all between 0 and 1). I also want to see the correlation between found clusters and some socio economic variables, so I am using

[R] Analysis of causal relations between rare (categorical) events

2015-09-25 Thread markshanks
Hi, I have only a relatively basic background in statistics (e.g., anova, regression), and the books on R I have read so far have focused on relatively common statistical analyses (e.g., outlier analysis, trend forecasting) and haven't helped me with the data problem I am facing. In short, imagin

[R] R studio and git

2015-09-25 Thread Glenn Schultz
Hello all, Maybe not the right place.  I have a project that is under git version control in R studio.  I just updated to recent OSX build and now the git tab in R studio is not available.  It seems I can still use version control with source tree but it is strange that I lost git version on t

Re: [R] How to get significance codes after Kruskal Wallis test

2015-09-25 Thread Kristina Wolf
Perhaps look into the function friedman.test.with.post.hoc() There is more information here: http://www.r-statistics.com/wp-content/uploads/2010/02/Friedman-Test-with-Post-Hoc.r.txt Note, this does not handle NA's though, and technically it is for blocked designs, but maybe it will lead you somewh

[R] How to calculate standard error of estimate (S) for my non-linear regression model?

2015-09-25 Thread Michael Eisenring
Hi all, I am looking for something that indicates the goodness of fit for my non linear regression model (since R2 is not very reliable). I read that the standard error of estimate (also known as standard error of the regression) is a good alternative. The standard error of estimate is descri

[R] How to get significance codes after Kruskal Wallis test

2015-09-25 Thread Michael Eisenring
Is there a way to get significance codes after a pairwise comparisons to a Kruskall wallis test? With significance codes I mean letter codes (a, b,c) that are assigned to treatments to indicate where differences are significant. With a traditional anova such a test can be performed using HSD.test

Re: [R] backsolve, chol, Matrix, and SparseM

2015-09-25 Thread Benjamin Tyner
Hi Martin, Thanks for the remarks and examples, and for confirming that I was indeed barking up the wrong tree with SparseM. A. I assume that is a typo and you meant to say, no need for backsolve(). B. Absolutely; however, in this case I am taking advantage of quadprog::solve.QP(..., factorized =

Re: [R] Multibyte strings

2015-09-25 Thread David Winsemius
On Sep 25, 2015, at 2:23 PM, Dennis Fisher wrote: > R 3.2.0 > OS X > > Colleagues, > > Earlier today, I initiated a series of emails regarding SASxport (which was > removed from CRAN). David Winsemius proposed downloading the source code and > installing with the following command: >

[R] Multibyte strings

2015-09-25 Thread Dennis Fisher
R 3.2.0 OS X Colleagues, Earlier today, I initiated a series of emails regarding SASxport (which was removed from CRAN). David Winsemius proposed downloading the source code and installing with the following command: install.packages('~/Downloads/SASxport_1.5.0.tar.gz', repos = NULL ,

Re: [R] Sampling the Distance Matrix

2015-09-25 Thread David Winsemius
On Sep 25, 2015, at 12:54 PM, Lorenzo Isella wrote: > Apologies for not letting this thread rest in peace. > The small script > > # > set.seed(1234) > > x <- rnorm(20) > y <- rnorm(20) > > > goodcls <- apply(mtxcomb , 2, function(idx) al

Re: [R] Accessing defunct package

2015-09-25 Thread David Winsemius
On Sep 25, 2015, at 11:23 AM, Dennis Fisher wrote: > David > > Thanks for suggesting this. Three issues of note: > > 1. I too work on a Mac. When I downloaded the last archived version (which > has a .gz extension), OS X automatically unzipped the file and removed the > .gz extension. I w

Re: [R] Sampling the Distance Matrix

2015-09-25 Thread Lorenzo Isella
Apologies for not letting this thread rest in peace. The small script # set.seed(1234) x <- rnorm(20) y <- rnorm(20) goodcls <- apply(mtxcomb , 2, function(idx) all( dist( cbind( x[idx], y[idx]) ) > 0.9)) mycomb <- mtxcomb [ , goodcls] #

Re: [R] Randomness tests

2015-09-25 Thread Giorgio Garziano
By "interested in any kind of deviation from randomness", I mean that I would like to apply all "randtests" R package randomness tests, if they give reliable results for {-1, 1} sequences. [[alternative HTML version deleted]] __ R-help@r-projec

Re: [R] Sampling the Distance Matrix

2015-09-25 Thread Lorenzo Isella
Absolutely right! Thanks to both David for their help. Cheers Lorenzo On Fri, Sep 25, 2015 at 01:54:54PM +, David L Carlson wrote: You defined x and y in your original email as: x<-rnorm(20) y<-rnorm(20) mm<-as.matrix(cbind(x,y)) dst<-(dist(mm)) - D

Re: [R] Accessing defunct package

2015-09-25 Thread Sarah Goslee
Not an R question, but tangentially relevant: On Fri, Sep 25, 2015 at 2:23 PM, Dennis Fisher wrote: > David > > Thanks for suggesting this. Three issues of note: > > 1. I too work on a Mac. When I downloaded the last archived version (which > has a .gz extension), OS X automatically unzipped

Re: [R] Accessing defunct package

2015-09-25 Thread Dennis Fisher
David Thanks for suggesting this. Three issues of note: 1. I too work on a Mac. When I downloaded the last archived version (which has a .gz extension), OS X automatically unzipped the file and removed the .gz extension. I was able to gzip the file so that I could execute your exact comman

Re: [R] Randomness tests

2015-09-25 Thread Bert Gunter
Yes, but I'll just note offlist that a test for "any deviation from randomness" is mathematically impossibile (ask on stackexchange for why if you like), so you may wish to think about the issue more carefully -- or consult a local statistician for advice. Cheers, Bert Bert Gunter "Data is not in

Re: [R] Quantmod several indicators

2015-09-25 Thread giorgio.garzi...@tin.it
Hi, Both following code examples plot bollinger band over the ticker main plotand Williams Power below the main plot. 1. chartSeries(YHOO,theme="white",TA = c(addBBands(200,2), addWPR(n=300))) 2. chartSeries(IB,theme="white",TA = c(addBBands(200,2))) addWPR(n=300) See: http://around-r.blogspot.i

Re: [R] Accessing defunct package

2015-09-25 Thread Ben Bolker
Dennis Fisher plessthan.com> writes: > Colleagues, > > In the past, I used a package: > SASxport > to output files to SAS’s XPT format. This was useful because FDA requests that data be submitted in that > format (even though they typically must reconvert to some other format before the

[R] ARCHETYPAL ANALYSIS IN R

2015-09-25 Thread euthymios kasvikis
To whom it may concern I am a postgraduate student in Aristotle University of Thessaloniki and a recent user of R also writing my thesis statement on a Sport Analytics tool based on Pr.Eugster's work on Archetypal Analysis. The problem is that after importing the data-set i want to work with, i can

Re: [R] Quantifying widths of polygons

2015-09-25 Thread David Winsemius
On Sep 25, 2015, at 7:20 AM, Alexander Shenkin wrote: > Hello all, > > I am working with data on tree crowns, and this data describes points > (verticies) around the polyhedron of the crown volume (think of the crown as > a single volume with vertices and faces). I can calculate maximum z dis

Re: [R] Accessing defunct package

2015-09-25 Thread David Winsemius
I think it might be even simpler at least for now. The error checking done by CRAN can be more rigorous than that done when an installation is done locally. I don't see a report in the current package checks listing of what error was identified, but experimentation is always an option. When I d

Re: [R] Randomness tests

2015-09-25 Thread Giorgio Garziano
Good suggestion, thanks. -Original Message- From: Jeff Newmiller [mailto:jdnew...@dcn.davis.ca.us] Sent: venerdì 25 settembre 2015 18:49 To: Giorgio Garziano; r-help@r-project.org Subject: Re: [R] Randomness tests You are way off topic for this list. Perhaps stats.stackexchange.com would

Re: [R] Randomness tests

2015-09-25 Thread Jeff Newmiller
You are way off topic for this list. Perhaps stats.stackexchange.com would be a better place to ask such a question. --- Jeff NewmillerThe . . Go Live... DCN:Basics: ##.#.

Re: [R] Randomness tests

2015-09-25 Thread Giorgio Garziano
I am interested in any kind of deviation from randomness. I would like to know if the fact that a time series can take values only from the set {-1, 1} restricts the type of randomness tests that can be done. -- Giorgio Garziano [[alternative HTML version deleted]] _

Re: [R] Quantmod several indicators

2015-09-25 Thread Giorgio Garziano
Hi, Both following code examples plot Bollinger bands over the ticker main plot and William's Percent below the main plot. 1. chartSeries(YHOO,theme="white",TA = c(addBBands(200,2), addWPR(n=300))) 2. chartSeries(IB,theme="white",TA = c(addBBands(200,2))) addWPR(n=300) In general, if you lik

Re: [R] Accessing defunct package

2015-09-25 Thread Jeff Newmiller
Obtain the source package and fix it? Most errors are relatively minor adjustments that just require reading the updated "Writing R Extensions" document to figure out. You might be unlucky, but I think the odds are in your favor. --

Re: [R] Accessing defunct package

2015-09-25 Thread Marc Schwartz
> On Sep 25, 2015, at 9:25 AM, Dennis Fisher wrote: > > R 3.2.0 > OS X > > Colleagues, > > In the past, I used a package: > SASxport > to output files to SAS’s XPT format. This was useful because FDA requests > that data be submitted in that format (even though they typically must > r

Re: [R] Quantifying widths of polygons

2015-09-25 Thread Bert Gunter
This is not the right list for such substantive questions. I suggest that you check out the "spatial" and perhaps the "Environmetrics" Task views on CRAN and/or post to a statistical site like stats.stackexchange.com or the R-sig-ecology mailing list instead. Cheers, Bert Bert Gunter "Data is

[R] Accessing defunct package

2015-09-25 Thread Dennis Fisher
R 3.2.0 OS X Colleagues, In the past, I used a package: SASxport to output files to SAS’s XPT format. This was useful because FDA requests that data be submitted in that format (even though they typically must reconvert to some other format before the data are used). It appears that t

[R] Quantifying widths of polygons

2015-09-25 Thread Alexander Shenkin
Hello all, I am working with data on tree crowns, and this data describes points (verticies) around the polyhedron of the crown volume (think of the crown as a single volume with vertices and faces). I can calculate maximum z distance between any 2 points (maximum depth) and maximum x/y dist

Re: [R] There is an error when I performed chisquare and postanova test...:(

2015-09-25 Thread John Kane
As Michael suggests carefully read the help file (type ?chisq.test to bring up the page) and match all the elements of your command with the elements listed on the help page. You do not need to use attach(chicken) given the code you are using and using attach(data) is generally considered bad

Re: [R] Sampling the Distance Matrix

2015-09-25 Thread David L Carlson
You defined x and y in your original email as: > x<-rnorm(20) > y<-rnorm(20) > > mm<-as.matrix(cbind(x,y)) > > dst<-(dist(mm)) - David L Carlson Department of Anthropology Texas A&M University College Station, TX 77840-4352 -Original Message- From: Da

Re: [R] There is an error when I performed chisquare and postanova test...:(

2015-09-25 Thread Michael Dewey
Let me repeat my advice Read the documentation for chisq.test and see if it has the parameters you think it has. On 25/09/2015 14:12, massmatics wrote: Yes, so far, we only learned how to use chisq.test when performed at 0.05 significance level But not when it is at 0.01 -- View this messag

Re: [R] There is an error when I performed chisquare and postanova test...:(

2015-09-25 Thread massmatics
Yes, so far, we only learned how to use chisq.test when performed at 0.05 significance level But not when it is at 0.01 -- View this message in context: http://r.789695.n4.nabble.com/There-is-an-error-when-I-performed-chisquare-and-postanova-test-tp4712770p4712772.html Sent from the R help mail

Re: [R] There is an error when I performed chisquare and postanova test...:(

2015-09-25 Thread Michael Dewey
Dear massmatics The list has a no homework policy but in this case I feel I can depart from that. Have you tried reading the documentation for chisq.test? On 25/09/2015 12:44, massmatics wrote: I am studying for a quiz and while I was solving two problems, an error occured. a) So the problem

[R] There is an error when I performed chisquare and postanova test...:(

2015-09-25 Thread massmatics
I am studying for a quiz and while I was solving two problems, an error occured. a) So the problem is we want to test if the block program makes a difference in retention at α = 1% (which is the significance level) This is the code I used to perform chisquare test. My solution: blockprogram <- matr

Re: [R] backsolve, chol, Matrix, and SparseM

2015-09-25 Thread Martin Maechler
Dear Ben, > Benjamin Tyner > on Thu, 24 Sep 2015 13:47:58 -0400 writes: > Hi I have some code which does (on a symmetric matrix 'x') > backsolve(chol(x), diag(nrow(x))) > and I am wondering what is the recommended way to > accomplish this when x is also sparse (