similar to: Contingency tables example

Hi, I guess that my problem has an obvious answer, but I have not been able to find it. Suppose I create a custom function, consisting of two beta-distributions: myfunction <- function(x) { dbeta(x,2,6) + dbeta(x,6,2) } How can I calculate the quantiles of myfunction? I have not seen any continous function treated in the docs, and applying the "quantile function" gives me an

Hello dear R-helpers, I'm working with R-2.15.2 on Windows 7 OS. I'm stucked with a merge of two data frames by characters. In each data frame I got two different list of names, that is my main-key to be merged. To figure out what I'm saying, I build up a modified "?merge" example, with errors by purpose: # Data for authors: authors <- data.frame( surname =

Hi, I am using the reldist package and having problems determining the p-value for the entropy value from the reldist function. I am able to properly determine the entropy value, but cannot figure out what function to use to find the p-value. I have tried using rpy, rpluy (which provides p-values for the polarization values) and investing the results from reldist(). Thus, far I cannot find the

Hi there! I´m a noob when it comes to R and I´m using it to run statisc analysis. With the code for ARIMA below I´m getting this error: Error in solve.default(res$hessian * n.used) :Lapack routine dgesv: system is exactly singular The code is: > s.ts <- ts(x[,7], start = 2004, fre=12) > get.best.arima <- function (x.ts, maxord=c(1,1,1,1,1,1)) + { + best.aic <- 1e8 + n <-

Command I am using for the plot: plot(Raw[][Plate==101]~well[][Plate==101], xlab="Well", ylab="Raw", main="Plate 101") I only want to label points on the graph where well equals B13, G13, K14 and N14 with the name of the well. Thank you for your help. -- View this message in context:

Dear R-helpers, I'm not a speciallist in writing complex functions, and the function still very rusty (any kind of suggestions are very welcome). I want to implement Gale-Shapley algorithm for R Language. It is based on http://www.jstor.org/stable/10.2307/2312726 Gale and Shapley (1962) , and it has evolved to

I've got a 20 slide Powerpoint set up against the expected day someone says "Mitch, could you talk for a short while to the residents on how to use statistics?" It looks good at this end, but I'd love to get at least one more pair of eyes on it. Any volunteers? -- Due to the recession, requests for instant gratification will be deferred until arrears in scheduled

Hello List, I have a client (a nursing home) that we are looking at installing a trixbox for. One of the features that they would really like is a customized, standard voicemail recording for each of the residents rooms. We are looking for something along the lines of a voicemail recording like this: "Welcome to (nursing home). You have reached room 5. Please leave a message after the

In the subdivision where I live, we have a well that time to time has problems. Currently, our trustees call me, I take a message, and then call the people on our "phone calling tree" and give them the message. They, in turn, pass the message on to the rest of the residents. We have a few problems with the way this works -- people don't bother to relay the message, or,

Hello fellow downstream residents, I see that r358546 accidentally deleted an entire subtree, which was reverted in r358552. This of course caused a big merge conflict in our local repo, and internally we've been debating tactics for dealing with it, hopefully without losing our original history. Has anyone else handled this in a way that they are happy with? We found a StackOverflow

Hello, We're using Samba 4.1.11 as domain controllers and over the past two weeks I've run into several issues with unrelated Windows software, the problems of which all point to Kerberos authentication and SPNs as being somehow involved. If there are many more issues it might start to get politically difficult *not* to blame the DCs, and I don't want to point fingers at Samba. Are

Hi all, Looks like I'm going to make the trip over to Astricon next month, but finances being what they are and since I'll be paying for my own flight from the UK, I'm trying to cut down on costs. The problem I've got is that the hotel is about 4 miles from the nearest public transportation (Brookhaven station) and the shuttle bus is only for hotel residents. Hiring a car seems

What is the status of i7 architecture support for CentOS-5? Do the latest updates support it? -- This e-mail may contain technical information which is controlled by the United States Government, Department of State, International Traffic & Arms Regulations (ITAR) (22 CFR 120-130) which requires an export license prior to sharing with foreign persons. Lacking such license, ITAR technical

Dear Listers: I have a categorical feature selection problem for random forest. Suppose I have a multiple-leveled category variable A, which has m=3 levels: red, green, and blue and the final target is binary classification. I want to evaluate its power in discrimination between 2 classes. We know rf splits multiple-leveled category variable by considering all combinations of its levels. So

v 2.0.1 (sooooh old!) on Win2k I think I know the answer to this but I can hope ... I have data for continuous variables (measures of residents) by a categorical variable in range (1,22), the units in which they live. I want to plot these data with a pair of boxplots one above another with same x-axis (1,22) using par(mfrow=c(2,1)) and then plotting first for the women then for the men.

I have a contingency table with a 0/1 variable (fwr) that flags a particular condition, a weight variable (cnt) and other variables (e.g. zz) characterizing that observation of the table. I am trying to use GLM. By converting the variables to factors and running the following regression: fwr1<-factor(fwr) zz1<-factor(zz1) res1<-glm(fwr1 ~ zz1,weights=cnt) Is that appropriate? Paul M.

Hello All, Is there any R function out there to turn a multi-way contingency table back to a flat file table of individual rows and attribute columns.? Thanks! marco --------------------------------- [[alternative HTML version deleted]]

Hi all, I have a 3x4 contingency table with row totals all being 100. I want to generate 3 x 4 tables from the null distribution. Which R function can do this? -- Thanks, Jim. [[alternative HTML version deleted]]

Dear R-philes: I am having an issue with exporting contingency tables with xtable(). I set up a contingency and convert it to a matrix for passing to xtable() as shown below. v.cont.table <- table(v_lda$class, grps, dnn=c("predicted", "observed")) v.cont.mat <- as.matrix(v.cont.table) Both produce output as follows: observed predicted uh uh~ uh 201

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]]