similar to: Analisys in Multidimensional contingency tables

How does one ideally handle and display multidimenstional contingency tables in R v. 2.6.2? E.g.: > prob1<- data.frame(victim=c(rep('white',4),rep('black',4)), + perp=c(rep('white',2),rep('black',2),rep('white',2),rep('black',2)), + death=rep(c('yes','no'),4), count=c(19,132,11,52,0,9,6,97)) > prob1 victim perp

Dear R-list, I know is going to be weird, but I need to do something about and I know the R-list is big group of good people. I?m moving to Philadelphia in August 2009 (my fianc? is there) and I need to find a job. I know to do so many things, but my skill is in science and research. I would like to know if you know a possible position near Philadelphia, New Jersey or New York. I was working in

Hello, I am searching for an R-Package that does an exentsion of the Mantel-Haenszel test for ordinal data as described in Liu and Agresti (1996) "A Mantel-Haenszel type inference for cummulative odds ratios". in Biometrics. I see packages such as Epi that perform it for binary data and derives a varaince for it using the Robbins and Breslow variance method. As well as another pacakge

Dear R help gurus, I have the following problem and I would be delighted if you could help me. >From a large (1500) cohort of patients we have been taking some measurements (ECG measurements, but its not important). The measurements are ordinal in 4 grades (Grade I-IV, grade IV being the most severe form). Every patients has been measured several times (usually once per year). The

Suppose that we have o 2-D contingency table where the row variable is nominal and the column one is ordinal. In SAS it is possible to compute the statistic named as row mean scores differ. How can we programmed it in R? (See also Aggresti (2002), Categorical Data Analysis, p. 302) With regards --------------------------------- [[alternative HTML version deleted]]

dear colleagues, I want to use hierarchical partitioning methods to evaluate variance contributions of differente predictors within a generalized linear mixed model framework. I don't know if this is possible, because I can not found anything about this issue. My error is poisson and I, at least, a random factor in the designed analysis. I acknowledge any idea about how deal with

Hello, Is there anywhere I can find some Bayesian Code for 2 by 2 tables or even the non-central hypergeometric distribution in R? Packages would be helpful but the actual coding in R is much better. Thanks, Thanks Jim [[alternative HTML version deleted]]

Hi everyone. I have these data: probClass<-seq(0,0.9,0.1) prob1<-c(0.0070,0.0911,0.1973,0.2949,0.3936,0.5030,0.5985,0.6869,0.7820,0.8822) prob2<-c(0.0066,0.0791,0.2358,0.3478,0.3714,0.3860,0.6667,0.6400,0.7000,1.0000) # which I'm plotting as follows: plot(probClass,prob1,xlim=c(0,1),ylim=c(0,1),xaxs='i',yaxs='i',type="n") lines(probClass,prob1)

Dear all, I would like to create a code for semiparametric Klein and Spady's estimator. For that I created a function that provides the log-likelihood function for each observation (so it is a function of betas and i, where i denotes the observation). Now, in order to maximize the log-likelihood function, I have to sum these log-likelihood functions for each i and so to get another function

Hi, I'm new to R and I've run into a problem that I'm not really sure how to express properly in the language. I've got a data table that I've read from a file containing some simple information about the performance of 4 algorithms. The columns are the name of the algorithm, the problem instance and the resulting score on that problem (if it wasn't solved I mark that

Hello everybody :-) I have some data that I want to model with a logistic regression, most of the independent variables are numeric and the only dependent is categorical, I was thinking that I could apply a logistic regression using glm but I wanted to deepen my knowledge of this so I tried to do some reading and found the "iris" dataset, now I would like to ask two things, first

I'm looking for data to use with nonparametrics exploration technics. I'm a undergraduate student in statistics at the Unicamp (Brazil), and i've to make a project (with orientation of a professor). I choose nonparametrics statistics to make the project. But in this project i must have to use a databank and present analisys. So, if somebody give me indications where i can find

Dear all, Is there any spectral analisys functionality available for R version 3.4.3? Series() functionality doesn't work in this version. Regards Danniel [[alternative HTML version deleted]]

Hello, In a previous email, Matthias mentioned that register pressure estimates in the machine scheduler are not absolute; they only account for the registers that are used in the block.I assume that he meant that registers that are live-through (both live-in and live-out) are not accounted for in register pressure calculations. If a register is either live-in or live-out but not both, it must be

Hi, i'm an italian student that use for the first time the mailing list. I need an help to use a function for elaborate a multinomial logit analisys. I'm making a paper on the USA's commons and their way of organize their service. The various tipologies of service are 4: public, contracting in, contracting out and no-profit. I have 10 different service for each common, and for each

Dear list members I want to analyze separately the months of a time series. In other words, I want to plot and fit models for each month separately. Taking the example of http://a-little-book-of-r-for-time-series.readthedocs.org/en/latest/src/timeseries.html births <- scan("http://robjhyndman.com/tsdldata/data/nybirths.dat") birthstimeseries <- ts(births, frequency=12,

Hello I have a rather simple code and for some reason it produces an error message. If someone can tell me why and how to fix it, I would be very greatful. Thank you in advance. ##### create data set.seed(10) n <- 200 # number of training points n.test <- 200 # number of test points p<-2 # dimension of input space z <-

From: Wade Maxfield <maxfield@ctelcom.net> > Many thanks to the moderator who pointed out errors and suggested the > correct information on this post. Over half the ideas are due to him. ;) > 1. Programs put data in local variables in functions. These variables > are on the computer stack. Feeding data to those variables (usually > string variables) causes the

Sorry, I pasted the wrong file in earlier... this is the correct one: pValCalculator(b, n=20, m=20) { ind <- 1:min(c(n,m)) prob <- (1-pnorm(b,sd=std*sqrt(ind))) prob1 <- sum((n-ind+1)*(m-ind+1)*prob) prob1 } inputData <- scan("/users/lvssso/projects/LAMA/output/pValLamaScore.tmp", list(block1 = "",block2 = "",width1 = 0,width2 = 0,alignment = 0,score

{again a message that was sent to owner-r-help (which is me, currently) why on earth ???!??!? reply to R-help or the original sender Brett Presnell; } I'm teaching a categorical data analysis course this term, and a minor "problem" has resurfaced that I have often thought about before. This applies equally to Splus I suppose, but my undergrads aren't using Splus. It