Re: [R] [FORGED] find numbers that fall in a region or the next available.

2016-02-03 Thread Alaios via R-help
That is a great tip thanks.That would indeed bring me points that are the closes to my area.. but if I am not wront that returns points that are part of a circle surface. It might be that I get a point that is just 50 meters outside of my map area. Is not that true? I would need after I find clo

Re: [R] [FORGED] Generate arrival of time based on uniform random number

2016-02-03 Thread Rolf Turner
This looks very much like a homework problem and this list has a "no homework" policy. cheers, Rolf Turner -- Technical Editor ANZJS Department of Statistics University of Auckland Phone: +64-9-373-7599 ext. 88276 On 04/02/16 17:08, smart hendsome via R-help wrote: Hi everyone, I have prob

[R] Generate arrival of time based on uniform random number

2016-02-03 Thread smart hendsome via R-help
Hi everyone, I have problem regarding to generate arrival of time based on uniform random number. Let day0 = 0.383, lambda = 0.2612 1) Generate uniform random number (already settled) using the code below: set.seed(1234) rand.no <- function(n,itr){   matrix(runif(n*itr, 0, 1), nrow=n, ncol=itr) }

Re: [R] LDheatmap

2016-02-03 Thread Val
Thank you Bert, Yes I looked a this one and I was looking for if any one has used it or not before? My data set is different what they are showing in the paper On Wed, Feb 3, 2016 at 4:00 PM, Bert Gunter wrote: > Have you looked here (found immediately by an internet search!)? > > https://

Re: [R] Create macro_var in R

2016-02-03 Thread Gabor Grothendieck
See Example 5. Insert Variables on the sqldf home page. https://github.com/ggrothendieck/sqldf On Wed, Feb 3, 2016 at 2:16 PM, Amoy Yang via R-help wrote: > First, MVAR<-c("population) should be the same as "population'". Correct? > You use tab[[MVAR]] to refer to "population" where doub

Re: [R] [FORGED] Re: determine the year of a date

2016-02-03 Thread Jim Lemon
Hi all, If you don't want Excel to surreptitiously "correct" dates that are not in mm/dd/ format, always specify international format -mm-dd in Excel. Jim On Thu, Feb 4, 2016 at 8:26 AM, peter dalgaard wrote: > Fortune candidate... > > -pd > > > On 03 Feb 2016, at 22:13 , Rolf Turner

Re: [R] LDheatmap

2016-02-03 Thread Bert Gunter
Have you looked here (found immediately by an internet search!)? https://cran.r-project.org/web/packages/LDheatmap/vignettes/LDheatmap.pdf Cheers, Bert Bert Gunter "The trouble with having an open mind is that people keep coming along and sticking things into it." -- Opus (aka Berkeley Breath

[R] LDheatmap

2016-02-03 Thread Ashta
Hi all, I am looking for an R package that calculates a pair wise LD (linkage disequilibrium) I came up with library(LDheatmap). has any one used this library? I would appreciate if I get a help how to use this library for my set of data.. My data set look like Geno file Name1 1 1 2 2 2 2 N

Re: [R] [FORGED] Re: determine the year of a date

2016-02-03 Thread peter dalgaard
Fortune candidate... -pd > On 03 Feb 2016, at 22:13 , Rolf Turner wrote: > > On 04/02/16 09:55, Boris Steipe wrote: >> Who said this was Excel? What did I miss? > > If data have been fucked up, it is odds-on that Excel is to blame. > > cheers, > > Rolf Turner > > -- > Technical Editor ANZJ

Re: [R] [FORGED] Re: determine the year of a date

2016-02-03 Thread Rolf Turner
On 04/02/16 09:55, Boris Steipe wrote: Who said this was Excel? What did I miss? If data have been fucked up, it is odds-on that Excel is to blame. cheers, Rolf Turner -- Technical Editor ANZJS Department of Statistics University of Auckland Phone: +64-9-373-7599 ext. 88276

Re: [R] determine the year of a date

