similar to: Fisher Test in R

Hi Everyone, I ran the same code in R and in R-studio, but got two different results. Does anybody know why this is occurring, and if there is a fix for this? and which is the "correct" program to use ? Some information about the code I am running: I am running the fisher test and it seems that the p-values are similar but not same for example, for an event A the p-value coming from

I would like to fit a hierarchical regression model from Witte et al. (1994; see reference below). It's a logistic regression of a health outcome on quntities of food intake; the linear predictor has the form, X*beta + W*gamma, where X is a matrix of consumption of 82 foods (i.e., the rows of X represent people in the study, the columns represent different foods, and X_ij is the amount of

Hi All: I am working with a dataset on Arsenic toxicity, and I am trying to calculate the 20th, 40th, 60th, 80th, and highest percentiles for a variable, dietary Moisture (variable name dMoist). The inbuilt function quantile(dMoist) would print 0, 25th, 50th, 75th, and 100th percentile. Does there exist a function that can calculate xth percentile (where x = 10th, 20th, ... etc) values?

Dear All, I wonder if someone can point me in the right direction here. I'm working with the rms library, R 2.9.2 under Windows XP. I'm trying to arrange two plots side by side for a colleague. mfrow or mfcol do not seem to work, however, so I am obviously missing something important. I know that there have been changes in the graphics from Design to rms, but am just not sure where to

Dear LLVM Users, The European LLVM User Group Meeting 2011 schedule is on-line at http://llvm.org/devmtg/2011-09-16/. Thank you every one that registered and submitted papers and posters, There are 5 talks planned, several side-sessions (one of them reserved for lightning talks) and more than 40 people already registered! The current proposed list of side sessions are: - OpenCL (and similar

What is the downside of sticking with 5.2

Hi R experts, I am interested on the effects of two dietry compunds on the growth of chicks. Rather than extracting linear growth functions for each chick and using these in an analysis I thought using ReML might provide a neater and better way of doing this. (I have read the pdf vignette("MlmSoftRev") and "Fitting linear mixed models in R" by Douglas Bates but I am not

fletcher said: > For any consensus to come about, > I think we need to agree on the > fundamental purpose and philosophy of > the consensus we claim to be interested in. it would be nice. > Otherwise many of these discussions will > continue to occur without much hope of > moving forward to any actual outcome/resolution. yep. > it's

* Reduce the amount of sleep you need * Cause wounds to heal faster * Lose weight while your sleeping * Become less winded when excersizing * Put color back in grey hair * Grow hair back where it had once fallen out * Tighten skin * Strengthen bones * Body builders - use this to build your muscles quicker ..........The List truly goes on and on.......... As seen on NBC, CBS, CNN, and Oprah! The

So, it's advisable to grow cane as near to the factory as possible. He is the fourth person to resign at NFA. They sold it to Russia at virtually give away prices. uuuu WATCH LITL ON WEDNESDAY DEC 13! Company: L INTL COMPUTERS INC (Other OTC:LITL.PK) Symbol: LITL Mon Close: $0.99 5-day Target: $5 L Reveals 8-Way Processing Professional Mobile Computing Design CALL YOUR BROKER TO PLACE

Hi all, I wanted to have a seasonality study like whether a particular month has significant effect as compared to others. Here is my data : 0.10499 0 0 1 0 0 0 0 0 0 0 0 0.00259 0 0 0 1 0 0 0 0 0 0 0 -0.06015 0 0 0 0 1 0 0 0 0 0 0 0.10721 0 0 0 0 0 1 0 0 0 0 0 0.03597 0 0 0 0 0 0 1 0 0 0 0 0.10584 0 0 0 0 0 0 0 1 0 0 0 0.02063 0 0 0 0 0 0 0 0 1 0 0 -0.03509 0 0 0 0 0 0 0 0 0 1 0 -0.03485 0 0 0

Hi, I'm having difficulty with getting a table to show with multiple rows and columns. Below is the commands that I've typed in and errors that I am getting. Thank you. Laura Table trying to enter: Diet: Binger-yes: Binger-No: Total: None 24 134 158 Healthy 9 52 61 Unhealthy 23 72 95 Dangerous 12 15 27 >

Otras variantes con y sin paquetes adicionales... > sapply(split(datIn$Gain, as.factor(datIn$Diet)), mean) d1 d2 d3 280 278 312 > by(datIn$Gain, datIn$Diet, mean) datIn$Diet: d1 [1] 280 -------------------------------------------------------------- datIn$Diet: d2 [1] 278 -------------------------------------------------------------- datIn$Diet: d3 [1] 312 > > library(dplyr) >

Me gusta la respuesta uqe has dado, pero si por ejemplo, alguno de los datos tiene datos faltantes, entonces devuelve NA. He probado con: sapply(split(datos$uno, as.factor(datos$dos)), mean(na.rm=TRUE)) pero da fallo. ¿Cómo se podría hacer para que devolviera además la media obviando los NA y que contara el numero de NA por categoria? > Date: Wed, 28 Oct 2015 00:13:45 +0100 > From: cof

Estimados Cuando existia epicalc, había una manera muy fácil de determinar la media de una variable (en esta caso Gain) por grupos, en este caso (Diet). ?Como se puede hacer ahora? Diet Gain 1 d1 270 2 d1 300 3 d1 280 4 d1 280 5 d1 270 6 d2 290 7 d2 250 8 d2 280 9 d2 290 10 d2 280 11 d3 290 12 d3 340 13 d3 330 14 d3 300 15 d3 300

Dear all, This question may be too basic quesition for this list, but if someone has time to answer I will be happy. I have tried to find out, but haven't found a consice answer. As an example I use "Pinheiro, J. C. & Bates, D. M. 2000. Mixed-effects models in S and S-PLUS. Springer, New York." page 225, where rats are fed by 3 different diets over time, which body mass has

Hello This is Elaine. I am using package lattice to generate boxplots. Using Richard's code, the display was almost perfect except the outlier shape. Based on the following code, the outliers are vertical lines. However, I want the outliers to be empty circles. Please kindly help how to modify the code to change the outlier shapes. Thank you. code package (lattice) dataN <-

Hi, In the examples from the ReShape package there is a simple example of using melt followed by cast that produces a smallish amount of output about the chicks database. Here's the code: library(reshape) names(ChickWeight) <- tolower(names(ChickWeight)) chick_m <- melt(ChickWeight, id=2:4, na.rm=TRUE) DietResults <- cast(chick_m, diet + chick ~ time) DietResults My challenge

Hi all , I have a dataframe (Hypertension) with following headers :- > Hypertension ID Hypertension(before drug A) Hypertension(On drug A) On drug B? Healthy diet? 1 160 90 True True 2 190

I purchased a Fisher Price Computer School program. It's a kid friendly USB keyboard (can be viewed here http://www.fisher-price.com/fp.aspx?st=10&e=ccslanding) I'm running Ubuntu 10.4 Lucid which also has the latest version of Wine. After installing I get this error Exception EExternalException in module ntdll.dll at 0002E3F8. External exception C0000025 Any ideas on this error. It