Re: [R] changing colors in filled.contour

2015-10-07 Thread William Dunlap
After mylevels <- seq(0,10,10) mylevels has length 2, hence you get 2 colors. Use length=10, Bill Dunlap TIBCO Software wdunlap tibco.com On Wed, Oct 7, 2015 at 3:40 AM, begrinner wrote: > Dear R community! > > I haven't worked with R before but I found it has

[R] how to suppress to bring value labels with memisc pkg

2015-10-07 Thread Yongnam Kim
Hello all, I know how to suppress to bring value labels from SPSS when importing sav file using foreign package by "use.value.labels = FALSE" in read.spss function. But, I don't know how to do the same thing when using memisc pkg, in particular, with the spss.system.file function. Many thanks!

Re: [R] Quantile Regression without intercept

2015-10-07 Thread Roger Koenker
> On Oct 7, 2015, at 3:46 PM, Preetam Pal wrote: > > So, what does weighted quantile regression even aim to achieve? Invariably, > this plane would not split the data set into the requisite fractions….. This officially ties this thread into a loop, since the OP wants

Re: [R] Error with Segmented package

2015-10-07 Thread David Winsemius
On Oct 7, 2015, at 6:50 AM, andrew haywood wrote: > Dear List, > > I am trying to run a simple pieewise regression using the segmented package. > > When running the following code > > library(segmented) > data = data.frame(x=c(50,60,70,80,90,100,110) , y= >

Re: [R] Quantile Regression without intercept

2015-10-07 Thread Preetam Pal
So, what does weighted quantile regression even aim to achieve? Invariably, this plane would not split the data set into the requisite fractions. -Original Message- From: "Roger Koenker" Sent: ‎06-‎10-‎2015 07:09 PM To: "Lorenz, David" Cc:

[R] RSM in R, optimize/minimize a response

2015-10-07 Thread simon.heintz
Good morning I'm trying to optimize (minimize actually) a response from a DOE with 4 factors. The 4 factors were built from a Latin Hypercube DOE design type. Then I proceeded this experiment on 28 different cases, so each case will include 2000 experiments. I would like to find the factor quartet

Re: [R] Measure the frequencies of pairs in a matrix

2015-10-07 Thread Hermann Norpois
Ok, this was misleading. And was not that important. My result matrix should look like this: 12 3 4 5 6 7 ... 1 p1 p2 2 p 3 4 p1 etc are the frequencies of the combinations 1 and 2 for instance do not appear in my example. So the values would be zero. Actually, this part is not

Re: [R] Installing pre-compiled R in Linux

2015-10-07 Thread Sasikumar Kandhasamy
Thanks a lot Mike. The Linux distribution we use is "Red Hat Enterprise Linux Server release 6.2". Also, couple of clarifications, 1. Do we have a R package compatibility matrix against the Linux kernel version? Or for the Red Hat Linux with kernel version 2.6.32-279, do you have any

[R] non-linear fit and statistical significance

2015-10-07 Thread Juan Munoz-Garcia
Dear R users, I’m using the constrOptim.nl to carry out non-linear fit with constraints between experimental and modelled data. As far as I understand the parameter “fval" indicates the goodness of this fit. My first question is: Is the fval parameter similar to commonly

Re: [R] Installing pre-compiled R in Linux

2015-10-07 Thread Michael Hannon
Hi, Sasi. Yes, I think that getting familiar with R in a "friendly" context before dealing with your embedded system is a good idea. Note that the non-GUI part of R is more or less identical across platforms. Hence, you could experiment with it on a Macbook, a Surface tablet, etc., etc., if

[R] ltraj contains irregular data

2015-10-07 Thread Wiebke Ullmann
Dear everyone, as far as I understood the setNA function in adhabitat is there to produce rows in a location dataset that have not been recorded by the GPS collar. I just can not figure out where the problem is. Hope somebody can help me. >daten$timestamp <-

[R] changing colors in filled.contour

2015-10-07 Thread begrinner
Dear R community! I haven't worked with R before but I found it has some nice abilities to create graphics. I made it to create a filled.contour with the settings I want but unfortunately I don't seem to be able to change colors. Instead of the standard colors, I'd like to get dark gray or black

[R] problem with dataframes

2015-10-07 Thread Pau Marc Muñoz Torres
Hello everybody, I am using the library rapidr of bioconductor, I am running it without problems as follows; library("RAPIDR") makeBinnedCountsFile("SRR611842_2.bam","SRR611842_2.bam.bai","binnecounts.r",mask=NULL,k=2) > SampleID=c("SRR611842_2") > DX=c("T21") > Gender=c("Male") > frame=

[R] additive proportional odds model in R

2015-10-07 Thread Neverstop
Hi all! I'm looking for a way to fit an additive proportional odds model in R since I need to analyse ordinal data. The polr() function of the MASS package of R allows only to fit a linear proportional odds model so I can't use it. The gam() function of the gam package of R allows to fit

[R] MSwM R package: model and state probability significance (p00+p11=1)

2015-10-07 Thread Maxim Fomin
Dear R community, It seems that MSwM package provides three model diagnostic tools: pooled residuals plot, QQ plot and ACFs plot (via plotDiag). However, I cannot find any statistic to test the whole model significance (and even for each reported regime only multiple R^2 is provided). After

[R] Error with Segmented package

2015-10-07 Thread andrew haywood
Dear List, I am trying to run a simple pieewise regression using the segmented package. When running the following code library(segmented) data = data.frame(x=c(50,60,70,80,90,100,110) , y= c(703.786,705.857,708.153,711.056,709.257, 707.4, 705.6)) model.lm = segmented(lm(y~x,data = data),seg.Z

[R] rpart cutpoint interpretation

2015-10-07 Thread Kishor Tappita
Hi All, I am trying to derive cutpoint/threshold with a poisson distributed dependent variable. I know how to interpret cutpoint with binary dependent variable based on direction. Can some on help me to intrepret cutpoint for poisson case with one independent variable with the derived threshold.

Re: [R] extract all dataframes from list

2015-10-07 Thread PIKAL Petr
Hi Or basically you can use for cycle. This can be as effective as lapply and if you plan to do some operations on your data frames maybe more manageable. Cheers Petr > -Original Message- > From: R-help [mailto:r-help-boun...@r-project.org] On Behalf Of maicel > Sent: Tuesday, October

Re: [R] Bug in auglag?

2015-10-07 Thread Rainer M Krug
Thanks a lot Ravi for the clarification and it makes sense - bug in my understanding acknowledged. I'll change the objective function as suggested. Thanks, Rainer Ravi Varadhan writes: > Dear Rainer, > This is NOT a bug in auglag. I already mentioned that auglag()

Re: [R] extract all dataframes from list

2015-10-07 Thread Gerrit Eichner
Hello, Maicel, if you only want to extract info from those data frames, maybe it is enough to attach the list to the search path using attach (and afterwards detach()); see ?attach. Otherwise the online help page of assign() might be a starting point for you. Hth -- Gerrit Hi Or

Re: [R] Installing different R versions

2015-10-07 Thread Loris Bennett
Dear Luca, Luca Cerone writes: > Dear all, > on one shared machine we have an older R version installed. Some packages > have known issues with that version that are fixed in newer R versions. > > Since that is a production machine with many jobs running we would like to

Re: [R] Polygon shade

2015-10-07 Thread Jim Lemon
Hi bgnum, You can try something like this: testdates<-as.Date(paste(1:30,"Sep",2015),"%d %b %Y") nemails<-sample(10:30,30,TRUE) library(plotrix) stackpoly(testdates,nemails,col="blue") Jim On Wed, Oct 7, 2015 at 12:20 AM, bgnumis bgnum wrote: > Hi All, > > I want to shade

[R] Installing different R versions

2015-10-07 Thread Luca Cerone
Dear all, on one shared machine we have an older R version installed. Some packages have known issues with that version that are fixed in newer R versions. Since that is a production machine with many jobs running we would like to keep things as they are. However I would also like to keep

Re: [R] vector graphics

2015-10-07 Thread Ivan Calandra
Thanks Jean for the tip. I'll try postscript() and devEMF::emf() and see if it works. Still, it's very complicated to export vector graphics in R... Ivan -- Ivan Calandra, PhD University of Reims Champagne-Ardenne GEGENAA - EA 3795 CREA - 2 esplanade Roland Garros 51100 Reims, France +33(0)3 26

Re: [R] Installing pre-compiled R in Linux

2015-10-07 Thread Michael Hannon
I don't think kernel compatibility is a significant issue for most applications. I can say for certain that I update the kernels on my linux boxes without having to reinstall R. There *are* R packages for RHEL and friends. Have a look at: https://cran.r-project.org/bin/linux/redhat/README

Re: [R] Installing pre-compiled R in Linux

2015-10-07 Thread Sasikumar Kandhasamy
Thanks a lot Mike & others for the valuable input. I was wrong on debian linux on the embedded box, it is monta vista embedded linux. I hope, there is no R core package available for monta vista distribution and i may need to cross compile from R source. Please correct me if i am wrong. Thanks &

Re: [R] Installing pre-compiled R in Linux

2015-10-07 Thread Sasikumar Kandhasamy
Thanks Mike. From source, i am able to compile and use R in my red hat linux box. I was trying to get hands on using R in linux box before trying R on the embedded box. My requirement is to run R (R scripts) on my embedded box which has customized debian linux (kernel version 2.6.32) in batch

Re: [R] Error with Segmented package

2015-10-07 Thread andrew haywood
Thanks David, I apologies for not posting the versions. I am running R version 3.2.2 (2015-08-14) Platform: x86_64-w64-mingw32/x64 (64-bit) Running under: Windows 7 x64 (build 7601) Service Pack 1 locale: [1] LC_COLLATE=English_United Kingdom.1252 [2] LC_CTYPE=English_United Kingdom.1252 [3]

Re: [R] additive proportional odds model in R

2015-10-07 Thread David Winsemius
On Oct 7, 2015, at 3:33 AM, Neverstop wrote: > Hi all! > I'm looking for a way to fit an additive proportional odds model in R since > I need to analyse ordinal data. > The polr() function of the MASS package of R allows only to fit a linear > proportional odds model so I can't use it. Can

Re: [R] Error with Segmented package

2015-10-07 Thread David Winsemius
On Oct 7, 2015, at 6:10 PM, andrew haywood wrote: > Thanks David, > > I apologies for not posting the versions. > > I am running > > R version 3.2.2 (2015-08-14) > Platform: x86_64-w64-mingw32/x64 (64-bit) > Running under: Windows 7 x64 (build 7601) Service Pack 1 > > locale: > [1]

[R] duplicate rows with rbind in a loop

2015-10-07 Thread Cara Fiore
Dear R users, I wrote a simple script to change the header lines in a fasta file that contains DNA sequences in a format: >header1 sequence1 >header2 sequence2 I am basically trying to replace the "header" in this file with a line from another file (taxonomy file). In order to do that I have to

Re: [R] [lattice::xyplot] Using (panel:.)abline with panel.superpose?

2015-10-07 Thread Duncan Mackay
Forgot to send to list -Original Message- From: Duncan Mackay [mailto:dulca...@bigpond.com] Sent: Thursday, 8 October 2015 08:44 To: 'Szumiloski, John' Subject: RE: [R] [lattice::xyplot] Using (panel:.)abline with panel.superpose? Hi John I only got grid lines on your # code 2 with

[R] Time Series Daily Frequency Part of a Year

2015-10-07 Thread zod jones
My data consist of daily sales figures for multiple products but only for a 3 month period each year (Oct., Nov., Dec). The goal is to forecast the daily sales figures for the following year *by day* (the following Oct., Nov., Dec.) So, I'd like to forecast what sales for product X will be on Nov.

Re: [R] Installing pre-compiled R in Linux

2015-10-07 Thread peter dalgaard
On 07 Oct 2015, at 13:05 , Jeroen Ooms wrote: >> On Tue, Oct 6, 2015 at 9:42 PM, Sasikumar Kandhasamy >> wrote: >>> Thanks a lot Mike. The Linux distribution we use is "Red Hat Enterprise >>> Linux Server release 6.2". > > On RHEL and CentOS the

Re: [R] Installing pre-compiled R in Linux

2015-10-07 Thread Jeroen Ooms
> On Tue, Oct 6, 2015 at 9:42 PM, Sasikumar Kandhasamy > wrote: >> Thanks a lot Mike. The Linux distribution we use is "Red Hat Enterprise >> Linux Server release 6.2". On RHEL and CentOS the easiest and most reliable way to get R and R packages is via EPEL. Simply add the

Re: [R] Measure the frequencies of pairs in a matrix

2015-10-07 Thread Boris Steipe
Still not sure I understand. But here is what I think you might mean: # Your data mat <- structure(c(5, 6, 5, 5, 4, 3, 6, 7, 4, 7, 5, 5, 5, 5, 6, 5, 5, 4, 3, 6, 7, 4, 7, 5, 5, 5, 6, 5, 4, 5, 5, 7, 5, 6, 3, 5, 6, 7, 6, 6, 5, 4, 5, 5, 7, 5, 6, 3, 5, 6, 7, 6), .Dim = c(26L, 2L)) # Create a square

Re: [R] vector graphics

2015-10-07 Thread Ivan Calandra
For the record, it seems that devEMF::emf() works fine on the Windows Vista SP2 machine with R3.2.2 and Illustrator CS4: the text is recognized as text and every point/line can be dissociated. Ivan -- Ivan Calandra, PhD University of Reims Champagne-Ardenne GEGENAA - EA 3795 CREA - 2

[R] F-test of equality of variances - but weighted?

2015-10-07 Thread Dimitri Liakhovitski
I would like to use an F-Test for Equality of Variances on a variable to compare two groups. Normally, this would be done with 'var.test'. However, the data need to be weighted (individual-level weights). R's package 'survey' is geared at running analyses with complex sampling weights. But,

Re: [R] Installing different R versions

2015-10-07 Thread Wolfgang Raffelsberger
check out the official document *R Installation and Administration* from https://cran.r-project.org/manuals.html There you'll find how to define a specific path for each installation. (Since a number of years I administrate multiple versions of R at different platforms, of course including Linux)

Re: [R] Installing pre-compiled R in Linux

2015-10-07 Thread stephen sefick
I couldn't tell from the OP's message what distribution they have installed. If it is redhat (or a derivative), the extra packages for enterprise linux (epel) has up-to-date R packages to install with yum. I have had some minor issues with installing A FEW (mostly GIS related) packages that were

Re: [R] Measure the frequencies of pairs in a matrix

2015-10-07 Thread David L Carlson
As with Boris, I'm not sure what you are looking for, but this may help > # To get all possibilities, create a grid > grd <- expand.grid(0:9, 0:9) > # Extract those with smaller first column values > grd <- grd[grd$Var1 <= grd$Var2,] > # Tabulate after pasting first and second column > grd2 <-

[R] [lattice::xyplot] Using (panel:.)abline with panel.superpose?

2015-10-07 Thread Szumiloski, John
Dear useRs, I recently had a query concerning how to customize the graphics parameters in lattice::xyplot to change not only the color but also the pch symbols, lty and lwd line parameters, etc., within each grouping variable in the plot.

Re: [R] Measure the frequencies of pairs in a matrix

2015-10-07 Thread William Dunlap
You could also call table() on the columns of the input matrix, first converting them to factors with levels 1:max. Then add together the upper and lower triangles of the table if order is not important. E.g., f2 <- function (mat) { maxMat <- max(mat) stopifnot(is.matrix(mat), all(mat