similar to: generic name of variables

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

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

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 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 R-listers, I have to explain "y", a continous variable, with two variables "t" and "p" which are "factor variables": t is [1:12] and p is [1:18] (a value for y for each combination t*p). OK, then I try to adjust models like: y=a_p + b*t (I mean as many intercepts as levels of p) or y=a_p + b_p*t (as many intercepts and coefficients for t as levels of

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

Why will this simple statement not work? I think I am following the documentation for if(cond) statements, and I have tried wrapping the cons.expr and alt.expr in {}, I get the same error. There is no example in the help file, and this is not covered in the Introduction to R, SimpleR or other tutorials I have looked into. mxx=max(cpx_list$nMV); mxy=max(trend_list$nMV); if (mxx>mxy) mxy=mxx

Dear R users,I'm a new user of R. I'm trying to do a LM test an there is this type of error: Error in t(mX) %*% mX : requires numeric/complex matrix/vector arguments. To be clear I write down the code in which mY ( 126,1 ) mX (126,1) mZ(126,1) are matrix. LMTEST <- function(mY, mX, mZ)#mY, mX, mZ must be matrices!#returns the LM test statistic and the degree of freedom{iT =

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"

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

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

Dear Researchers, I wrote two function to fit a circle using noisy data. 1- the fitCircle() is derived from MATLAB code of * zhak Bucher* from the link http://www.mathworks.com/matlabcentral/fileexchange/5557-circle-fit/content/circfit.m 2- the CircleFitByPratt() from MATLAB code of *Nikolai Chernov *from the link

Hi all, I have developed a full featured Windows IAX phone based on LIBIAX library . It is now in a prerelease version (0.9.0) and you can download it for free from my web page: http://www.laser.com/dante or http://www.geocities.com/tdanro Some of the features are: - registering with Asterisk PBX; - can use any audio device as ring device (including PC speaker), independent of the play device;

FYI -------- Original Message -------- Subject: RE: AOL Mail Date: Wed, 25 Jun 2008 13:47:06 -0500 From: Dustin Davis <dustin@txls.com> To: "'Dale Schroeder'" <dale@BriannasSaladDressing.com> References: <419cf800ca3991a24f4b9829f3806aba@briannassaladdressing.com> <145401c8d604$ccfb1780$66f14680$@com>

Hi, I'm new here and I would like to know what's the best place to ask for some structural questions or help? If that's the right place, I'll send another email with the details and my question. Thanks in advance! Sven -- Sven SAULEAU +33 6 28 69 51 44

Of course Sven, questions regarding LLVM belong here (clang questions should be directed to cfe-dev). Developers on this list are usually very helpful when it comes to answering questions that are clearly stated and show that the person asking has done a reasonable amount of investigation and is just stuck. It usually helps to start with a succinct description of what it is you are trying to

https://bugzilla.samba.org/show_bug.cgi?id=8666 Summary: --debug=all9 fail Product: rsync Version: 3.1.0 Platform: All OS/Version: All Status: NEW Severity: normal Priority: P5 Component: core AssignedTo: wayned at samba.org ReportedBy: chris at onthe.net.au QAContact: rsync-qa at

Hello, I am almost new in R, so perhaps my question will be silly. I try to use R for analyzing capture-recapture data in epidemiology. A cancer registry has different sources of patients. We know in each list, patients already known in all other list. The aim is to use capture-recapture models for estimating the number of patients unknow of all the sources. Because no order in sources, one