similar to: bootstrap pValue in DClusters

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

Hi, As Frank M. Howell noticed (and probably other users), the Windows binaries for DCluster I put in my web page are not working... I have compiled the source code again and know it does. Please, download it again, and sorry for the inconvenience. The URL is http://matheron.uv.es/~virgil/Rpackages/DCluster DCluster is a package that implements some methods for the detection of spatial

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

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 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

When I try to load the DCluster package I get the following error message (R 2.5.1; PowerBook G4; Mac OS X 10.4.10) (Load needed package) Lade nötiges Paket: spdep (error: package spdep could not been loaded) Fehler: Paket 'spdep' konnte nicht geladen werden (additional warning message) Zusätzlich: Warning message: (There is no package spdep in:.....) es gibt kein Paket 'spdep'

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

Hi, I have just put the source code (and a zip file for Window$ users) of package DCluster in my website: http://matheron.estadi.uv.es/~virgil/Rpackages/DCluster/ DCluster is a package that contains routines for the detection of spatial clusters of diseases (Openshaw's GAM, Besag and Newell, Kulldorff and Nagarwalla, Stone's Test and others). A full description can be found here:

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

Dear R users, I am glad to announce the release of a new package for R: DCluster. It implements a number of methods for the detection of clusters of disease. It includes methods to test Poisson extra-variation (Chi-square test and Potthoff-Whittinghill's test), Spatial Autocorrelation (Moran's I and Geary's c), general clustering (Whittermore's statistic and Tango's

Dear R users, I am glad to announce the release of a new package for R: DCluster. It implements a number of methods for the detection of clusters of disease. It includes methods to test Poisson extra-variation (Chi-square test and Potthoff-Whittinghill's test), Spatial Autocorrelation (Moran's I and Geary's c), general clustering (Whittermore's statistic and Tango's

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

Hello, I'm a student in biology, writting from Strasbourg, France. I need to use the logiciel R to analyse biological results. I have ADE4 since 2001, but technology's now largely evoluted ! That's why I'm just trying to get the new version of : - ADE4_1.3-3.zip ; - R 2.0.1.ter ; - DCluster_0.1-3.zip. However, I don't know HOW TO INSTALL R on my computer ? I don't

hi, I'm looking to compare two area under ROC values for different classifiers on the same data -- is there an r function to do this? Thanks! [[alternative HTML version deleted]]

Hi, I'm trying to evaluate a Monte Carlo p-value (using truncgof package) on a left truncated sample. >From an empirical sample I've estimated a generalized pareto distribution parameters (xi, beta, threshold) (I've used fExtremes pkg). I'm in doubt on what of the following command is the most appropriate: Let: x<-sample t<-threshold xt<-x[x>t] xihat<-gpdFit(x,

Dear R people, I have a naive question: after fitting "lm" to a data, I can't extract the pvalue corresponding to a specific covariate in a direct way. Could anyone give me a hint? Thank you very much. Frank

Hi all, I have some t-test values, and I am trying to obtain the associated p-values. Is 'pt' the right command? I wonder why 1) it returns different values for x and -x, and 2) how to obtain a 2-sided p-value. example [R version 2.0.1, WinXP]: #if t=2.23 (df=10), the expected p-value is 0.05 for 2-sided and 0.025 for 1-sided t-test pt (2.23,10) [1] 0.9750788 pt (-2.23,10) #or pt