similar to: mgcv 'bam' : prediction levels for random effects

Hello, I'm trying to use the bam() function in the R mgcv package for a large set of grouped binary data. However, I have found that this function does not take data in the format of cbind(numerator, denominator) on the left hand side of the formula. As an example, consider the following dat1 <- data.frame(id=rep(1:6, each=3), num=rbinom(18, size=10, prob=0.8), den=rbinom(18, size=5,

Dear useRs, This is the first time I post to this list and I would appreciate any help available. I've used the excellent mgcv package for a while now to investigate geographical patterns of language variation, and it has has always worked without any problems for me. The problem below occurs using R 2.14.0 (both 32 and 64 bit versions in Windows and the 64 bit version in Unix) and mgcv (both

Dear useRs, I reran an analysis with bam (mgcv, version 1.7-17) originally conducted using an older version of bam (mgcv, version 1.7-11) and this resulted in the same estimates, but much lower standard errors (in some cases 20 times as low) and lower p-values. This obviously results in a larger set of significant predictors. Is this result expected given the improvements in the new version? Or

Do you have the mgcv package installed (I think it's part of the standard distro, though) /loaded? ziP is there, not in BAM. Other than that, sorry, no clue. 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 Mon, Nov 6,

Dear all, I am trying to use 'bam' to run the following generalized additive model: ### m <- bam(result ~ factor(city) + factor(year) + lnpopulation + s(lnincome_pc) + ,data=full_df,na.action=na.omit,family=ziP(theta = NULL, link = "identity",b=0)) But getting the following error: ### Error in bam(result : extended families not supported by bam The documentation for

Hi all: I've got a data of ".bam" which is created from my partner under linux sysyem. My system is window xp, and I wanna know how to read the .bam file. Many thanks! My best [[alternative HTML version deleted]]

Hello, I am trying to update the mgcv package on my Linux box and I keep getting an "Out of memory!" error. Does anyone know of a fix for this? Below is a snippet of the message that I keep getting: Thank you. Geoff ** R ** inst ** preparing package for lazy loading ** help *** installing help indices >>> Building/Updating help pages for package 'mgcv' Formats:

Hello everyone. I am doing a logistic gam (package mgcv) on a pretty large dataframe (130.000 cases with 100 variables). Because of that, the gam is fitted on a random subset of 10000. Now when I want to predict the values for the rest of the data, I get the following error: > gam.basis_alleakti.1.pr=predict(gam.basis_alleakti.1, +

Hello R users, I have a quick question I was hoping to get your input on. I am new to R and the smooth statistical regression world, and am trying to wrap my mind around the issues concerning using splines for mixed effect modeling. My question is the following: in the ?gamm? function, generalized additive mixed models can be estimated by including random components. These can be explicitly

mgcv: Constructing probabilities for binomial GAM with crossed random intercepts and factor by variable Hello, (I'm sorry if this has been discussed elsewhere; I may not have been looking in the right places.) I ran a binomial GAM in which "Correct" is modelled in terms of the participant's age and the modality in which the stimulus is presented (written vs spoken).

Hi, I also had some problems using chan_oh323 together with gnugk. * <-> gnugk <-> h323-phone When I called the phone and hang up, befor the phone was picked up, the h323-phone continued ringing. The same, when the h323- and some sip-phones were called, and the sip-phone picked up the call first. (It is annoying, when you are talking to someone at the phone and the phone on the

Hello, I'm trying to set up an alternate user with a different inbox using passwd-file and userdb extra fields and the extra field settings don't seem to be taking effect. I've used auth_debug and mail_debug and I'm unclear how I can further troubleshoot this. Is it possible to have inbox=yes in one namespace globally and then override this for one user by setting inbox=no for

Hello, I'm having a problem renaming and sorting the underlying factor of a ggplot2 based plot. Here's my code: ---8<---------- > delta <- ggplot(subset(data, Model==c("dyn", "dl4", "dl3")), aes(x=Problemsize, y=Fitness)) + geom_smooth(size=1, color="black", fill=alpha("blue", 0.2))+ geom_point(size=0.5,

Hi I'm missing something here but I cannot figure out what. What I can see is that the same code works when I load it via source(...) yet fails when I execute it after loading the package I have built (which includes the code. Below is a transcript of my R session. First I load the code from a file, using source(). Then I execute it fine. Then I remove the function object, I load the

Hi, Is the stars command in the base package or do I need to download? I am looking to make star/radar/spider charts. Thanks! Britt Britt Aronovich Marketing Analyst BAM (Brooklyn Academy of Music) ... Peter Jay Sharp Building 30 Lafayette Ave. Brooklyn, NY 11217-01486 ... P: 718.724.8038 E: baronovich@BAM.org<mailto:baronovich@BAM.org> BAM.org Facebook.com/BAMstage

I have * working on my Sony Vaio PCG-FX120 laptop. I am trying to get * to recognize my internal PCI Intel modem as an FXO port. I have modified wcfxo.c in order to identify the PCI modem properly. Based on output from dmesg, wcfxo didn't recognize the modem until I inserted the proper vendor and device IDs into wcfxo.c and re-compiled. Note that the error message from "modprobe

Hi all, We have been using SuSE7.2 and samba (client and server) without any problems for a few months. Now, we migrated some workstations (samba clients) to SuSE8.0 and we get segfaults while accessing files mounted through samba. The segfault happens on the CLIENT side. Here is our network configuration: - SuSE7.2 server with default samba server package. - Few SuSE7.2 workstations

Hi, I have a irregularly spaced time series dataset, which reads in from a .csv. I need to convert this to a regularly spaced time series by filling in missing rows of data with NAs. So my data, called NtuMot, looks like this (I've removed some of the additional rows for simplicity).... ELEID date_time height slope 1 2009-06-24 00:00:00

Hi all, I have installed a TE110P in a BOX but when I load zaptel module I can't see any device in /proc/zaptel. And led of the card is green. My zaptel.conf is the next: span=1,1,0,ccs,hdb3,crc4 bchan=1-15,17-31 dchan=16 loadzone=es defaultzone=es and cat /proc/pci throguh next: PCI devices found: Bus 0, device 0, function 0: Host bridge: Intel Corp. 82845 845

Hi, The problem is most likely, you need to call a R CMD BATCH with your arguments and the R-script inside of a shell script that you submit to your qsub. Unfortunately we don't use qsub anymore so can't test it, but it should be as follows: R-script eg. test.R: > ##First read in the arguments listed at the command line > args=(commandArgs(TRUE)) > > ##args is now a list of