similar to: Intercept by level in GLM

Dear R-list, I think it's really a newbie question but ... I try to model age-period-cohort models with polynoms in each effect. I have for each level of age and period some cases -k- and persons-years -py-. Models are A_xP_y: log(k_age,period/py_age,period)=f(polynom(x),age)+f(polynom(y),period). For exemple I try to adjust an A2P3 :

Dear all, I want to use GraspeR (on R 1.6.0) and have to install tcltk library. When I try to do this, I get back this error message: Error in firstlib(which.lib.loc, package) : TCL_LIBRARY is not set Error in library(pkg, character.only = TRUE) : .First.lib failed But I have a tcl/tk package (Active TCL 8.4.1.0) correctly install on my PC (Win95) and in the autoexec.bat the

Dear R-list, I try to adjust GLMM on incidence cancer data. Without random effect, in GLM the command is, for example with sex effect, glm(Observed~sex+offset(log(Expected)),family=poisson) because the observed are Poisson distribued with parameter Expected*incidence rate. But know I want to introduce random effect (for example spatial effect) and it seems to me that the "offset" does

Dear R-list, I wish a very happy new year and send you a little question: I have different variables which names are m1, m2, m4, .., m10, ... and want to obtain for example mean of each of them without typing each mean(m1), mean(m2), .... What is the solution for decomposing names in mXX? With my best regards, EA. ============================================ Erik-André SAULEAU SEAIM Hôpital du

Hello R-List I use DClusters package (I work in a cancer regestry). I have 2 questions about it: 1-how is it possible to get back the bootstrap pValue? I mean the pValue of the calculated statistic with respect of the distribution of this statistic under the null hypothesis. 2-how is it possible to test an overdispersion in the poisson model? for choosing a best model I need this mesure of

R-users, i have some time series and use HEGY seasonal UR tests in the uroot package. But i have problems with interpretation of t and F statistics. Is there any critical values for these tests in R? is anybody ever write a bootstrap of these tests? thanks for advance, with all my bests, erik. ================================= Erik-André SAULEAU SEAIM Centre Hospitalier 87, Avenue

Dear R-list, I try to use StructTS for some univariate time series, with a type="BSM". I have two questions: 1-in the "fitted" value of the object, it seems to me that for getting back fitted values, one have to add level, slope and seasonal (or "level" contains already the two other ones) 2-i am really circumspect about the results of "KalmanForecast"

Dear R-list, i have a problem, in the framework of simulations, i want to vectorize for earning time: a variable, say X, has values on intervals and an other variable, say Y, has values on other intervals. For example Inf Sup X 0 2 1 2 4 2 4 6 3 and Inf Sup Y 1 3 1 3 5 2 5 7 3 i want to create a matrix like this Inf Sup

Can someone point me to the appropriate functions for fitting multiple statistical models that are coupled to each other. The data are measurements of salinity s and temperature t at stations id and pressures at p as well as surface elevations h at stations id. The problem is, for any new station, to estimate s at all p, given t at all p and given h. If h is ignored, the for each p, there

Sorry to bother you all with what is probably a simple question. I am attempting my first go at a simple AGI application using PHP (Getting Asterisk to SAY PHONETIC ABC). I have dabbled with PHP but I am by no means a professional standard developer. My script seems to execute ok, and I can see asterisk playing the sounds but my phone goes from ringing to busy, and I don't hear the

Sorry to re-post this but I'm sure it's something simple that someone has found before. To summarise: Dial plan answers the phone AGI script executes AGI debug in console show phonetics ABC - However no audio at all on the phone and this step is less than 1 second. Dial plan Busy Phone hangs up. Total time less than a second. This is such a simple AGI script do I need

Hello, Apologies if this question has already arised, hope you can help me to the find the solution to this or point the place to look at. I have a multidimensional piece-wise regression linear problem, i.e. to find not only the regression coefficients for each "interval" but also the beginning and ends of the intervals. To simplify it to the one dimensional case and two intervals,

Hi I am still having problems with using the stl function, when I read the csv file into R into a file called tkr and use dim(tkr) the result is 132 1 which is fine. When coerce it into a trime series using ts either: tstkr <- ts(t(tkr), deltat=1/12) or tstkr <- ts(c(tkr), deltat=1/12) and use the stl function I get the following error: Error in

Un texte encapsul? et encod? dans un jeu de caract?res inconnu a ?t? nettoy?... Nom : non disponible Url : https://stat.ethz.ch/pipermail/r-help/attachments/20070220/cb89f312/attachment.pl

Hi R-helpers, I'm new to time series modelling, but my requirement seems to fall just outside the capabilities of the arima function in R. I'd like to fit an ARMA model where the variance of the disturbances is a function of some exogenous variable. So something like: Y_t = a_0 + a_1 * Y_(t-1) +...+ a_p * Y_(t-p) + b_1 * e_(t-1) +...+ b_q * e_(t-q) + e_t, where e_t ~ N(0, sigma^2_t),

Hi, i want to create an object which contains different classes: for example i have some time series and test if ARIMA models are best than HoltWinters models: for each of my time serie i want to collect in an unique object which model was the best: for some it will be an HoltWinters class and for some other an Arima class. is there any solution?? thanks in advance, with my best, erik sauleau

Hello, Here is a patch for protocol_subnet.c with two modifications : - in tunnelserver mode, tinc must check subnets in the ".../hosts/owner" config file, not in "c->config_tree" (which is the configuration of the meta-connection from which we receive the ADD_SUBNET message). - this checking can be made before the check of the owner, especially before any

I have been experimenting with different possibilities for an "introduction" page for my packages. That is, a good place to tell users about the most important things in a package, and where to start. Recently there was a discussion about this, and a suggestion to use <foo>-package.Rd, and also a function that generates a skeleton document. My problem with this suggestion is

Hi, When looking into a Flac file with an hexadecimal editor, I observe that the different tags (date, tittle, etc...) are, one time all displayed in capital letters (DATE=), on an other time all displayed in minus letters (date=), or...a mix of both : some tags in capital letters, others in minus letters. For tagging (and encode/decode) I use the Flac Frontend listed on your site. When testing

All: I cvsup'ed earlier this evening and am still able to reproduce this panic at will. Command line, panic backtrace, dmesg, ldd output, and mplayer version follow. If more information is needed, just let me know. % mplayer foo.mov [a quicktime file] [plays for awhile, and then panics] # gdb -k -s /usr/obj/usr/src/sys/LICHEN/kernel.debug -e /var/crash/kernel.2 -c /var/crash/vmcore.2