Re: [R] constrast in edge R

2015-01-14 Thread Jeff Newmiller
You seem to be writing in a jargon that is not typical on this list, and I cannot identify what package you are using. I suspect that you need to be asking your question in the Bioconductor support area. --- Jeff Newmiller

Re: [R] (two way) ANCOVA and subsequent multiple comparison tests

2015-01-14 Thread Richard M. Heiberger
You found the HH package. That is a good start. Look at the ?MMC help page, and specifically at the entire maiz example. For an example with two factors and a covariate and with multiple comparisons look at the apple example in file system.file("scripts/hh2", package="HH") The example is in c

Re: [R] merging pre-sorted data frames

2015-01-14 Thread Alex Fun
package dplyr's full_join, left_join, right_join, inner_join are also comparable in speed to data table. The syntax is also more like merge's. On Thu, Jan 15, 2015 at 2:17 PM, Mike Miller wrote: > Thanks, Jeff. You really know the packages. I search and I guess I > didn't use the right terms.

Re: [R] Install R on linux Mint

2015-01-14 Thread Jeroen Ooms
On Wed, Jan 14, 2015 at 7:40 PM, Jeroen Ooms wrote: > Try 'precise' instead of 'utopic' I meant 'trusty'. The version should match the underlying ubuntu: http://www.linuxmint.com/oldreleases.php __ R-help@r-project.org mailing list -- To UNSUBSCRIBE an

Re: [R] Install R on linux Mint

2015-01-14 Thread Jeroen Ooms
Try 'precise' instead of 'utopic' On Wed, Jan 14, 2015 at 7:31 PM, John Sorkin wrote: > I am trying to install R on Linux mint 17.1. I followed the instructions > found on CRAN and got a messages about unmet dependencies. I detail below the > steps I took: > > I added deb http://lib.stat.cmu.ed

[R] Install R on linux Mint

2015-01-14 Thread John Sorkin
I am trying to install R on Linux mint 17.1. I followed the instructions found on CRAN and got a messages about unmet dependencies. I detail below the steps I took: I added deb http://lib.stat.cmu.edu/R/CRAN/bin/linux/ubuntu utopic/ to /etc/apt/sources.list As can be seen below: john-OptiPlex-G

Re: [R] merging pre-sorted data frames

2015-01-14 Thread Mike Miller
Thanks, Jeff. You really know the packages. I search and I guess I didn't use the right terms. That package seems to do exactly what I wanted. Mike On Tue, 13 Jan 2015, Jeff Newmiller wrote: On Tue, 13 Jan 2015, Mike Miller wrote: I have many pairs of data frames each with about 15 mill

[R] Using portfolio.optim in tseries (portfolio optimization question)

2015-01-14 Thread Ernie Stokely
I apologize if this question is posted in the wrong place. I am using portfolio.optim to run an optimization on a stock portfolio. As I understand modern portfolio theory, to run a mean-variance optimization of the allocation for a portfolio, you must specify an expected return. The examples at

Re: [R] Probably dumb question about lists

2015-01-14 Thread Rolf Turner
I will add a "[.tile.list" method in the next release of deldir so that the appropriate attributes are preserved. This won't happen for a month or so though; I have some other pressures on me at the moment. cheers, Rolf Turner On 15/01/15 04:52, Raphael Päbst wrote: Thanks Boris! I know th

Re: [R] Help installing packages with dependencies for R, behind corporate firewall

2015-01-14 Thread Prof Brian Ripley
On 14/01/2015 21:37, jose.nunez-zul...@barclays.com wrote: Hello R-users, I have no practical experience with the R language itself but I've been tasked to install it behind a corporate firewall. Basic installation seems sane but when my user tries to install a custom library like this: insta

[R] (two way) ANCOVA and subsequent multiple comparison tests

2015-01-14 Thread lindsay hanford
Hello R Community! I am an intermediate-level R user and I am trying to figure out how program a two-way (Group and MCCB score) ANCOVA analysis and subsequent post-hoc analysis. My factors: Group (2 levels) Score (2 levels) Covariates: Age (continuous) ICV (continuous) Sex (M/F) Y: GMV(continuous)

[R] Help installing packages with dependencies for R, behind corporate firewall

2015-01-14 Thread jose.nunez-zuleta
Hello R-users, I have no practical experience with the R language itself but I've been tasked to install it behind a corporate firewall. Basic installation seems sane but when my user tries to install a custom library like this: install.packages("ggplot2") Installing package into '/home/myuser/