2016-02-03 Thread Boris Steipe
Who said this was Excel? What did I miss? B. On Feb 3, 2016, at 3:52 PM, Ista Zahn wrote: > On Wed, Feb 3, 2016 at 3:43 PM, Jeff Newmiller > wrote: >> Not exactly. >> >> If you have been "Excel"ed and have both formats in one column then you >> really have a data cleanup task to do. > > No

Re: [R] determine the year of a date

2016-02-03 Thread Ista Zahn
On Wed, Feb 3, 2016 at 3:43 PM, Jeff Newmiller wrote: > Not exactly. > > If you have been "Excel"ed and have both formats in one column then you > really have a data cleanup task to do. No need to reinvent the wheel. Just do library(lubridate) mdy(c("09/01/15", "09/01/2015")) Best, ista The be

Re: [R] determine the year of a date

2016-02-03 Thread Jeff Newmiller
Not exactly. If you have been "Excel"ed and have both formats in one column then you really have a data cleanup task to do. The best route to resolution is to fix the source of the data (e.g. go back into Excel and reformat the column with consistent formatting). If you need to do this on an

Re: [R] Create macro_var in R

2016-02-03 Thread David Winsemius
> On Feb 3, 2016, at 11:57 AM, Amoy Yang via R-help > wrote: > > Right! the following works to r but not sqldf. > MVAR <- "population" > tab[[ MVAR ]] > sqldf("select tab[[MVAR]] from tab") If ypu need a string to pass to sqldf, then use the `paste` function to build a single string: sqldf(

Re: [R] determine the year of a date

2016-02-03 Thread Boris Steipe
As the documentation says, that's exactly what the expansion will do by default. B. On Feb 3, 2016, at 2:31 PM, carol white wrote: > yes, some of them are like 09/01/15 and some others are 09/01/2015 and all of > them range between 2015 and 2016. The goal was to convert 09/01/15 to > 09/01/2

Re: [R] Create macro_var in R

2016-02-03 Thread Amoy Yang via R-help
Right! the following works to r but not sqldf. MVAR <- "population" tab[[ MVAR ]] sqldf("select tab[[MVAR]] from tab") On Wednesday, February 3, 2016 1:18 PM, Amoy Yang via R-help wrote: First, MVAR<-c("population) should be the same as "population'". Correct? You use tab[[MVAR]] to re

Re: [R] Create macro_var in R

2016-02-03 Thread MacQueen, Don
Repeat many times: R does not work the same way as SAS, do not expect it to If I needed to construct an SQL statement in which the name of a field is provided by the value of another variable, I would consider this: key <- 'pop' sql.stmt <- paste("select grade, count(*) as cnt, min(",key,") a

Re: [R] Install.package CAIC

2016-02-03 Thread JLucke
You should consider the package "caper", which provides a CAIC. From the caper manual "The caper package implements the methods originally provided in the programs CAIC (Purvis and Rambaut, 1995b) and MacroCAIC (Agapow and Isaac, 2002)." Hope this helps. Gabriela Wofkova Sent by: "R-help"

Re: [R] [FORGED] find numbers that fall in a region or the next available.

2016-02-03 Thread David L Carlson
Look at the point.in.polygon() and over() functions in package sp. - David L Carlson Department of Anthropology Texas A&M University College Station, TX 77840-4352 -Original Message- From: R-help [mailto:r-help-boun...@r-project.org] On Behalf Of Alaios

Re: [R] determine the year of a date

2016-02-03 Thread carol white via R-help
yes, some of them are like 09/01/15 and some others are 09/01/2015 and all of them range between 2015 and 2016. The goal was to convert 09/01/15 to 09/01/2015. I don't know if the fact that they range between 2015 and 2016 make the task easier. Best wishes Carol On Wednesday, February 3, 20

Re: [R] Create macro_var in R

2016-02-03 Thread Amoy Yang via R-help
First, MVAR<-c("population) should be the same as "population'". Correct? You use tab[[MVAR]] to refer to "population" where double [[...]] removes double quotes "...", which seemingly work for r-code although it is tedious in comparison direct application in SAS %let MVAR=population. But it does

Re: [R] Install.package CAIC

2016-02-03 Thread Boris Steipe
I don't see a "CAIC" package on CRAN or bioconductor. Are you sure you got that right? Perhaps you meant another package? Or a function in another package? E.g. ape has a read.caic() function ... See https://cran.r-project.org/web/views/Phylogenetics.html B. (Please also read the posting guide

Re: [R] determine the year of a date

2016-02-03 Thread Boris Steipe
Sorry - messed up example: corrected here... d <- "7/27/77" strptime(d, format="%m/%d/%y") # "1977-07-27 EDT" x <- strptime(d, format="%m/%d/%y") strftime(x, format="%m/%d/%Y") # "07/27/1977" On Feb 3, 2016, at 1:29 PM, Boris Steipe wrote: > It seems your autocorrect is playing tricks on y

Re: [R] Create macro_var in R

2016-02-03 Thread Jeff Newmiller
Are you perhaps needing to (re-)read the discussion on indexing in the "Introduction to R" document that comes with the software? (That is a common deficiency...) It looks to me like you want something like MVAR <- "population" tab[[ MVAR ]] -- Sent from my phone. Please excuse my brevity. O

Re: [R] Create macro_var in R

2016-02-03 Thread Bert Gunter
I am not sure what you want, but: 1) If you have not already done so, go through an R tutorial or two. For some suggestions: https://www.rstudio.com/resources/training/online-learning/#R R is quite different than SAS. 2) If I misunderstand, perhaps ?within is what you are looking for. Cheer

Re: [R] Create macro_var in R

2016-02-03 Thread ruipbarradas
Hello, You can't use tab$MVAR but you can use tab[[MVAR]] if you do MVAR <- "population" (no need for c()). Hope this helps, Rui Barradas   Citando Amoy Yang via R-help : > population is the field-name in data-file (say, tab). > MVAR<-population takes data (in the column of population) rat

Re: [R] Create macro_var in R

2016-02-03 Thread Duncan Murdoch
On 03/02/2016 1:23 PM, Amoy Yang wrote: population is the field-name in data-file (say, tab). MVAR<-population takes data (in the column of population) rather than field-name as done in SAS: %let MVAR=population; In the following r-program, for instance, I cannot use ... tab$MVAR...or simply M

Re: [R] determine the year of a date

2016-02-03 Thread Boris Steipe
It seems your autocorrect is playing tricks on you but if I understand you correctly you have a two digit year and want to convert that to a four digit year? That's not uniquely possible of course; by convention, as the documentation to strptime() says: On input, values 00 to 68 are prefixed

Re: [R] OLS Regression on subset of data

2016-02-03 Thread Bert Gunter
Did you not read about the "subset" argument of lm() ? Cheers, Bert Bert Gunter "The trouble with having an open mind is that people keep coming along and sticking things into it." -- Opus (aka Berkeley Breathed in his "Bloom County" comic strip ) On Wed, Feb 3, 2016 at 10:09 AM, Preetam Pal

Re: [R] Problems loading lme4

2016-02-03 Thread Jeff Newmiller
R 3.3.0 is not a released version of R. It appears that you have obtained a version of lme4 that is ahead of its time. This would have been a perfect time for you to have followed the Posting Guide mentioned below (in particular the inclusion of complete version information such as the sessionIn

Re: [R] Create macro_var in R

2016-02-03 Thread Amoy Yang via R-help
population is the field-name in data-file (say, tab). MVAR<-population takes data (in the column of population) rather than field-name as done in SAS:  %let MVAR=population; In the following r-program, for instance, I cannot use ... tab$MVAR...or simply MVAR itself since MVAR is defined as "popu

[R] OLS Regression on subset of data

2016-02-03 Thread Preetam Pal
Hi, I am performing OLS regression on a dataset involving variables Y, X1 and X2, each having 500 observations. I want to perform the regression only on a subset of the dataset (say, using observations 50 till 350). Is there any way in which I can pass the entire dataset, but somewhere select thes

Re: [R] Create macro_var in R

2016-02-03 Thread Duncan Murdoch
On 03/02/2016 12:41 PM, Amoy Yang via R-help wrote: There is a %LET statement in SAS: %let MVAR=population; Thus, MVAR can be used through entire program. In R, I tried MAVR<-c("population"). The problem is that MAVR comes with double quote "" that I don't need. But MVAR<-c(population) di

Re: [R] Create macro_var in R

2016-02-03 Thread Henrik Bengtsson
Don't know what `population` is, but a simple assignment MVAR <- population may provide what you need. Note, no c(). For example, > foo <- rnorm > foo(3) [1] -0.08093862 -0.87827617 1.52826914 /Henrik On Wed, Feb 3, 2016 at 9:41 AM, Amoy Yang via R-help wrote: > There is a %LET statement

[R] Create macro_var in R

2016-02-03 Thread Amoy Yang via R-help
There is a %LET statement in SAS: %let MVAR=population; Thus, MVAR can be used through entire program. In R, I tried MAVR<-c("population"). The problem is that MAVR comes with double quote "" that I don't need. But MVAR<-c(population) did NOT work out. Any way that double quote can be remov

[R] Install.package CAIC

2016-02-03 Thread Gabriela Wofkova
​H ​ello, I am a begginer in R and I need to incorporate phylogeny into my data. So i tried to instal package CAIC. > install.packages("CAIC") > install.packages("CAIC", repos="http://R-Forge.R-project.org",type="source";) > install.packages("CAIC", repos="http://R-Forge.R-project.org";) But it

[R] Problems loading lme4

2016-02-03 Thread Sarah Brennenstuhl
Hi there! I having a problem loading lme4. When I try to load the package I get the following error message: Error : object 'sigma' is not exported by 'namespace:stats' In addition: Warning message: package 'lme4' was built under R version 3.3.0 Error: package or namespace load failed for 'lme4'

[R] Keep factor order in lmer

2016-02-03 Thread Pascal A. Niklaus
I would like to keep a specific order of fixed effects in a model passed to lmer. In particular, I would like to prevent that interactions are automatically moved after all main effects. In aov and lme, this is possible with terms(..., keep.order=TRUE). Unfortunately, I have not found a way to

Re: [R] find numbers that fall in a region or the next available.

2016-02-03 Thread S Ellison
> Dear all,I have GPS coordinates (one vector for longitude and one for > latitude: > GPSLong and GPSLat) of small are that is around 300meters X 300 meters > (location falls inside UK).At the same time I have two more vectors (Longitude > and Latitude) that include position of food stores again t

Re: [R] [R-pkgs] New Package: stripless (V. 1.0)

2016-02-03 Thread Michael Friendly
Hi Burt This looks like an interesting package. However, you should know that 'strucplot()' might not be the best choice for your generic, because vcd::strucplot() is now well-established as the general name for functions plotting structured multi-way frequency tables, like mosaic plots and

[R] determine the year of a date

2016-02-03 Thread carol white via R-help
Hi,might be trivial but how to determine the year of a date which is in the %m/%d/%y format and those whose year is century should be modified to ISO so that all date will have with year in ISO? Regards, Carol [[alternative HTML version deleted]]

[R] vertical lines in riverplot

2016-02-03 Thread Brandstätter Christian
Dear List, I want to draw a Sankey-diagramm in R, for which I would use the riverplot package. It would be really nice to be able to 1. draw vertical lines (="edges") between nodes on the same x-axis (T03-T04 in the example below) and 2. to invert the plot. Now the "edge" between T03 and T04 is st

Re: [R] [FORGED] find numbers that fall in a region or the next available.

2016-02-03 Thread Alaios via R-help
Thanks. I am using distm of the geoshere package.I still wonder if there is a package that can tell me if a gps coordinate or not falls inside my area that is defined as: bbox <- c(min(PlotPoints[, 1])-0.001, min(PlotPoints[, 2])-0.001, max(PlotPoints[, 1])+0.001, max(PlotPoints[, 2])+0.001) Pl

Re: [R] Use predict() after lmer() (library lme4)

2016-02-03 Thread Thierry Onkelinx
Dear Marc, This question is more suited for R-Sig-mixed models to which I'm forwarding it. Your manual predictions for out3 are wrong. Here are the correct manual predictions. fixed <- model.matrix(~effect, data = dataF) %*% fixef(out3) random <- rowSums(model.matrix(~Sector, data = dataF) * ran