similar to: COZIGAM removed from CRAN

Hi there, My question is about the 'theta' parameter in specification of a NB GAMM. I have fit a GAM with an optimum structure of: SB.gam4<-gam(count~offset(vol_offset)+ s(Depth_m, by=StnF, bs="cs")+StageF*RegionF, family=negbin(1, link=log), data=Zoop_2011[Zoop_2011$SpeciesF=='SB',]) However, this GAM shows heterogeneity in the

Hi there, I'm using the lattice package to create an xy plot of abundance vs. depth for 5 stages of barnacle larvae from 5 species. Each panel of the plot represents a different stage, while different loess smoothers within each panel should represent different species. However, I would like depth to be on the y-axis and abundance to be on the x-axis, because this is more intuitive as an

Hi I have the following barchart to which I want to add error bars. library(lattice) barchart(Change~fTreat,groups=Process,change, auto.key=list(points=FALSE,rectangles=TRUE), panel=function(x, y,...){ panel.barchart(x,y,origin = 0,...); panel.abline(h=0,col="black",...); } ) I have tried

Dear All I am attempting some power analyses, based on simulated data. My experimental set up is thus: Bleach: main effect, three levels (control, med, high), Fixed. Temp: main effect, two levels (cold, hot), Fixed. Main effect interactions, six levels (fixed) For each main-effect combination I have three replicates. Within each replicate I can take varying numbers of measurements (response

I making an xyplot and the y label is too long and needs to be in two rows, but when I brake it there is a huge gap between the last text string and the expression, and I can't get rid of it. Any ideas? Data: structure(list(Temp = c(8L, 8L, 8L, 8L, 8L, 8L, 12L, 12L, 12L, 12L, 12L, 12L), CO2 = c(380L, 380L, 380L, 750L, 750L, 750L, 380L, 380L, 380L, 750L, 750L, 750L), Treat = structure(c(3L,

Greetings everyone, I am running MDS on a very large dataset (12 x 25071 - 12 model runs with 25071 output values each), and also on a very much reduced version of the dataset (randomly select 1000 of the 25071 output values). I would like to look at similarities/dissimilarities between the 12 model runs. When I use metaMDS on the full dataset, I get a warning message: Warning message: In

I making an xyplot and the y label is too long and needs to be in two rows, but when I brake it there is a huge gap between the last text string and the expression, and I can't get rid of it. Any ideas? Data: structure(list(Temp = c(8L, 8L, 8L, 8L, 8L, 8L, 12L, 12L, 12L, 12L, 12L, 12L), CO2 = c(380L, 380L, 380L, 750L, 750L, 750L, 380L, 380L, 380L, 750L, 750L, 750L), Treat = structure(c(3L,

Dear All I am having problems to install the package COZIGAM. Despite I have the update version of R (3.0.2) when I try to load the package I cannot find it in the list of available packages. If I try to install it from a local zip file I do get this error message: Error in read.dcf(file.path(pkgname, "DESCRIPTION"), c("Package", "Type")) : cannot open the

Hello, I'd like to implement a regression model for extremely zero-inflated continuous data using a conditional approach, whereby zeroes are modelled as coming from a binary distribution, while non-zero values are modelled as log-normal. So far, I've come across two solutions for this: one, in R, is described in the book by Gelman & Hill (http://www.amazon.com/dp/052168689X), where

Dear List, I posted this in R-mixed and did not receive any feedback. I might post it in the wrong place. I re-post in R-help and hope to receive any suggestions and\or thoughts regarding data analysis. The objective of the study is to investigate effects of soil properties on insect outbreaks. There are four study fields (or sites). Data were collected from 1996 through 2009. Below is

 Hello everyone, Apologies in advance, as this is partially a stats question and partially an R question.  I have been using a GAM to model the activity level of bats going into and coming out from a forested edge.  I had eight microphones set up in a line transect at each of eight sites, and I am hoping to construct a model for each of 7 species.  My count data has a reverse J-shaped skew and

Hello, I would like to retrieve data stored in MySQL database, so I installed RMySQL package. I can successfully connect with the my database using the following code > dvr<-dbDriver("MySQL") > con2<-dbConnect(dvr,group="exbardiv") > mysqlDescribeConnection(con2) <MySQLConnection:(972,0)> User: mmorag Host: localhost Dbname: exbardiv

I posted this last week but haven't heard anything. I'm not sure if this is because nobody knows the answer (can't believe that!) or I'm missing something obvious in the documentation and people are thinking "Read The Fine Manual". Whatever the reason, if anyone has any insights into this problem I'd be very grateful for their comments. We're using Samba 3.0.23b

Classification: NOT PROTECTIVELY MARKED Solaris 9 Samba 3.4.5 I know this isn't the sort of query that gets much response but I'd be really grateful of any advice people can offer. I'm getting really fed up with Samba as I've never been able to make it work properly. Either I'm missing something basic (probably) or it just doesn't behave in the way I think it should!

Hi, I am trying to calculate a distance matrix on a binary data frame using dist.binary() {ade4}. This is the code I run and the error I get: > sjlc.dist <- dist.binary(as.data.frame(data), method=2) #D = (a+d) / (a+b+c+d) Error in if (any(df < 0)) stop("non negative value expected in df") : missing value where TRUE/FALSE needed I don't know if the problem are the

Hi Can anybody tell me where i can download or find out what the latest version of samba for Tru64 V5.1a Pk 3. Currently we are running on Samba version 2.2.7 pre 1 and i need to know if any upgraded versions will run with no issues on any version of windows. Many Thanks Jason ********************************************************************** The information in this E-Mail is

We are encountering some problems compiling a working version of Samba 3.4.1 with ADS support on Solaris 9. We've run autogen.sh and configure successfully (config.log attached), followed by make and make install. However, despite having openldap and Heimdal Kerberos (1.2.1) installed, and including the configure switches --with-ads --with-ldap --with-krb5=/usr/heimdal, when we try a net join

I was curious to read the following in the "Guidelines for Rd files" at http://developer.r-project.org/Rds.html: "For consistency, aim to use British (rather than American) spelling. (NB: British spelling often uses -ize as in `capitalize'. The view of spell -b of British spelling is a flawed American one. There are English/Scottish differences

I am currently using the function daisy in package cluster to create a dissimilarity matrix because my multivariate dataset contain missing data and variables of various types including factors, symmetric and asymmetric binary and quantitative. This is a step prior to using pco within ecodist. There is a warning which comes twice ">NAs introduced by coercion" I've used

Speaking of optimization and speeding up R calculations... I mentioned last week I want to speed up calculation of generalized inverses. On Debian Wheezy with R-2.15.2, I see a huge speedup using a souped up generalized inverse algorithm published by V. N. Katsikis, D. Pappas, Fast computing of theMoore-Penrose inverse matrix, Electronic Journal of Linear Algebra, 17(2008), 637-650. I was so