similar to: About ade4 and overlaying points

hello, I work since a few time on R and i wanted to know how to obtain the Wald chi square value when you make a binary logistic regression. In fact, i have the z value and the signification but is there a script to see what is the value of Wald chi square. You can see my model below, Best regards, S??verine Erhel [Previously saved workspace restored] > m3 = glm(reponse2 ~ form +

Hello I am working on radio tracking data, with a short programme someone gave me and ... which should, supposedly, work ... In this programme, there is the function : getareahr(kern, levels = 95). But i cannot find any 'getareahr' in R ... could anyone help me? thanks! Agn??s

Greetings, I'm posting this OT query here because of out very international membership! In the French sentence "Les taux de tirage sont calcul??s de mani??re ?? ce que la dispersion soit inf??rieure ?? 5 % dans chaque strate." it would seem intended that the "dispersion" is to be calculated in a specific way (unstated) -- otherwise, how to ensure that it shall be

R-masters, I have a problem that I have been working on for a while and it seems that there may be a simple solution that I have yet to figure out, so I thought that I would venture to post to the help list. Let's say there was a data.frame with three vectors, two that are factors identifying the data, and one that holds the frequency of occurrence (the events are binary, yes or no). I

Hi all, Does anybody know the easiest way to import excel data into R? I copied and pasted the excel data into a txt file, and tried read.table, but R reported that Error in read.table("data_support.txt", sep = " ", header = T) : more columns than column names Thanks! Ling

I am trying to model data in which subjects are followed through time to determine if they fall, or do not fall. Some of the subjects fall once, some fall several times. Follow-up time varies from subject to subject. I know how to model time to the first fall (e.g. Cox Proportional Hazards, Kaplan-Meir analyses, etc.) but I am not sure how I can model the data if I include the data for those

Hi everyone, I'm posting again, since I haven't got an answer (yet :( ). According to R help, manova does not support the inclusion of the Error() term in the formula call. I have repeated measures data for two dependent variables, so this means I can't account for subject variance in time?. Any lights? Thanks in advance, Rodrigo Abt, Department of Economic Studies, SII, Chile.

Dear R-helpers, I have found a slightly annoying problem when trying to plot lines on graphs. I first created my data using tapply, thus:- y1=as.vector(fit1$coef$random$id) x1=tapply(o1,id,median,na.rm=T) x2=tapply(o2,id,median,na.rm=T) #then I plot the data, thus:- plot(x1[x2==0],y[x2==0]) #if I now try to fit the linear regression, R 'hangs up'

Hi, Is it possible to recreate "smoothed" data sets in R, by performing a PCA and then reconstructing a data set from say the first 2/3 EOFs? I've had a look in the help pages and don't seem to find anything relevant. Thanks in advance, Laura Laura Quinn Institute of Atmospheric Science School of Earth and Environment University of Leeds Leeds LS2 9JT tel: +44 113 343 1596

Does anyone know how to extract residuals in lmer? Here's the error I get: > crop.lme=lmer(response~variety*irrigation*pesticide+(1|rep)+(1|rep: pesticide)+(1|rep:pesticide:irrigation), crop.data) > qqnorm(crop.lme) Error in qqnorm.default(crop.lme) : y is empty or has only NAs > resid(crop.lme) NULL Thanks! --Jake

Hi, I have a 3 panel xyplot with different variables in the y axis. I'm trying to insert a newline after "Width (cm)," in the ylab argument as in the example below. My goal is to have the y axis label broken into two lines, split after the string just mentioned. plotfun <- function() { fakedf <- data.frame(A = sample(1:100, 50), B = rnorm(50),

Dear Fellows, How can I do to proced a step wise regression in R, if it´s possible ? Thanks, Walmir __________________________________________________________________________ Acabe com aquelas janelinhas que pulam na sua tela. AntiPop-up UOL - É grátis! http://antipopup.uol.com.br/ [[alternative HTML version deleted]]

Hello, Does anyone know of a "slick" way to get R to convert a string which is all upper case to a string where the first letter in each word is upper case and all others are lower case? I suspect the solution is to begin by parsing the string, convert the appropriate letters to upper and lower case using "toupper" and "tolower", and then to paste the pieces back

Dear list, I am brand new to R and using Dalgaard's (2002) book Introductory Statistics with R (thus, some of my terminology may be incorrect). I am fitting regression models and I want to use Hurvich and Tsai's AICC statistic to examine my regression models. This penalty can be expressed as: 2*npar * (n/(n-npar-1)). While you can obtain AIC, BIC, and logLik, I want to impose the AICC

Dear useRs We have a data-set (comma delimited) with 12Millions of rows, and 5 columns (in fact many more, but we need only 4 of them): id, factor 'a' (5 levels), factor 'b' (15 levels), date-stamp, numeric measurement. We run R on suse-linux 9.1 with 2GB RAM, (and a 3.5GB swap file). on average we have 30 obs. per id. We want to aggregate (eg. sum of the measuresments under

Hi, I'm trying to do a box-whisker plot of two columns of a data frame, a list of category names in one column vs. some numerical values in the other. The plot itself works fine, but only a few points of the x-axis ( the category names ) are labelled. I think that this is because the category names are too long. Is there any way to force R to label each x-axis value, preferably at a

Hello, I am a beginner in R programming and recently heard about this mailing list. Currently, I am trapped into a simple problem for which I just can't find a solution. I have a huge dataset (~81,000 observations) that has been analyzed and the final result is in the form of 0 and 1(one column). I need to write a code to process this column in a little complicated way. These 81,000

What is the easiest way to change within vector of strings each letter after a space into a capital letter? E.g.: c( "this is an element of the vector of strings", "second element" ) becomes: c( "This Is An Element Of The Vector Of Strings", "Second Element" ) My reason to try to do this is to get more readable abbreviations. (A suggestion would be to

Exactly where does R search for foo.R if I type source("foo.R")? Only from current working directory (same as getwd()), from all directories specified by e.g. $PATH? Thanks. Daehyok Shin

Exactly where does R search for foo.R if I type source("foo.R")? Only from current working directory (same as getwd()), from all directories specified by e.g. $PATH? Thanks. Daehyok Shin