similar to: Sample size

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

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

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

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

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

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

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

Hello, R for Dummies. How can I exclude the first 1000 values of a vector (length 12000)? More generally all the values up to the ith? Thanks for your help, Dr Olivier Charansonney Cardiologue Centre Hospitalier Sud-Francilien, Corbeil-Essonnes, France [[alternative HTML version deleted]]

Dear R users, I have just upgraded R from 2.11 to 2.12 on Ubuntu 9.04 (see more informations at the end) from the cran apt-get repository. One of the new things concerning the install.packages() function is stated here : install.packages() and remove.packages() with lib unspecified and multiple libraries in .libPaths() inform the user of the library location used with a message

Hi Since 2006, I used in my departement a Samba solutions based on NT4 style PDC, 2 BDCs and some files servers, desserving one hunderd persons. The backends for passwords and idd are a master and two slave openldap. Now i have to integrate a much larger organisation, an University Hospital, running with Windows AD. For political reasons, I should not maintain DCs anymore, but I will still

Hello, I have just completed first experiments in using R, especially creating and using ARIMA models, e.g. # create a model as in example(arima) fit <- arima(USAccDeaths, order = c(1,1,1),seasonal = list(order=c(1,1,1))) # use the model to generate a prediction dp<-predict(fit, n.ahead = 24) plot(dp$pred) # view the prediction Now I created a combined chart of the original and

Hi, My departemental PDC/BDCs will be removed and the remaining linux file servers and workstations will joined the institutional Windows AD as member servers and workstations. I have the rights to add workstations and servers to the AD, but I will lost users administration. All that is Ok. Nervertheless, groups administration at the departement level is still usefull for me and my

Hello, It seems that domain user can access share when they are specified in "valid list" but not when "valid list" use local group definition. First if added the domain user "duser" to the group "lgroup" in /etc/group Then i defined a samba share and add the domain user "duser in the "valid list" [lgroup] comment =