similar to: Fw: Using if, else statements

Hi all, Again I have searched the net and so on, without finding an answer to this surely simple problem. A short bit of code would be appreciated. I have a object named `data' with the following column headings. Date, maxitemp, minitemp, admissions, d.o.w. Where d.o.w. is day of the week, written "Sun" "Mon" etc. I just need to scale the Monday admissions by

Hi all, Thanks for the help with my previous post. I have just two more questions for the minute. I think I said in a previous post that I like to use the terminal, i.e. run rterm.exe. On exiting the terminal, I am asked if I want to save the workspace. If I hit y (yes), the workspace is just saved as .rdata in my working directory, does anyone know how I can name it directly from the

Estimados: Yo tengo un archivo de texto Unicode y no funciona read.table con encoding Latin-1 o UTF-8 estando en un entorno Windows. Si no es con ninguno de los anteriores, ¿como lo importo? ¿Cómo puedo saber que codificación tiene? Gracias, Sebastián. El día 22 de abril de 2010 07:00, <r-help-es-request en r-project.org> escribió: > Envíe los mensajes para la lista R-help-es a

Hola todos: acabo de actualizar R a R 2.10.1, y cuando actualice Rcmdr (Versión 1.5-4) me cargo la versión en inglés, no me acuerdo que hice antes para tener Rcmdr en español (trabajo con windows XP). alguien podría recordarme lo que hay que hacer? Gracias Gabriela. ______________________________ Lic. María Gabriela Cendoya Magíster en

Dear all, I am having difficulties using the seemingly-simple function lag. I have a dataframe with several weather variables (maxitemp, windspeed, rainfall etc), and the response variable (admissions). The dataset is fairly large (1530 observations). I simply want to model the response against a lag of a couple of the explanatory variables, say maxitemp and rainfall. I would like to look at

Hi all R users. I'm just starting to use R, and I'm still a little lost. I'd like to analyse some data following the path analysis aproach. I'd like to know if there is some package for implement this kind of analysis. Many thanks Rodrigo Sala Posgrado en Producción Vegetal E.E.A. INTA Balcarce Ruta Nacional 226 km 73.5 (7620) Balcarce, Buenos Aires, Argentina [[alternative

Hi R-Helpers: I?m trying to use the function GLMM from lme4 package, (R-1.8.1, Windows 98),and I get the following error: > pd5 = GLMM(nplant~sitio+ + fert+ + remo+ + sitio:fert+ + remo:sitio+ + remo:fert+ + remo:fert:sitio + data=datos, + family=binomial, +

Hola, mi pregunta es la siguiente: Tengo un modelo genético aditivo, donde los individuos estan emparentados. Quiero que la matriz de parentesco sea la matriz de coeficientes de las correlaciones entre genotipos. En SAS, con PROC MIxed podia agregarla como un archivo y usarla en la sentencia random...No encuentro como hacerlo en R, alguna idea?? En resumen, puedo imponerle a los efectos

Hi R-Helpers: I?m trying to fit the same linear model to a bunch of variables in a data frame, so I was trying to adapt the codes John Fox, Spencer Graves and Peter Dalgaard proposed and discused yesterday on this e-mail list: for (y in df[, 3:5]) { mod = lm(y ~ Trt*Dose, data = x, contrasts = list(Trt = contr.sum, Dose = contr.sum)) Anova(mod, type = "III") } ## by John Fox or for

Hi: Does R have a function as gsorth is SAS, that perform a the Gram-Schmidt orthonormal factorization of the m ?n matrix A, where m is greater than or equal to n? That is, the GSORTH subroutine in SAS computes the column-orthonormal m ?n matrix P and the upper triangular n ?n matrix T such that A = P*T. or any other version of Gram-Schmidt orthonormal factorization? I search the help, but I

Hi: I'm using function lmer from package lme4, and I get this message: " There were 12 warnings (use warnings() to see them)" So I checked them: Warnings 1 to 11 said: 1: optim returned message ERROR: ABNORMAL_TERMINATION_IN_LNSRCH in: "LMEoptimize<-"(`*tmp*`, value = structure(list(maxIter = 50, ... and Warning 12 said: 12: IRLS iterations for glmm did

This is the Department of Computer Science at the University of Illinois at Urbana-Champaign. Your inquiry has been received. If you wish to reach another department or campus you must contact them directly. We can not forward your message to other departments or campuses. Please do not send email to bcicone@cs.uiuc.edu. This address is no longer valid for inquiries. SPECIAL NOTES: DEADLINE FOR

I''m trying to get my routes to work properly. Here''s how I''m defining one of my routes: map.downloads '':category/downloads'', :controller => ''downloads'', :action => ''list_downloads_for_category'', :requirements => {:category => DownloadableFileCategory.category_regex} And here''s the method

I fit the model fit<-lm(thresh~cond*Ne) where thresh is the reponse cond is a factor with levels a, b, and c Ne is a continuous indep var I think of this full model as having three lines: thresh as a function of Ne for each condition. Thus we have slopea, slopeb, slopec, inta, intb, intc. lm output my params ------------------------- (Intercept) inta condb intb - inta condc

Dear users I got this msg accessing s share: server:~# smbclient '\\localhost\Adm1' -U gabriela Added interface ip=192.168.0.254 bcast=192.168.0.255 nmask=255.255.255.0 Password: Domain=[EQUATORIAL] OS=[Unix] Server=[Samba 2.0.5a] smb: \> put lixo2 ERRDOS - ERRnoaccess (Access denied.) opening remote file \lixo2 The smb.conf says [Adm1] comment = Adm-Financeiro path

I'm trying to fit a Cox proportional hazards model to some hospital admission data. About 25% of the patients have had at least one admission, and of these, 40% have had two admissions within the 12 month period of the study. Each patients has had one of 4 treatments, and one of the treatment groups has had no admissions for the period. I used:

The print method for lme *prints out* the fitted correlation matrix for the random effects. Is there any way to get these values as an object in R? I have examined the components of the lme object (called "junk" in the example below) and the components of summary(junk) without finding these numbers. (How I did this: I dumped the entire lme object to a text file and then used egrep to

Hello, I have a general help question that some of you might be able able to help. I would like to design a study (a two group comparison) based on reduction on events(say hospital admissions). In a previous study hospital admission rate of 140 admissions per 72 patients (over a 4 month period) has been observed. That is rate is about 1.9. In order to see the admission rate reduction of 50% over

R / S Programmer - New York Description: Kaplan Test Prep & Admissions is looking for a R / S Programmer to join their research based in New York City. We are looking for highly motivated individuals to work in a customer-focused environment. This is a unique opportunity to develop with a leading educational company with a diverse series of statistical applications in education/learning,

Estimados me ha surgido una duda con un simple stem and leaf el mismo es el siguiente: > stem.leaf(Datos$ej1.21) 1 | 2: represents 0.12 leaf unit: 0.01 n: 40 LO: 0.72 0.85 3 10 | 9 11 | 4 12 | 4 5 13 | 7 7 14 | 710 (*) 9 15 | 18 16 16 | 3447899 (6) 17 | 045599 (***) 18 18 | 2568 14 19 | 23710 (**) 10 20 | 389 7 21