search for: scrucca

...tution: > plot(1:10, main = substitute("Monotonic Multigamma run (" theta == t1 * ").", list(len = len, t1 = t1))) which works fine. How can I have more than one value substitution, together with mathematical notation and text? Thanks in advance for any reply. Luca Scrucca +-----------------------------------------------------------------------+ | Dr. Luca Scrucca | | Dipartimento di Economia, Finanza e Statistica | | Sezione di Statistica tel. +39-075-5855226 | |...

...-matrix and the response: lqs(y ~ x1 + x2, method="lms", nsamp="exact", x.ret=T, y.ret=T) Error in model.frame(formula, rownames, variables, varnames, extras, extranames, : variable lengths differ Does anyone know what's going on? Regards, with thanks in advance, Luca Scrucca +-----------------------------------------------------------------------+ | Dr. Luca Scrucca | | Dipartimento di Scienze Statistiche tel. +39 - 075 - 5855278 | | Universita' degli Studi di Perugia fax. +39 - 075 - 43242...

...".") { owd <- getwd() on.exit(setwd(owd)) setwd(path) } to be FALSE and then non changing of current directory is done. Is this the desired behavior or some bug fix is required? Luca +-----------------------------------------------------------------------+ | Dr. Luca Scrucca | | Dipartimento di Economia, Finanza e Statistica | | Sezione di Statistica tel. +39-075-5855226 | | Universit?? degli Studi di Perugia fax. +39-075-5855950 | | Via Pascoli - C.P....

...ain the above ordering (in the example an anti-clockwise ordering) such that I can use polygon() to connect the points? I searched previous messages but I did not find any relevant to this problem. Thanks, Luca +-----------------------------------------------------------------------+ | Dr. Luca Scrucca | | Dipartimento di Scienze Statistiche tel. +39 - 075 - 5855278 | | Universit佮 degli Studi di Perugia fax. +39 - 075 - 5855950 | | Via Pascoli - C.P. 1315 Succ. 1 | | 06100 PERUGIA (ITAL...

...lem ? I hope to be able to install the RXLisp package, I loved the dynamic graphics capabilities of Xlisp-Stat and it would be very nice to have them available within R. Thank you very much for your help. Luca +-----------------------------------------------------------------------+ | Dr. Luca Scrucca | | Dipartimento di Scienze Statistiche tel. +39 - 075 - 5855278 | | Università degli Studi di Perugia fax. +39 - 075 - 5855950 | | Via Pascoli - C.P. 1315 Succ. 1 | | 06100 PERUGIA (ITAL...

...ed") In this case you must simply use an S chart: > obj <- qcc(a,type="S") Some sparse notes: - don't use a[,] to select a whole matrix, simply a, i.e. qcc(a) - don't use <<- for assignment (it has another meaning...), <- is enough I hope this help. Luca Scrucca -------------------------------------------------- Luca Scrucca Dipartimento di Economia, Finanza e Statistica Sezione di Statistica Universit? degli Studi di Perugia Via A. Pascoli, 20 06123 PERUGIA (ITALY) Tel. +39-075-5855233 Fax: +39-075-5855950 E-mail: luca at stat.unipg.it Web page: http...

Dear list, This is a question about overdispersion and the ML estimates of the parameters returned by the glm.poisson.disp (L. Scrucca) and glm.nb (Venables and Ripley) functions. Both appear to assume a negative binomial distribution for the response variable. Paul and Banerjee (1998) developed C(alpha) tests for "interaction and main effects, in an unbalanced two-way layout of counts involving two fixed factors, when data...

Dear R People: Is there a Quality Control in R book that would be accessible for undergraduates, please? I'm teaching a QC course in the fall semester and would be happy to have a Use R book or something like that, please. Thanks in advance, Erin -- Erin Hodgess Associate Professor Department of Computer and Mathematical Sciences University of Houston - Downtown mailto: erinm.hodgess at

Hi All, I am running a scatter plot and trying to add a best fit line. I use an abline function, but get no line drawn over the points. I also get no error. I arm using V 2.10.0 on Windows 7. Here is my code, including the SAS transport file import: require (foreign) require (chron) require (Hmisc) require (lattice) clin <- sasxport.get("y:\\temp\\subset.xpt") attach(clin)

Hi, I am performing glm with binomial family and my data show slight overdispersion (HF<1.5). Nevertheless, in order to take into account for this heterogeneity though weak, I use F-test rather than Chi-square (Krackow & Tkadlec, 2001). But surprisingly, outputs of this two tests are exactly similar. What is the reason and how can I scale the output by overdispersion ?? Thank you,

Hello, I have a problem regarding calculation of Cumulative Incidence Function. The event of interest is failure of bone-marrow transplantation, which may occur due to relapse or death in remission. The data set that I have consists of- lifetime variable, two indicator variables-one for relapse and one for death in remission, and the other variables are donor type (having 3 categories), disease

violating.runs I read from the news cran.r-project.org/doc/Rnews/Rnews_2004-1.pdf that the criteria for the violating is 5 but 1)I cannot find "5" in the code of the function. Where is the "5" ? 2)What is the easiest way to change it ? 3)Is there any more criterias made somewhere ? Yours sincerelly, Tommi Viitanen

Dear R-help, I'm very sorry to ask 2 questions in a week. I am using the package 'crr' and it does exactly what I need it to when I use the dataset a. However, when I use dataset b I get the following error message: Error in drop(.Call("La_dgesv", a, as.matrix(b), tol, PACKAGE = "base")) : system is computationally singular: reciprocal condition number =