similar to: index of dispersion

Dear list, I calculated the kappa-accuracy for two differnt classifications. How can I test now the two kappa-value for significance? thanks, Mark ..................................................................... Markus Schwarz Wissenschaftliche Mitarbeiterin Eidg. Forschungsanstalt WSL Forschungsprogramm Musterland Z??rcherstrasse 111 CH-8903 Birmensdorf Telefon +41-44-739 22 87 Fax

Hello, I have two classifications. How can I compare the overall-accuracy of these classifications to each other? Is there a possibility within R to test if the achieved overall-accuracy for the classifications are differing significantly? Additionaly, are the McNemar-Test and Broker-Test implented in a package of R? Thank in advance for your help, Markus

I have the following problem: I measured co2 on 6 paired sites (one grubbed and one non-grubbed fence per site -> grubbing = treatment). These measurements I repeated 15 times over 2 years. So, now my problem is how to analyze these data. I tried the following model: mod1_CO2<-lme(co2~treatment+time,random=~1|site,data=CO2_t1_t15) I think with the random effect I included the paired

Dear R colleagues, I am trying to create a new column in a data frame, which converts values and NA's from another column into binary format. Essentially I need the NA's to become 1 and the rest to be 0. The code I wrote is returning the following error message: Error in if (mort[i, 4] != NA) mort[i, 8] <- 0 else if (mort[i, 4] == : missing value where TRUE/FALSE needed

I' ve seen that several people are looking for a function that creates a barplot with an error indicators (I was one of them myself). Maybe you will find the following code helpful (There are some examples how to use it at the end): # Creates a barplot. #bar.plot() needs a datavector for the height of bars and a error #indicator for the interval #many of the usual R parameters can be set:

Dear all, I would like to multiply two matrixes with the different dimension column by column. Let make an example: If I have two matrixes "X" and "Y"as follow: X<- matrix(1:12, nrow=4, ncol=3, dimnames=list(c("A","B","C","D"), c("stage1","stage2","stage3"))) Y<- matrix(1:28, nrow=4, ncol=7,

Dear R users, I have a dataframe with two columns: first column is date data (e.g. 1/1/2000 with character format: daily data from 1/1/1970 till 31/12/2003) and second column is temperature value. Now I'd like to calculate mean for each month in a year (i.e. May 2001, June 1997) and mean for each month in all of years. As the number of days in some months is different from others I could

Dear r-users! I would like to cross two data frame which have the same row number but different in the number of column. Can anybody help me for this case ? Thanks a lot in advance -------------------------------------------------------------------------------- Majid Iravani PhD Student Swiss Federal Research Institute WSL Research Group of Vegetation Ecology Z?rcherstrasse 111

Dear R-users I want to save a list with characters a point plot and a Venn diagram in a single pdf page. I am successful to do this when I use a character list and two point plots. However when I try to replace the first point plots with my Venn diagram (built with Vennerable package, compute.Venn() and plot.Venn()) the Venn plot will not position at the right place in the pdf. I guess there

Dear list members, I have got another problem. I imported an SPSS file with the Memisc package using the following commands: mz <- spss.system.file("myspssfile.sav") mz <- subset(mz,select=c( bsex,balt,xurb,dtaet,kartab,bgeb,boseit,bgeblan,xnuts2,kausb,xerwstat, asbper,asbhh,ajahr,aquartal,bstaat,xwieoft,gew1,apkz,bpkzm,bpkzv)) Afterwards I checked the measurements of the

Hello Cracks, we want to import 10 x 365 files containing each 7 variables with 1211 rows. The filenames end with numbers from 1 to 365. We tried a for loop to import using read.table but we couldn't create the r objects within the loop. We want each file in an own r-objekt (matrix). Does anyone know what can be done? Thanks, Felix -- Felix Tiefenbacher Eidg. Institut f?r Schnee- und

Hello everybody 1. If there is/ever will be a function to perform "duncan multiple range test" please inform me at once. 2. Is there a way to create a grouping as duncan does from TukeyHSD output? My experimental design contained 62 genotypes, so the pairwise comparison is not that usefull and clear to look at. How would I do that? Is there an other test that does give grouping as

I want to use the plot(model) function to generate Tukey-anscomb and Q-Q plots of a lm(). I manage to change all labels but the main one which apparently is neither main or sub. So far I have tried as par setting: cex (changes symbol size within the plot), cex.main (no effect), cex.sub (no effect) cex.lab (changes label size), cex.axis (changes axis label size). What I would like to change is the

Dear all, I've got a list L <- list(L1 = list (foo = "bar" , SL = NULL ) , L2 = list ( foo = "bar" , SL = list (SSL1 = list (DF = data.frame(val = 21, foo = "bar") , DFOO = list(foo = "foo", bar = "bar") ) ,

Dear R-users with a bit of grief I had to repeat an extensive analysis because I did not suspect (and therefore did not read the documentation) that round was implemented as "for rounding off a 5, the IEC 60559 standard is expected to be used, 'go to the even digit'", resulting in round(1.5) = 2 round (2.5) = 2. As a non-mathematician I am both puzzled and intrigued by

Dear List adding a table to a plot using addtable2plot(plotrix) does not seem to work when using a logarithmic axis. The table is then reduced to one line. Is there an argument to indicate that a log-scale is used, or an alternative to add a bunch of information to a plot? I'm using R 2.6.1 on WinXP SP2. thanks Markus -- Markus Didion Wald?kologie Forest Ecology Inst. f.

Hi I am using nlsList from the nlme package to estimate group specific coefficient sets to get a first impression of the structure of coefficients in relation to potential predictor variables. In some cases I got coefficients in other there are conversion problems due to bad start values I guess, resulting in typical error messages as known from nls: Error in nls(formula = formula, data = data,

Hi, is there a possibility to calculate the different "types" of residuals directly using the multinom function from MASS as it is possible for the functions gam, glm using type="deviance" or "working" or "pearson" or "response"? I tried it but got always the "response" type, I guess. thanx Matthias

Hi, I am searching for methods to compare regression models with an ordered categorical response variable (polr versus multinom). The pattern of predictions of both methods (using the same predictor variables) is quite different and the AIC is smaller for the multinom approach. I guess polr has more strict premises for the structure of the response variable, which methods can be used to test for

Dear all, I would like to create a stacked area chart to show the development of species biomass over time. Since it is intended for publication I need to prepare in black and white. I have tried to modify the "stackedPlot" function submitted to this list by Christian Lasarczyk on Tue, 16 Aug 2005 using shading lines rather than colors to distinguish between the species.