similar to: Simulation

I wish to simulate: Ho: Odds Ratio =1 H1 Odds Ratio <> 1 Can any of the R users show me how to generate this data for say 20% from Ho and 80% from H1? Thanks, Jim [[alternative HTML version deleted]]

Hi, I would like to create a 3 dimensional barplot of 16 odds ratios that demonstrate an interaction between two variables (CD14 and CD23). Here is an example of what the graph I would like to produce (http://www.xfy.com/manual/dev/developer/1.4/spec/image/chart/type_3d.gif): Here is are the odds ratios I would like to graph: Odds ratios CD14-1 CD14-2 CD14-3 CD14-4 CD23-1 1 0.61

Hi! I have a dataset with three columns -the first two refer to x and y coordinates, the last one are odds ratios. I'd like to plot the data with x and y coordinates and the odds ratio shown as a fourfold plot, which I prefer to do using the stars function. Unfortunately the stars option in symbols is not as cool as the stars function on its own, and now i can't figure out how to do it!

Hello everyone, Can someone tell me exactly how to generate data from a null distribution for the fisher exact test? I know I have to use the hypergrometric but exactly what commands do I use? Jim [[alternative HTML version deleted]]

#See David Williams' book "Weighing the odds", p286 y <- c(1.21, 0.51, 0.14, 1.62, -0.8, 0.72, -1.71, 0.84, 0.02, -0.12) ybar <- mean(y) ylength <- length(y) ybarv <- rep(ybar, ylength) x <- 1:ylength plot(x,y,asp=1,xlab="position",ylab="ybar",type="n",ylim=c(-1,1)) segments(x[1], ybar, x[ylength], ybar) segments(x,ybarv,x,y)

Hello, I am planning of building a list of lists specifically, my first list is some what of the sort: lidta <- list(m, p, r, s, q, A, B) where A and B are matrices that may be of different number of rows . The number of rows in matrix A and matrix B depends on the the values of m. The question is I don;t know how to put all the 1000 or so of these lists into a 'mega' list. Can you

Hello, Anyone who knows a fast and accurate algorithm for generating draws from an underdispersed Poisson distribution. Or even better, if there is a package containing such an implementation. Thanks Rene

Rhelpers: I'm trying to make a barchart of a 2-group dataset (barchart(x~y,data=data,groups=z,horizontal=FALSE)). My problem is that I can't, for the life of me, seem to get rid of the inter-bar space -- box.ratio set to 10000 doesn't do much. Any ideas? I'd ideally want zero space between the bars. Thanks! --j -- Jonathan A. Greenberg, PhD Assistant Project Scientist

Hi, I need to produce an ordinary scatter plot and it is vital that the aspect ratio equals 1. I set the axis as: plot(x, y, type="n", asp=1, ,ylim=c(-80,70),xlim=c(0,100)). The problem is that I get some 'additional' blank plot area (basically, the lower bound of xlim becomes quite negative). The xlim is not the range of the x-data, but I need the 0 there for further

How do I make a picture that is a horizontal strip? I tried > plot(x=c(1,2,3,4),y=c(1,1,1,1)) #works but screen image is square. > pdf("ratio.pdf",height=1,width=6) > plot(x=c(1,2,3,4),y=c(1,1,1,1)) I got the following error message: Error in plot.new() : figure margins too large Is it possible to produce an on-screen picture that is a horizontal strip? (I use Mac Os X

Hi, How to mark points on x axis of a graph keeping x axis as constant and changing y from y1 to y2 respectively. I want to highlight the area from y1 to y2. Any suggestions Thank you Jeet [[alternative HTML version deleted]]

Hi all, I want to calculate in each row a ratio based on number in the current row and the previous row. Is there a way to do this without for-loops because that is extremely slow. A B [1] 2 2 [2] 2 3 [3] 4 5,5 ... B2 = A2 + 0.5*B1 Thanks in advance. Best Marcus

The tkrplot package is not working in version 2.13.0 for windows. I contacted the maintainer who unfortunately does not have easy access to a windows computer and says that it is working on the other platforms. I traced the problem down to the line in the .First.lib function: .Tcl(paste("load", file, "Rplot")) With file being "C:/Program

hi, Excuse me asking three questions in a row for a day, but I had collected those questions as I'm still experimenting with R. This one is how do you compose plots with alot of data in one graph. First what I currently do is after i generated all the data I need to plot.. let say 5-6 arrays I do a plot() and then I created a function to scale the other data to approx fit the initial range so

Hi, I have problem with simulating. This is my task... Suppose that there are N persons some of whom are sick with influenza. The following assumptions are made: * when a sick person meets a healthy one, the chance is á that the latter will be infected * all encounters are between two persons Write a function which simulates this model for various values of N (say, 10 000) and á (say, between

Hello there, I have a fairly simple request (I hope!) I have produced a filled contour plot like this: library(grDevices) library(gplots) library(plotrix) filled.contour(contour, axes=F, frame.plot=TRUE, color=terrain.colors, ylab= "Length Along Flume (m)", key.title = title(main="Velocity\n(m/s)"), key.axes = axis(4, seq(0, 0.6, by = 0.1)), asp=2, plot.axes = {

# how would I code in R to look at the letter of the alphabet # in the second column and create a indicator column for the # corresponding letter? data(InsectSprays) InsectSprays$spray

Oh sorry. An example say lots more than words. The data below, when submitted to twoord.plot return the mentioned error. Rain are bars and Salt lines, the bars appear and the error occurs with the salt data. Point Rain Salt Fev/03 365.6 13 Mar/03 235 18 Abr/03 115.1 18 Mai/03 47.4 18.75 Jun/03 112 15 Jul/03 156.8 17 Ago/03 66.1 15 Set/03 149.8 14 Out/03 167,1 11.5 Nov/03 269.3 17.5 Dez/03 283.7

Hello all!!! 1) I am wondering is there a way to generate random numbers in R for Bivariate Uniform distribution? 2) Does R have  built-in function for generating random numbers for any given bivariate distribution. Any help would be greatly appreciated !! Good day! Haneef Anver [[alternative HTML version deleted]]

hello, i did several binomial tests and noticed for one sparse dataset that binom.test(1,1,0.5) gives a p-value of 1 for the null, what i can't quite grasp. that would say that the a prob of 1/2 has p-value of 0 ?? - i must be wrong but can't figure out the right interpretation.. best, kay ----- ------------------------ Kay Cichini Postgraduate student Institute of Botany Univ. of