similar to: mgcv::gam error message

Hi all, My question concerns 2 error messages; one in the gam package and one in the mgcv package (see below). I have read help files and Chambers and Hastie book but am failing to understand how I can solve this problem. Could you please tell me what I must adjust so that the command does not generate error message? I am trying to achieve model selection for a GAM which is required for

Or with the Bioconductor IRanges package: df <- with(input, DataFrame(station, year=IRanges(from, to), record)) expand(df, "year") DataFrame with 24 rows and 3 columns station year record <character> <integer> <character> 1 07EA001 1960 QMS 2 07EA001 1961 QMC 3 07EA001 1962 QMC 4

Hi, A benchmarking study with an additional (data.table-based) solution. Enjoy! ;) Cheers, Denes -------------------------- ## packages ########################## library(dplyr) library(data.table) library(IRanges) library(microbenchmark) ## prepare example dataset ########### ## use Bert's example, with 2000 stations instead of 2 d_df <- data.frame( station =

On 11/29/2017 04:15 PM, T?th D?nes wrote: > Hi, > > A benchmarking study with an additional (data.table-based) solution. I don't think speed is the right benchmark (I do agree that correctness is!). For the R-help list, maybe something about least specialized R knowledge required would be appropriate? I'd say there were some 'hard' solutions -- Michael (deep

Hi Martin, On 11/29/2017 10:46 PM, Martin Morgan wrote: > On 11/29/2017 04:15 PM, T?th D?nes wrote: >> Hi, >> >> A benchmarking study with an additional (data.table-based) solution. > > I don't think speed is the right benchmark (I do agree that correctness > is!). Well, agree, and sorry for the wording. It was really just an exercise and not a full

On 11/26/2017 08:42 PM, jim holtman wrote: > try this: > > ########################################## > > library(dplyr) > > input <- tribble( > ~station, ~from, ~to, ~record, > "07EA001" , 1960 , 1960 , "QMS", > "07EA001" , 1961 , 1970 , "QMC", > "07EA001" , 1971 , 1971 ,

try this: ########################################## library(dplyr) input <- tribble( ~station, ~from, ~to, ~record, "07EA001" , 1960 , 1960 , "QMS", "07EA001" , 1961 , 1970 , "QMC", "07EA001" , 1971 , 1971 , "QMM", "07EA001" , 1972 , 1976 , "QMC", "07EA001" , 1977 ,

Dear Mr. Simon Wood, dear list members, I am trying to fit a similar model with gam from mgcv compared to what I did with BayesX, and have discovered the relatively new possibility of incorporating user-defined matrices for quadratic penalties on parametric terms using the "paraPen" argument. This was really a very good idea! However, I would like to constraint the coefficients

dplyr may have something for this, but in base R I think the following does what you want. I've shortened the name of your data set to 'd'. i <- rep(seq_len(nrow(d)), d$YEAR_TO-d$YEAR_FROM+1) j <- sequence(d$YEAR_TO-d$YEAR_FROM+1) transform(d[i,], YEAR=YEAR_FROM+j-1, YEAR_FROM=NULL, YEAR_TO=NULL) Bill Dunlap TIBCO Software wdunlap tibco.com On Sat, Nov 25, 2017 at 11:18 AM,

On Tue, Sep 15, 2015 at 12:31:25PM +0900, Timo Sirainen wrote: > On 12 Sep 2015, at 04:45, Bradley Giesbrecht <pixilla at macports.org> wrote: > > > > Is the Dovecot Object Storage plugin still available for purchase? > > > > It's available, but at least for now we're only selling it to big customers. Well how big? Would a univeristy with 25.000 users

To David W.'s point about lack of a suitable reprex ("reproducible example"), Bill's solution seems to be for only one station. Here is a reprex and modification that I think does what was requested for multiple stations, again using base R and data frames, not dplyr and tibbles. First the reprex with **two** stations: > d <- data.frame( station =

Hi; I would like to fit a contiuous right skewed dist. to a set of data. So, I want to use MLE or Robust estimation of the parameters, then use A-D or K-S test or Robust alternative for the goodness-of -fit. Also I would rank the fitted distributions according to the p-value of the test. Does anyone know how I do this in R. any backages of Robust estimation and for goodness-of-fit? Thank you

Hi, I tried to load a .RData object on unix system using R, it gives error: Error: restore file may be empty -- no data loaded In addition: Warning message: file 'junk3.RData' has magic number '' Use of save versions prior to 2 is deprecated This happens only for using MY user account for the Unix system. I tried to use a friends's user account to load the same data

Dear List, I have two questions about how to do predictions using lrm, specifically how to predict the ordinal response for each observation *individually*. I'm very new to cumulative odds models, so my apologies if my questions are too basic. I have a dataset with 4000 observations. Each observation consists of an ordinal outcome y (i.e., rating of a stimulus with four possible

Bert wrote ... Bill's solution seems to be for only one station. No, it works for any number of stations. Bill Dunlap TIBCO Software wdunlap tibco.com On Sun, Nov 26, 2017 at 11:10 AM, Bert Gunter <bgunter.4567 at gmail.com> wrote: > To David W.'s point about lack of a suitable reprex ("reproducible > example"), Bill's solution seems to be for only one

Hi all, I need to run several simple linear regressions at once, using the following data. Response variables: Bird species (sp 1, sp2, sp3...spn). Independent variable: Natprop - proportion of natural area. covarate: Effort = hours). One single linear regression would be: lmSp1 <- lm(sp1~ natprop + effort). However, I need to run this linear regression for all bird species that I have

I have a returned tibble of station operational record similar to the following: > data.collection # A tibble: 5 x 4 STATION_NUMBER YEAR_FROM YEAR_TO RECORD <chr> <int> <int> <chr> 1 07EA001 1960 1960 QMS 2 07EA001 1961 1970 QMC 3 07EA001 1971 1971 QMM 4 07EA001 1972 1976 QMC 5

> On Nov 25, 2017, at 11:18 AM, Hutchinson, David (EC) <david.hutchinson at canada.ca> wrote: > > I have a returned tibble of station operational record similar to the following: > >> data.collection > # A tibble: 5 x 4 > STATION_NUMBER YEAR_FROM YEAR_TO RECORD > <chr> <int> <int> <chr> > 1 07EA001 1960

Dear Sirs! My name is Mohamad TALIAA and I am supposed to be the BUDDY for Viola Rossini at the Univeristy of Vienna! The ERASMUS Bureau at Vienna just provided me with a wrong email address! I reported to them the mistake but I just was told to look for it maybe on the internet or change the extension at the end. Now I just came to your site as i looked for Viola Rossini at google.at I also

On 25 Sep 2015, at 17:56, Hans Morten Kind <Kind at adm.uib.no> wrote: > > On Tue, Sep 15, 2015 at 12:31:25PM +0900, Timo Sirainen wrote: >> On 12 Sep 2015, at 04:45, Bradley Giesbrecht <pixilla at macports.org> wrote: >>> >>> Is the Dovecot Object Storage plugin still available for purchase? >>> >> >> It's available, but at