Re: [R] Cross correlation between two time series over nested time periods?

On 2015-05-14 , at 02:11, Tim via R-help wrote: Hello Tim, Re: > I have two time series > > > Calculate and plot cross correlation between two time series over nested time > periods. Each point in either time series is for a week (not exactly a > calendar week, but the first week in a c

[R] Dropping predictor variables based on adjusted R square

Hello experts, I have recently (1month) started using R. Earlier I was using SAS to work on analytic assignments. In SAS there is an option - forward selection, backward selection, step wise selection where in it removes the least impacting predictor variable from the set of variables based on a

Re: [R] Dropping predictor variables based on adjusted R square

Hi Shiv82, For a start, look at the "step" function in the stats package. There are a number of functions in other packages that offer different variable selection procedures. Jim On Thu, May 14, 2015 at 9:35 PM, Shivi82 wrote: > Hello experts, > > I have recently (1month) started using R. Earl

[R] Need help with lm function on MAC OS X. R version - 3.2.0

I Have a data frame named BSE and CP is my independent variable and here;s the error I get if I try to run an lm function any idea whats wrong P.S - all my data is in numeric except company which is a factor. I have 2700 row and 450 variable P.P.S - I have no missing data, I have 0 in Empty fiel

[R] Counting consecutive events in R

Hi, I have the following dataframe structure(list(Type = c("QRS", "QRS", "QRS", "QRS", "QRS", "QRS", "QRS", "QRS", "QRS", "QRS", "QRS", "QRS", "RR", "RR", "RR", "PP", "PP", "PP", "PP", "PP", "PP", "PP", "PP", "PP", "QTc", "QTc", "QTc", "QTc", "QTc", "QTc", "QTc", "QTc", "QTc", "QTc", "QTc", "QTc"

Re: [R] binding two lists of lists of dataframes together

Hi David, Understood. I've done a lot of work around this structure, but if there isn't anyway to change it back to the original structure then I'll work around it. Could you possibly show me how to add colnames to list3a? c("id","WgtBand","Wgt","Held","LID","Issuer","Bid","Offer") these wo

[R] create a vector from several data frames

Dear r-users, suppose that I have 20 data frames df1, df2, ..., df20 (one for each different location) with the same column names and column types (the first column contains a date, the others are numeric) like day tmax tmin 2015-05-10 20 10 2015-05-11 21 12 2015-05-12 17 9 2015-05-13 24 13 2015

Re: [R] Counting consecutive events in R

Assuming I understand the problem correctly, you want to check for runs of at least length five where both Score and Test_desc assume particular values. You don't care where they are or what other data are associated, you just want to know if at least one such run exists in your data frame. Here's

Re: [R] binding two lists of lists of dataframes together

On May 14, 2015, at 7:15 AM, Vin Cheng wrote: > Hi David, > > Understood. I've done a lot of work around this structure, but if there > isn't anyway to change it back to the original structure then I'll work > around it. I suppose there might be a way to make it back into a data.table, but

Re: [R] error using the huge R package

On May 13, 2015, at 3:45 PM, Juan Fernandez wrote: > Dear List > > I am trying to do an association network using some expression data I have, > the data is really huge: 300 samples and ~30,000 genes. I would like to apply > a gaussian graphical model to my data using the huge R package. > >

Re: [R] binding two lists of lists of dataframes together

Thanks David! Have a good one! > Subject: Re: [R] binding two lists of lists of dataframes together > From: dwinsem...@comcast.net > Date: Thu, 14 May 2015 07:58:43 -0700 > CC: r-help@r-project.org > To: newrnew...@hotmail.com > > > On May 14, 2015, at 7:15 AM, Vin Cheng wrote: > > > Hi Davi

Re: [R] create a vector from several data frames

If you combine all of the df's into a list, e.g. dfn <- paste0("df", 1:20) df <- lapply(dfn, get) names(df) <- dfn and if "target" is the day you want in the same date/time format as the day variable in the data frames: sapply(df, function(x) x[x$day==target, "tmax"]) will return a named vecto

Re: [R] Plotting times at night and getting plot limits correct

I might do it this way using Jim's sample data: epoch <- as.POSIXct( "1970-01-01" ) # any date you like dta <- data.frame( Timestamps = epoch + as.difftime( ifelse( Times >= 5/24 , Times

[R] Determinant

Hello, I am Vasilis Chasiotis and I am a candidate Ph.D. student in Aristotle University of Thessaloniki in Greece. I have the following problem. I want to check if the square of a number ( for example the square of 1.677722e+29 ) is an integer. The problem is that this number is the cal

[R] mixed-and weights

Hi,   I´m beginner with R and I hope that someone can help me.   I need to know how make command weights in mixed-effect model wit random intercept. Command weights=xxx is not function in nlme package. But I think that this command is commonly used for linear models.   Should I use a d

Re: [R] Counting consecutive events in R

I normally use rle() for these problems, see ?rle. for instance, k <- rbinom(999, 1, .5) series <- function(run) {

[R] specific package to laply

Hi, I have a batch jobs problem for parallel code without sudo. When I try to send a package to the different nodes, as follows: .libPaths( c(.libPaths(), "/my/first/library", "/my/second/library") ) library(foreach) library(iterators) library(parallel) librar

Re: [R] specific package to laply

I suggest you post using plain text to minimize communication problems on this list. I use the clusterEvalQ and cluster export functions to setup the slave processes before I start processing. --- Jeff Newmiller

Re: [R] Determinant

Of course the number _is_ an integer. It seems you are asking whether that integer can be exactly _represented_ on a computer? That depends on your processor (eg. 32/64 bit) and the size of the number; alternatively you could calculate with arbitrary precision with the Rmpfr package. For more i

[R] Error: Gradient function might be wrong - check it! for system of equations using Optimx

Hello,I am estimating a system of 3 equations using a nonlinear GMM.The first and second equations consists of a linear part and a nonlinear part. The nonlinear parts is specified as a logit and the logit as a whole has a coefficient attached to it. For the first two equations, it is specified a

Re: [R] Determinant

On May 14, 2015, at 1:44 AM, chasiot...@math.auth.gr wrote: > Hello, > > I am Vasilis Chasiotis and I am a candidate Ph.D. student in Aristotle > University of Thessaloniki in Greece. > > I have the following problem. > > I want to check if the square of a number ( for example the square of

Re: [R] Determinant

On May 14, 2015, at 12:50 PM, David Winsemius wrote: > > On May 14, 2015, at 1:44 AM, chasiot...@math.auth.gr wrote: > >> Hello, >> >> I am Vasilis Chasiotis and I am a candidate Ph.D. student in Aristotle >> University of Thessaloniki in Greece. >> >> I have the following problem. >> >> I

Re: [R] Cross correlation between two time series over nested time periods?

Thanks! Here the period of my time series B is a proper subinterval of the period of A. Does ccf(A,B) requires A and B span the same period? If A and B don't span the same period, what does ccf do? When moving B along the period of A by a lag, does ccf(A,B) calculate the cross correlation betwe

[R] Reading a tiff image

Dear R users, I am new to R and trying to learn raster image processing in R. I have a tiff image with pixel values and I am trying to read this image in R so that I can make some calculations and process the image. Could you please tell me how to read a tiff image in R? Thanks! -- Regards, *

[R] Reading access file

Hello everybody. I have a access file to read in R but I can’t to do this. I used Hmisc package, but it doesn’t work. Someone has the commands to read this kind of file? I attached the access file. Thanks. Silvano. __ R-help@r-project.org mailing

[R] Deparse substitute assign with list elements

Hello, When I use function `foo` with list elements (example 2), it defines a new object named `b[[1]]`, which is not what I want. How can I change the function code to show the desired behaviour for all data structures passed to the function? Or is there a more appropriate way to sort of "pa

Re: [R] Reading access file

Hi Silvano, your attachment did not go through (only a couple of formats are supported). Please try to give your code next time and what error you get but first read this, it might help: http://www.statmethods.net/input/dbinterface.html best, daniel Fel

Re: [R] Reading a tiff image

Hello, please set your mailing program to plain text instead of HTML if you post to this list. If you google read tiff in R you end up with the package called 'tiff'. (http://cran.r-project.org/web/packages/tiff/tiff.pdf) In case you do not know how to install and load a package, please consult

Re: [R] Reading a tiff image

There are many easy ways to answer this kind of question yourself. One possibility: go to www.rseek.org and search for "read tiff file". That will quickly suggest the tiff package, which you can install from CRAN, and also the rtiff package, ditto. It will also link to previous discussions of tha

Re: [R] Reading access file

mdbConnect<-odbcConnectAccess("C:\\temp\\demo.mdb"); sqlTables(mdbConnect); demo<-sqlFetch(mdbConnect, "tblDemo"); odbcClose(mdbConnect); rm(demo); On Thu, May 14, 2015 at 6:31 AM, silvano wrote: > Hello everybody. > > I have a access file to read in R but I can’t to do this. > > I used Hmisc pa

Re: [R] Deparse substitute assign with list elements

You could use a 'replacement function' named 'bar<-', whose last argument is called 'value', and use bar(variable) <- newValue where you currently use foo(variable, newValue). bar <- function(x) { x + 3 } `bar<-` <- function(x, value) { bar(value) } a <- NA bar(a) <- 4 a # [1] 7 b <- list

Re: [R] Determinant

Note that when using double precision arithmetic you cannot tell if that number is different from other numbers within 10 billion of it -- they all have the same representation. > 1.677721599849e+29 == 1.677721599849e+29 + 1e10 [1] TRUE You need to use other, time consuming, methods (e.g.

Re: [R] Reading access file

Be sure to use the 32bit version of R with the code below. This is a limitation of the ODBC driver for "mdb" files. If you use the odbcConnectAccess2007() function to connect to "accdb" files then you can use 64bit R. --- Jef

[R] Interpreting GLM Interaction Contrasts in R (using glht)

Hello, I am trying to do a BACI analysis on count data, and I am having trouble interpreting the output from multcomp::glht. I don't understand how the contrast's coefficients are related to effect size (if at all??). I have 5 treatment conditions (one is a control), and I have counts from

[R] Using iconv from glibc

The R Installation and Administration manual section A.1 states that glibc should provide a suitable iconv function, but I can't get R's configure script to accept/validate iconv on a platform I need to support using glibc 2.20. Is glibc is actually compatible (or if gnu libiconv is essentially

Re: [R] Using iconv from glibc

This belongs on R-devel. Read the Posting Guide, which also warns you to post in plain text (for less confusion and a better reception). --- Jeff NewmillerThe . . Go Live... DCN: