similar to: Problem with effects package

Dear R People, I'm working on a plot function to produce the graph shown below. One of the features my supervisor does not like is that the labels in the Y- axix are shown vertically. Is there a way to change that to horizontally ? This is the axix() function I used : axis(2,at=seq(nstrat),labels=snames) The labels (snames) are: SEV1:HU:C, SEV1:MU:C, SEV1:MU:L, SEV1:RT:C and SEV1:RT:L,

Hi friends - a javascript question.... Does the prototype object add some way for a class to find out the variable name of its instance variable? Or does javascript have some built in way to do this? This is confusing to describe - but if I define a variable as some object: var myVariable= new Widget(''fdfa''); can some built-in method inside the Widget class be

Hello Group, I could use an advice on how SEXP handles work. My aim is to implement a system where I initially set a few global variables that are used for communication between C and R code. Then I do some work with R code and periodically call a function of my own that will update the system state. Such a design is useful for many purposes (for GUIs to name one). I am not entirely sure that R

Hi, Could anyone tell me what is wrong: > length(unique(mydata$myvariable)) [1] 2 > and in t-test: (...) Error in t.test.formula(othervariable ~ myvariable, mydata) : grouping factor must have exactly 2 levels > I re-checked the code and still don't get what is wrong. Moreover, there is some strange behavior: /1 It seems that the error is vulnerable to NA'a, because it

Hi - How can I 'manually' reproduce the results in 'pred1' below? My attempt is pred_manual, but is not correct. Any help is much appreciated. library(splines) set.seed(12345) y <- rgamma(1000, shape =0.5) age <- rnorm(1000, 45, 10) glm1 <- glm(y ~ ns(age, 4), family=Gamma(link=log)) dd <- data.frame(age = 16:80) mm <- model.matrix( ~ ns(dd$age, 4)) pred1 <-

Hi, I am quite new to R - but quite expierience in programming. Nonetheless I have some problemes in accessing information of the "summary" object. Here is what I do: model <- lm ( y ~ myVariable ) "summary(model)" gives me an object which has a lot of information about the regression. Now I'd like to access programmatically the level of significance which is

Often it is useful to keep a "codebook" to document the contents of a dataset. (By "dataset" I mean a rectangular structure such as a dataframe.) The codebook has as many rows as the dataset has columns (variables, fields). The columns (fields) of the codebook may include: ? variable name ? type (character, factor, integer, etc) ? variable label

I generated maps with the function symbols (graphics). These are basic maps generated with : symbols(x,y,circles=myvariable) where x et y are spatial coordinates corresponding to replicates of "myvariable". I would associate legend to this kind of maps, is it possible? Regards, Marion. -- Marion Jacquot Laboratoire de Chrono-environnement UMR UFC/CNRS 6249 USC INRA Universit? de

Hello everybody. I want to make a menu with the incoming calls, I want that when someone calls the Asterisk play a record (in gsm) and them the caller must choose a option (1,2 or 3). if he choose 1 it will redirect to 101 extension if he choose 2 it will redirect to 102 extension if he choose 3 it will redirect to 103 extension my extensions.conf is this one: [default] exten =>

Estimados amigos, Estoy familiarizándome con los modelos lineales generalizados en R. Estoy interesado en realizar un análisis lig linear y me gustaría saber cuáles son o como extraer los valores correspondientes al chi cuadrado en el análisis para cada grupo y para las interacciones. Desde ya muchas gracias y disculpas si la pregunta es muy básica, adjunto los comandos que estoy utilizando. Si

Estimados amigos, Estoy familiarizándome con los modelos lineales generalizados en R. Estoy interesado en realizar un análisis lig linear y me gustaría saber cuáles son o como extraer los valores correspondientes al chi cuadrado en el análisis para cada grupo y para las interacciones. Desde ya muchas gracias y disculpas si la pregunta es muy básica, adjunto los comandos que estoy utilizando. Si

Estimados amigos, Estoy familiarizándome con los modelos lineales generalizados en R. Estoy interesado en realizar un análisis lig linear y me gustaría saber cuáles son o como extraer los valores correspondientes al chi cuadrado en el análisis para cada grupo y para las interacciones. Desde ya muchas gracias y disculpas si la pregunta es muy básica, adjunto los comandos que estoy utilizando. Si

Hi! I would like to perform an F-Test over more than one variable within a generalized mixed model with Gamma-distribution and log-link function. For this purpose, I use the package mgcv. Similar tests may be done using the function "anova", as for example in the case of a normal distributed response. However, if I do so, the error message "error in eval(expr, envir, enclos) :

Ted Harding said: > I can get the estimated RRs from > RRs <- exp(summary(GLM)$coef[,1]) > but do not see how to implement confidence intervals based > on "robust error variances" using the output in GLM. Thanks for the link to the data. Here's my best guess. If you use the following approach, with the HC0 type of robust standard errors in the

Howdy, I am trying to run a Cobol Compiler through wine. I am running RH7.3 with a current wine release. My directory structure is /usr/share/wine-c/COBOL/. I have been successful in running other applications through wine, but not ones that require a PATH= entry. My problem is: it's looking for an entry in the PATH environment variable that includes 'C:\COBOL' which stores

Dear all, I am running a logistic regression and this is the output: glm(formula = educationUniv ~ brncntr, family = binomial) Deviance Residuals: Min 1Q Median 3Q Max # ???? ????? ?? ???????? -0.8825 -0.7684 -0.7684 1.5044 1.6516 Coefficients: Estimate Std. Error z value Pr(>|z|) (Intercept) -1.06869 0.01155 -92.487 <2e-16 *** brncntrNo

I am getting a repeated error when I try to run a logistic regression in R 2.8.1 >(glm(prop1~x1,data=glm1,family=binomial("logit"),weights=nt1)) Error in model.frame.default(formula = prop1 ~ x1, data = glm1, weights = nt1, : invalid type (list) for variable 'x1' x1 is multistate categorical (3 categories). 2 of the categories have 12 observation, one has 9. Is this what

R version 1.3.0 OS: SunOS 5.7, but I think the same problem occurs with Windows An incorrect formula seems to be used to calculate the pearson residuals for a generalized linear model with a binomial response. Here is a simple program which gives (a) the pearson residuals calculated directly, (b) the pearson residuals from glm, and (c) the deviance residuals from glm. The first and last

This is a question regarding the 'regTermTest' function in the 'survey' package. Imagine Z as a three level factor variable, and code ZB and ZC as the two corresponding dummy variables. X is a continuous variable. In a 'glm' of Y on Z and X, say, how do the two test specifications test.terms = c("ZB:X","ZC:X") # and test.terms = ~ ZB:X + ZC:X in

Hi, all I have a question on mgcv and ns. Now I want to compare the results from glm, gam and ns. Take a simple model y~x for example. glm1 = glm(y~x, data=data1) gam1 = gam(y~s(x), data=data1) ns1 = glm(y~ns(x),data=data1) In order to confirm the result from glm1 is consistent to those from gam1 and ns1, I want to define degree=1 in mgcv and ns. I am wondering if there is somebody can give me