similar to: Offset for Poisson GLM

Hello, I am synthesising some sales data over a twelve month period, and then trying to use the "predict" function, firstly to generate a thirteenth month forecast with upper and lower 95% confidence limits. So far so good But what I then want to do is add the upper sales value at the 95th confidence limit to the vector of thirteen months and their respective sales to create a

I have realised that I should have used "detach" before attaching another dataframe, but even when I do this it's still giving me lots of lines, rather than just one: My code: m<-runif(1,0,1) m mres<-m*(seq(1,12)) mres ssd<-rexp(1,1) ssd devs<-rep(0,length(mres)) for(i in 1:length(mres)){devs[i]<-rnorm(1,0,ssd)} devs plot(-10,-10,xlim=c(1,24),ylim=c(0,20000))

Hello, This is probably very simple but I am trying to select a region from the OIv2 SST dataset that crosses 0 longitude. The longitude variable is arranged from 0 to 360 and I get an error message when trying to go from 330 to 30 for example. The data are from a NOAA netCDF file. *temp = open.ncdf("c:/documents and settings/theoni/desktop/sstoi_v2.nc") print(temp) y1 = get.var.ncdf(

Hello to all, I am trying to use the spherical wavelets package on my altimetry data. I am having some problems now and then. Currently I am having an error of "error in matrix(0,KK,JJ) : too many elements specified" there is a thread of possibly same error (https://stat.ethz.ch/pipermail/r-help/2005-January/064308.html) but it seems like it doesnt really help me. More

When I compile my C fibonacci example fib.c with 'clang -O3 -c -emit-llvm -o fib-clang.bc fib.c&& llvm-dis fib-clang.bc' I get fib-clang.ll that has some degree of inlining in it. But when I get an equivalent to fib.c file fib.ll and run it through opt with the command 'llvm-as fib.ll&& opt -O3 fib.bc -o fib-opt.bc&& llvm-dis fib-opt.bc' resulting

Hello All, I am looking for some help with was is probably very simple: I want the y-axis title '$B&D(B13C$B"s(B VPDB' and similarly for nitrogen on a secondary axis. I am able to get the delta notation and per mil notation separately but one messes up the other with the inclusion of quotations: ylab=expression("{delta}^13*C \211 VPDB") gives the per mil just fine but

Signed-off-by: Christoph Hellwig <hch@lst.de> Index: linux-2.6/fs/ocfs2/cluster/tcp.c =================================================================== --- linux-2.6.orig/fs/ocfs2/cluster/tcp.c 2007-05-06 13:51:17.000000000 +0200 +++ linux-2.6/fs/ocfs2/cluster/tcp.c 2007-05-17 15:00:14.000000000 +0200 @@ -261,14 +261,12 @@ out: static void o2net_complete_nodes_nsw(struct o2net_node

R-helpers, I am hoping to find someone who uses both R and program Stata for GLMs. I am a beginner R user, finding my own way through; learning code etc. at the same time as learning the statistics I need to complete my project. What I have is the code from Stata and am trying to reproduce the same analysis in R - my program of choice. . glm count md ms rf sg, family(poisson)

G'Day, I'm a little confused why the c function has the followng effect on classes - is this a feature ? My workround [ class(cc) <- c("POSIXt", "POSIXct") ] seems to do the job. Many thanks Bernie McConnell "R version 1.3.1, 2001-08-31" on NT > aa <- as.POSIXct("2001-09-23") > bb <- as.POSIXct("2001-09-24") > cc

R-users, Can an offset term be included in a Poisson model? I get an error message when trying that: >r3o <- glm(tax ~ areal + offset(o), family=poisson) Error in (if (is.empty.model(mt)) glm.fit.null else glm.fit)(x = X, y = Y, : inner loop 1; can't correct step size In addition: Warning message: Step size truncated due to divergence in: (if (is.empty.model(mt))

Hi, Sorry if this is an obvious one, but is there a simple way to modify the lm function to test whether a slope coefficient is significantly different from 1 instead of different from 0? Thanks, Martin Martin Biuw SEA MAMMAL RESEARCH UNIT Gatty Marine Laboratory School of Environmental and Evolutionary Biology University of St Andrews Fife, Scotland KY16 8LB UK phone +44 (0)1334 462630 fax

Hello, Once again, I posted a message without a subject line. Sorry.... here is the question again. Is there a simple way to modify the circ.mean function in the CircStats package to include a vector of weights to obtain a weighted average angle? Thanks! Martin -- Martin Biuw Sea Mammal Research Unit Gatty Marine Laboratory, University of St Andrews St Andrews, Fife KY16 8PA Scotland Ph:

Hello all, Does anyone know of any already written functions for carrying out stochastic dynamic programming in R? Thanks! Martin -- Martin Biuw Sea Mammal Research Unit Gatty Marine Laboratory, University of St Andrews St Andrews, Fife KY16 8PA Scotland Ph: +44-(0)1334-462637 Fax: +44-(0)1334-462632 Web: http://smub.st.and.ac.uk

Hello, I''m trying to fit a linear mixed effects model of the form: lme(y ~ x * Sex * Year, random=x|subject) where Sex and Year are factors with two and three levels respectively. I want to compare the fixed effects for each level to the overall mean, but the default in R is to compare to the first level. This can be changed by adding the term -1 to the righthand side of the model

Hello, I'm using the nls function and would like to increase the number of iterations. According to the documentation as well as other postings on R-help, I've tried to do this using the "control" argument: nls(y ~ SSfpl(x, A, B, xmid, scal), data=my.data, control=nls.control(maxiter=200)) but no matter how much I increase "maxiter", I get the following error

I myself am not that well-versed in c and opengl - especially opengl - I can read and write basic c code but that's my limit. The problem only appears when using an offset that's not zero.If you want to debug the issue you could very well use the dolphin-emu and apitrace. You can find one apitrace here http://markus.members.selfnet.de/dolphin/dolphin-emu.trace (~48M) The only 2 reasons I

Hi I am trying to join a out-of-the box win2k3 AD domain controller it's been + forest prep for r2 domain samba Version 3.3.2 on Freebsd 6.3-RELEASE openldap-sasl-2.3 heimdal 0.6.3 adserver = AD DC server i installed (win2k3 box) domain = my domain name /etc/resolv.conf search domain.net nameserver adserver contents of /usr/local/etc/smb.conf [global]

Some how I''m getting NameError: uninitialized constant when accessing a child models that inherits form the parent model. However, if I access the parent model first, then the child model is resolved. What am I doing wrong? My steps are below... Thanks, Francis. > ./script/generate model Blah type:string class CreateBlahs < ActiveRecord::Migration def self.up create_table

Thomas, OK, thanks for the update. I don't have a working machine with NV graphics that supports UBOs. Could you compile and run a test program for me if I wrote it? Though it may appear to not work for UBO offsets other than 0, there may be another pattern here, such as multiple of page size or something crazy like that. If we can find more information, it may help other developers. Patrick

Hi, It looks like KDE's file manager Dolphin doesn't support MTP under CentOS 7. When I plug in an MTP device like my Nikon camera or an Android phone, Digikam shows its content OK, but Dolphin seems to ignore it. For now, my custom CentOS 7 + KDE desktop meets about 99 % of my needs. Now I've decided to tackle the remaining 1 %. Yesterday for example, I managed to correct the