Re: [R] Approximation of a function from R^2 to R

2015-01-14 Thread Greg Snow
Yes, there are several. Which is best and which subset to suggest depends on what you are trying to do, what your inputs look like (do you have the function, but want a simpler approximation? or do you have observations/datapoints?) If you can give us more detail about what you have to work with

[R] constrast in edge R

2015-01-14 Thread Noha Osman
Hey everyone, I have data frame called subtest as following: RNA.LATER.MEN2B_S5 RNA.LATER.ROSA_S4 RNA.MEN2B.1_S2 RNA.MEN2B.2_S3 RNA.ROSA_S1 1 13707 13866 12193 12671 10178 2 0 0 0 0 1 3 7165 5002 1256 1341 2087 6 8537 16679 9042 9620 19168 10 19438 25234 15563 16419 16582 16 3 3 11 3 5 I would li

[R] Approximation of a function from R^2 to R

2015-01-14 Thread rala
Hello, I am a beginner and I wanted to know if there is a way to approximate a function that has 2 input variables and one output. -- View this message in context: http://r.789695.n4.nabble.com/Approximation-of-a-function-from-R-2-to-R-tp4701805.html Sent from the R help mailing list archive a

[R] stringr::str_split_fixed query

2015-01-14 Thread David Barron
I'm puzzled as to why I get this behaviour with str_split_fixed in the stringr package. > stringr::str_split_fixed('ab','',2) [,1] [,2] [1,] "" "ab" > stringr::str_split_fixed('ab','',3) [,1] [,2] [,3] [1,] "" "a" "b" In the first example, I was expecting to get [,1] [,2] [1,

Re: [R] regular expression question

2015-01-14 Thread MacQueen, Don
Good point, John. Illustrates the danger of assuming there are no "perverse cases". -Don -- Don MacQueen Lawrence Livermore National Laboratory 7000 East Ave., L-627 Livermore, CA 94550 925-423-1062 From: John McKown mailto:john.archie.mck...@gmail.com>> Date: Wednesday, January 14, 2015 at 8

Re: [R] regular expression question

