similar to: discrepancy in the result of R and SAS on same data in logistics regression

Any thoughts on the following I'd be most grateful - I'm sure there is an easy and quick way to do this but I'm having a mental block this evening. Essentially, I'm trying to replace missing data in my dataset with reference values based on age and sex. So an example dataset is set.seed(1) X = data.frame(age=rnorm(10, 50, 10), sex=rbinom(10, 1, 0.5), A=rnorm(10), B=rnorm(10))

Hola, A lo mejor la solución pasa por ver si Statgraphics puede exportar su fichero a otro formato de cualquiera de los que R sí que puede leer. No sé si además del fichero, tienes acceso a Statgraphics.... Saludos, Carlos Ortega www.qualityexcellence.es El 12 de junio de 2013 13:26, Roberto <rv.tec@zoho.com> escribió: > Hola, es la primera lista en la que participo, espero no meter

Dear R users, Is it possible to have an overlap histogram plot? For example: stuff <- data.frame(Mode = c("Land", "Air"), AgeGroup = c("Young", "Adult", "Old"), Value = sample(1:300)) histogram( ~Value | AgeGroup * Mode, data = stuff, auto.key = TRUE) Instead of having 2 * 3 panel, I want to have just 3 panel, overlapping the "Mode"

I'm new to R - and trying to create a plot similiar to the spider plot at http://www.statgraphics.com/eda.htm#radar . I can't figure out several things... most of which I would think would be straightforward.... How can I change the lines for each series plotted instead of creating a filled area? How can I get the labels for each of the radial axes at the outside of the plot

With R 1.8.1 running in Fedora Core 1 Linux, I am having some trouble recoding and ploting some factor variables. First, can I give you some example data? Here is a column of names for age groups: agegroups <- c( "15-19", "20-24", "25-29","30-34", "35-39",

Hola, es la primera lista en la que participo, espero no meter mucho la pata y aprender rápido, se ruega comprensión. Simplemente me gustaría saber la forma de importar datos desde un archivo ".sf3' de Statgraphics. Creo que ese es un programa bastante difundido, pero no lo tengo y en principio no quiero tenerlo porque no es libre. Pero encuentro referencias a todo (piInfo, Minitab,

Hola Jorge: Disculpa la tardanza pero me han tenido liado en otros menesteres. Yo no soy ni muchísimo menos un experto ni en series temporales ni en R de hecho retomo el tema después de muchos (demasiados) años aparcado y dedicándome a labores de programación pura y dura. Respecto a la diferencia de resultados con R y Statgraphics, no conozco el proceso de selección del modelo ARIMA que hace

Hi, The problem is that I would like to do an interaction plot with intervals based on Tukey's honestly significant difference (HSD) procedure, but I do not know how to do it in R. I have 3 factors "A", "B" and "C" and a response variable "response". I would like to study a model where there are main effects and second order interaction effects.

Hello! I am trying to run a GLMM using LME4, and keep getting the warning message: "In mer_finalize(ans) : false convergence (8)" I am quite new to R, and in looking into this thus far, it appears that there are a variety of reasons why this might occur, such as needing to standardize some parameters or if all subjects in one combination of parameters all have the same outcome. I also

Hi, I have what I believe is a repeated-measures dataset that I'm trying to analyze using lme(). This is *not* homework, but an exercise in my trying to self-teach myself repeated-measure ANOVA for other *real* datasets that I have and that are extremely similar to the following design. I'm fairly sure the dataset described below would work with lme() -- but it'd be great if anybody

Hello. I am attempting to duplicate a negative binomial regression in R. SAS uses generalized estimating equations for model fitting in the GENMOD procedure. proc genmod data=mydata (where=(gender='F')); by agegroup; class id gender type; model count = var1 var2 var3 /dist=NB link=log offset=lregtm; repeated subject=id /type=exch; run; Since my dataset has several observations for

Hello! I would like to add names to a table. The table presents median values, numbers of values (n) and percentage(%) for a given agegroup. Thus I would like to add a label above these three variables for each category. But if I try to use colnames I get a message telling me that the number of labels must be equal to the number of columns in the table. I hope someone could understand what I mean!

Hello everyone, I use Bob Muenchen's approach for reading in "in-stream" (to use SAS parlance) delimited data within a script. This works great: mystring <- "id,workshop,gender,q1,q2,q3,q4 1,1,f,1,1,5,1 2,2,f,2,1,4,1 3,1,f,2,2,4,3 4,2, ,3,1, ,3 5,1,m,4,5,2,4 6,2,m,5,4,5,5 7,1,m,5,3,4,4 8,2,m,4,5,5,5" mydata <- read.table( textConnection(mystring),

Hi R users: I want to know if any of you had used Stata or Statgraphics. What are the advantages and disadvantages with respect to R on the following aspects? 1. Statistical functions or options for advanced experimental design (fractional, mixed models, greco-latin squares, split-plot, etc). 2. Bayesian approach to experimental design. 3. Experimental design planing options. 4. Manuals

Hi all, Completely new to this list and completely new to R so, please, accept my apologies in advance if any question is too basic or out of place. Thanks. After reading some docs, googling around and searched the mail lists I've only found spare references to accessing databases. Among these references, few are related with JDBC and none of them include both MySQL and JDBC access. Our

Hi, I'm trying to teach myself R, and by the way, re-learning statistics using Crawley's "Statistics: an introduction using R". I've reached the regression chapter, and when it deals with non-linear regresion using the nls library I face the following problem: I follow the steps--- >deer<-read.table("c:\\temp\\jaws.txt",header=T) ---data available at

Hello all, I've been able to use R very successfully to run simple statistics and generate the plots I require. I've been evaluating Sweave, and have hit upon a small problem that I don't seem to be able to workaround. Sweave runs very well for single file latex documents, but I have a complex thesis made up of several parts and chapters. These are arranged with a master latex file

Hello I just installed wine without any problems, and executed a statistics programs called statgraphis. Everything went fine, I could do the calculations but then I wanted to save the result, the way I am used to it when I use it in windows. However when saving the result in a rtf file, I am asked where and I chose drive Y, which seems to be my home directory. Now some strange messages appear on

Hello there! I have a really basic question. I have extracted some data from a big database. I have ordered the data in to different agegroups in the rows. Now I want to to have labels of the columns which is going to be divided into men , women and totals. How do I give names to a column? My data looks like this for the moment; 0-39 12 34 43 40-49 14 23 23 50-59 16

I'm creating a stacked bar chart using the likert command in the HH package. My data are in a data frame, with two numeric variables and a categorical variable, I can't get likert to use the column containing the categorical variable as a my y axis label. Here is a quick example: library(HH) #my data are: