similar to: kappa-accuracy and test for signifcance

Hi list I'm looking for statisitic measurements describing the pattern of points within a given polygon. Is there a function calculating the Index of Dispersion and/or are there other functions summarising an observed pattern? Markus Schwarz ..................................................................... Markus Schwarz Wissenschaftliche Mitarbeiterin Eidg. Forschungsanstalt WSL

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

Buenas, ¿alguien sabe qué puede ser este error? data es una df y las columnas 23 y 24 son las que quiero comparar; están bien. testK<-kappa2(data[,c(23,24)], "equal") Error in sort.list(y) : 'x' must be atomic for 'sort.list' Have you called 'sort' on a list? -- Dr Manuel Mendoza Department of Biogeography and Global Change National Museum of Natural

Hi, I've got two vectors with ratings from two coders, like this: x<-c("red", "yellow", "blue", "red") #coder number 1 y<-c("red", "blue", "blue", "red") #coder number 2 I want to find Cohen's Kappa using the wkappa function in the psych package. The only example in the docs is using a matrix, which

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 R-user, Did anybody know how to calculate the kappa statistics of two vectors for the agreement? for example: a<-c(0,1,0,0,1,0) b<-c(0,1,1,0,0,1) I know the percent of agreement is 3/6, but how to get the kappa? [[alternative HTML version deleted]]

Hello, I?m looking for a way to measure the goodness of fit of my model with Cohen?s Kappa (scaling between 0 and 1). The kappa function does not give the results I?m looking for. Heres the code: z<-glm(x~y,binomial) kappa(z, exact = T) Does anyone know more? many thanks Christian -- _______________________________________________ Search for businesses by name, location, or phone number.

Greetings: I'm running R-1.9.1 on Fedora Core 2 Linux. I tested a proportional odds logistic regression with MASS's polr and Design's lrm. Parameter estimates between the 2 are consistent, but the standard errors are quite different, and the conclusions from the t and Wald tests are dramatically different. I cranked the "abstol" argument up quite a bit in the polr

Hello! I have question about Debian package management with puppet. I''m wondering is there sane way to make puppet respects packages pinning? i.e., if I have several repos for one package, let''s say it is "nginx" which can be found in lenny & lenny-backports repos. I''ve created pinning file like: Package: nginx Pin: release a=lenny-backports Pin-Priority:

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, I have 50 raters and 180 cases which are rated as malignant, probably malignant, probably benign, and benign. I want to compare all the raters but I want a weighted kappa to penalize differences between malignant and benign more than differences between malignant and probably malignant. I only found the weighted option in the 2-rater functions (kappa2, wkappa, cohen.kappa). 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 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

The irr package seems to install correctly: > install.packages("irr") trying URL 'http://cran.us.r-project.org/src/contrib/irr_0.61.tar.gz' Content type 'application/x-tar' length 13848 bytes opened URL ================================================== downloaded 13Kb * Installing *source* package 'irr' ... ** R ** data ** help >>>