2015-01-14 Thread John McKown
On Wed, Jan 14, 2015 at 10:03 AM, MacQueen, Don wrote: > I know you already have a couple of solutions, but I would like to mention > that it can be done in two steps with very simple regular expressions. I > would have done: > > s <- c("lngimbintrhofixed","lngimbnointnorhofixed","test", >

Re: [R] regular expression question

2015-01-14 Thread MacQueen, Don
I know you already have a couple of solutions, but I would like to mention that it can be done in two steps with very simple regular expressions. I would have done: s <- c("lngimbintrhofixed","lngimbnointnorhofixed","test", 'rhofixedtest','norhofixedtest') res <- gsub('norhofixed$', '',s) r

Re: [R] ""error 1 in extracting from zip file"

2015-01-14 Thread John McKown
On Wed, Jan 14, 2015 at 4:32 AM, Emilio Gianicolo wrote: > Hi, > > I've some troubles in running this program > http://rpubs.com/adam_dennett/8955. > > I create this object: > temp <- tempfile(fileext = ".zip") > > But, after downloading this file > download.file(" > http://epp.eurostat.ec.europa

Re: [R] Rstudio and GIT

2015-01-14 Thread Karim Mezhoud
Hi, I think you need to open your project after In RStudio, go to project options, then to the Git/SVN panel. Change the “Version control system” from “None” to “Git”: http://r-pkgs.had.co.nz/git.html#git-status Karim Ô__ c/ /'_;kmezhoud (*) \(*) ⴽⴰⵔⵉⵎ ⵎⴻⵣⵀⵓⴷ http://bioinformatics.tn/

Re: [R] Probably dumb question about lists

2015-01-14 Thread Raphael Päbst
Thanks Boris! I know that I am just leaving out some tiles in the plot, but I specifically want to drop those tiles for the plot and not calculate a tesselation without these points. So I'll go with this solution and hope it makes sense in the end. Thanks again! Raphael On 1/14/15, Boris Steipe

[R] error 1 in extracting from zip file

2015-01-14 Thread Emilio Antonio Luca Gianicolo
Hi, I've some troubles in running this program http://rpubs.com/adam_dennett/8955. I create this object: temp <- tempfile(fileext = ".zip") But, after downloading this file download.file("http://epp.eurostat.ec.europa.eu/cache/GISCO/geodatafiles/NUTS_2010_60M_SH.zip";, temp) when I try to unz

Re: [R] Psych package: why am I receiving "NA" for many of the factor scores?

2015-01-14 Thread William Revelle
Dear Elizabeth, Factor scores in the fa function are found by multiplying the standardized data by the factor weights using matrix multiplication. This will give scores only for subjects with complete data. However, if you want, you can create them yourself by standardizing your data and then

[R] sparse matrix from vector outer product

2015-01-14 Thread Philipp A.
Hi, creating a matrix from two vectors a, b by multiplying each combination can be done e.g. via a %*% t(b) or via outer(a, b) # default for third argument is '*' But this yields a normal matrix. Is there an efficient way to create sparse matrices (from the Matrix package) like that? Right

[R] Rstudio and GIT

2015-01-14 Thread AURORA GONZALEZ VIDAL
Dear R users, I have a  very specific question. I want to know how to create a local git repository from an exisitng file (with some documents inside) just like we do when typing git init but from Rstudio. I tried selecting FIle-->New Project-->Existing Directory--> and I select the file but I

[R] ""error 1 in extracting from zip file"

2015-01-14 Thread Emilio Gianicolo
Hi, I've some troubles in running this program http://rpubs.com/adam_dennett/8955. I create this object: temp <- tempfile(fileext = ".zip") But, after downloading this file download.file("http://epp.eurostat.ec.europa.eu/cache/GISCO/geodatafiles/NUTS_2010_60M_SH.zip";, temp) when I try to un

Re: [R] Probably dumb question about lists

2015-01-14 Thread Boris Steipe
plot.tile.list() expects its argument to be of class tile.list, and to have an attribute "rw", both of which are not conserved after subsetting. You can do ... if (!require(deldir)) { install.packages("deldir") library(deldir) } x <- rnorm(10) y <- rnorm(10) del <- deldir(x, y) tl <- tile.li

Re: [R] Probably dumb question about lists

2015-01-14 Thread Bert Gunter
1. Please in future specify the package (deldir here) that contains the functions you refer to. 2. **Always** first try ?str before posting queries like this, as this will often reveal the problem. str(tl[1:800]) 3. I would **guess** (ergo could well be wrong) that "[" is not preserving the clas

Re: [R] Probably dumb question about lists

2015-01-14 Thread PIKAL Petr
Hi I do not know what is deldir, tile.list and what is structure of tl and do not have time to search. what is result of tl[1:20]? Do you expect such result? maybe plot(lapply(tl, "[", 1:800)) shall work. Cheers Petr > -Original Message- > From: R-help [mailto:r-help-boun...@r-proje

Re: [R] using poly() to predict

2015-01-14 Thread Prof Brian Ripley
On 14/01/2015 14:20, Stanislav Aggerwal wrote: This method of finding yhat as x %*% b works when I use raw polynomials: x<-1:8 y<- 1+ 1*x + .5*x^2 fit<-lm(y~poly(x,2,raw=T)) b<-coef(fit) xfit<-seq(min(x),max(x),length=20) yfit<-b[1] + poly(xfit,2,raw=T) %*% b[-1] plot(x,y) lines(xfit,yfit) But

[R] Probably dumb question about lists

2015-01-14 Thread Raphael Päbst
Hello everybody! I feel very stupid right now but suspect it has something to do with tiredness. I am trying to drop the last couple of Elements from a list and this doesn't work as expected. My code looks something like this: del <- deldir(x, y) tl <- tile.list(del) plot(tl) Now, I only want to

[R] using poly() to predict

2015-01-14 Thread Stanislav Aggerwal
This method of finding yhat as x %*% b works when I use raw polynomials: x<-1:8 y<- 1+ 1*x + .5*x^2 fit<-lm(y~poly(x,2,raw=T)) b<-coef(fit) xfit<-seq(min(x),max(x),length=20) yfit<-b[1] + poly(xfit,2,raw=T) %*% b[-1] plot(x,y) lines(xfit,yfit) But it doesn't work when I use orthogonal polynomials

Re: [R] R 3.1.2 mle2() function on Windows 7 Error and multiple optimal solutions

2015-01-14 Thread Ben Bolker
Mosab Alqurashi KSU.EDU.SA> writes: > > ?Hello, > I want to get the Maximum Likelihood Estimate (MLE) for an Extended > weibull with three parameters. I have two problems. First, some > times if I change the starting point for one of the parameters I got > an error message. Second, sometime wh

Re: [R] Avoid R shell process echoes of piped inputs from a different programming language on Linux

2015-01-14 Thread Duncan Murdoch
I don't know the answer to your question. If nobody else answers it, I'd suggest asking in R-devel instead, and include some simple code that illustrates the problem, e.g. using "echo" to pipe something into an R process and explaining what part of the output you want to suppress. Duncan Murdoch

Re: [R] two-sample KS test: data becomes significantly different after normalization

2015-01-14 Thread Andrews, Chris
Your definition of p-value is not correct. See, for example, http://en.wikipedia.org/wiki/P-value#Misunderstandings -Original Message- From: Monnand [mailto:monn...@gmail.com] Sent: Wednesday, January 14, 2015 2:17 AM To: Andrews, Chris Cc: r-help@r-project.org Subject: Re: [R] two-sam

[R] Avoid R shell process echoes of piped inputs from a different programming language on Linux

2015-01-14 Thread marcel Austenfeld
Dear R developers,   i have rather a developer question.   >From an external application (different programming language)  i contact R >through pipes (process pipes -> exec).   On Linux i use the R (bash) script to start the R process which will be available in a custom shell.   If i pipe a

Re: [R] par function for barchart

2015-01-14 Thread PIKAL Petr
Hi par is for base graphics, barchart is grid graphics, you can not easily mix those two together. You can myd<-melt(mydata) Using REGIONE as id variables barchart(REGIONE ~ value|variable ,data=myd) or with ggplot p<-ggplot(myd, aes(x=REGIONE,y=value)) p+geom_bar( stat="identity")+facet_grid(

Re: [R] two-sample KS test: data becomes significantly different after normalization

2015-01-14 Thread Martin Maechler
> Monnand > on Wed, 14 Jan 2015 07:17:02 + writes: > I know this must be a wrong method, but I cannot help to ask: Can I only > use the p-value from KS test, saying if p-value is greater than \beta, then > two samples are from the same distribution. If the definition

Re: [R] list of vectors which are part of an initial vector

2015-01-14 Thread Stefano Sofia
Thank you Jeff for your tips, I found what I was looking for. Stefano Da: Jeff Newmiller [jdnew...@dcn.davis.ca.us] Inviato: lunedì 12 gennaio 2015 9.13 A: Stefano Sofia; r-help@r-project.org Oggetto: Re: [R] list of vectors which are part of an initial ve

[R] par function for barchart

2015-01-14 Thread Massimiliano Tripoli
Dear R-users, I can't put 1 x 2 pictures on one plot by using barchart function of lattice package. I often use par function but I'm not able to do it. This is my code and data: mydata <- structure(list(REGIONE = c("ABRUZZO", "BOLZANO - BOZEN", "CALABRIA", "CAMPANIA", "EMILIA-ROMAGNA", "FRIULI

[R] R 3.1.2 mle2() function on Windows 7 Error and multiple optimal solutions

2015-01-14 Thread Mosab Alqurashi
?Hello, I want to get the Maximum Likelihood Estimate (MLE) for an Extended weibull with three parameters. I have two problems. First, some times if I change the starting point for one of the parameters I got an error message. Second, sometime when I change the starting point for one of the p

Re: [R] Links to vignettes in the base packages

2015-01-14 Thread Martin Maechler
> Richard Cotton > on Wed, 14 Jan 2015 09:23:42 +0300 writes: > Contributed packages have copies of their vignettes on CRAN (in their > package page at /web/packages/). > Since base packages no longer have a page here, I can't find a web link to them. > I'm aware th

[R] Links to vignettes in the base packages

2015-01-14 Thread Richard Cotton
Contributed packages have copies of their vignettes on CRAN (in their package page at /web/packages/). Since base packages no longer have a page here, I can't find a web link to them. I'm aware that I can find the vignette via browseVignettes() or vignette("vignettename", package = "packagename")

Re: [R] Help on Principal Component Analysis in R

2015-01-14 Thread Brandstätter Christian
Hi, it could also be that during reading your csv the matrix became a character matrix. Try str(x); this should show you if your values are seen as numeric. Also, if you have descriptive columns (e.g. date, labels, etc.) they must be removed before calculating PCR. With best regards Brandstätter