similar to: Wrinting integers in a matrix faile

Hello everyone, I am using the following program to get the p-value of some numbers (column 'LR' of the data.dat file). I want to write the 1st and 2nd column of the output file (data.out) as an integer while the program change them. Could anybody please tell me how I can write the code which writes the values of the first two columns as integer? For your convenience, I have attached the

Hello all, I am doing a Likelihood Ratio (LR) test in my simulation and I have a vector LR values (each with 1 degree of freedom) at the end of my simulation. Can anybody tell me how I can write a 'R' code which gives me the p-value for each of those LR values. Thanks [[alternative HTML version deleted]]

Hello, I have used the following data to draw my barplot: BL LR Q 36.35 1.00 1.92 36.91 4.00 0.00 25.70 6.00 0.00 34.38 3.00 1.92 05.32 0.50 0.00 BL<-c(36.35, 36.91, 25.70, 34.38, 05.32) LR<-c(1.00, 4.00, 6.00, 3.00, 0.50) Q<-<(1.92, 0.00, 0.00, 1.92, 0.00) barplot(dt$LR, main='LR Value', col='orange', border='black',

Hello all, I have this simple barplot code: ifn <- "id.dat" dat <- read.table(ifn) ofn <- "id.png" bitmap(ofn, type = "png256", width = 30, height = 30, pointsize = 30, bg = "white",res=50) par(mar=c(5, 5, 3, 2),lwd=5) par(cex.main=1.6,cex.lab=1.6,cex.axis=1.6) names(dat)<-c("NumberOfPeople","Average")

Hello All: If I do have: x = (2, 4, 5, 5, 6, 4, 5, 2, 1) y = (9, 11.5, 12.5, 13, 14, 19, 20, 21, 22) I wanted to find a simple function in R which calculates the averages of X and Y for every 4 unit increase of Y. The results should look like: x = (4, 6, 3) #where (2+ 4+ 5+ 5)/ 4 = 4 and ... y = (12, 14, 20.5) #where (9+11.5+12.5+13)/ 4 = 12 and ... Thanks, Mohsen

Dear all, I want to design a histogram and I need to have the frequency at certain points. For example I have the following 2 columns: *X Y* 0.1 25 0.4 22 0.45 11 0.55 21 I want the chart to have 4 columns. First column is from 0.0-0.1 (on X) and frequency is 25. Next colum is wider and form 0.1-0.4 with 22 frequency. Next column is narrow with 11 frequency and the last column is

Hello All: I am having the following data file named "data.dat" Number of people Number of pets Age of trees ifn <- "data.dat" dat <- read.table(ifn) colnames(dat)<-c("Variables") I want R to read all these in a string but when I ask R to read these, it gives me error because there more than one word in each line. Any comments to solve this problem?

Hello All: Can anybody tell me what is the problem with my program please. I have an error message as appears below. My program is: ifn <- "Jul08_09.LM" data <- read.table(ifn) ofn <- "Jul.png" bitmap(ofn, type = "png256", width = 30, height = 30, pointsize = 30, bg = "white",res=50) par(mar=c(5, 5, 3, 2),lwd=5)

Hello All: I am having the following error message when I increase the resolution ("res=30" to "res=350") of my graph: Error in plot.xy(xy, type, ...) : ignoring SIGPIPE signal Calls: plot -> plot.default -> plot.xy Execution halted Here is the code: ofn <- "MAY11.png" bitmap(ofn, type = "png256", width = 21, height = 30, pointsize = 30, bg =

Hello All: I was wondering why many plot options are not working with mhtplot. I could not find much info about this function on the web too. Any comments, any detailed manual ... Thanks, MJK [[alternative HTML version deleted]]

Hello All, I am drawing a graph having 18 small graphs inside using par(mfrow = c(6,3)) command. My problem is how to specify the margins of the whole 18 graphs. I used par(mar=c(6.5, 6.5, 1.5, 1.5)) for each graph separately already but it does not left any margins for the 'mtext()' for the margins of the whole 18 graphs. Any comments please. Thanks, MJK [[alternative HTML version

Hi, Does anyone know how to do a coord_trans() in which the y-axis is tranformed into (for example) -1000/y? Thanks, _____________________________ Professor Michael Kubovy University of Virginia Department of Psychology USPS: P.O.Box 400400 Charlottesville, VA 22904-4400 Parcels: Room 102 Gilmer Hall McCormick Road Charlottesville, VA 22903 Office: B011

Hello All: Could anybody tell me please what is the meaning of this error message: ESP Ghostscript 815.02: Unrecoverable error, exit code 1 Error in plot.xy(xy, type, ...) : ignoring SIGPIPE signal Calls: plot -> plot.default -> plot.xy Execution halted My code has to go to 18 loops for three times and in each loop it has to read a file which is 100-300 Megs. It always crashes at loop 14

Hello there, I have the following info in MyData: 01 0.40 02 0.40 03 0.40 04 0.35 05 0.34 06 0.33 I want to read them using the following code: ifn <- "MyData" dat <- read.table(ifn) names(dat1)<-c("Code","M") names(dat2)<-c("Code","M") I want to have the first three lines in

Dear R-helpers, How do I transform v <- c('insd-otsd', 'sppr-unsp') into c('insd--otsd', 'sppr--unsp') ? _____________________________ Professor Michael Kubovy University of Virginia Department of Psychology USPS: P.O.Box 400400 Charlottesville, VA 22904-4400 Parcels: Room 102 Gilmer Hall McCormick Road Charlottesville, VA 22903

How does one define PDFs as yet undefined in R, such as the ex- gaussian, the sum of two RVs, one exponential, one Gaussian? The PDF would then be the convolution of an exponential PDF, dexp(), and a normal, dnorm(). Kindly cc me in your reply to r-help. Thanks, _____________________________ Professor Michael Kubovy University of Virginia Department of Psychology USPS: P.O.Box 400400

Hi, I would like to use the effect() function (actually a slightly modified version of it) on the output of the lmer() function in the lme4 package. But the effects package requires the nlme pacvkage, which is incompatible with lme4. Workaround? ______________________________________________ Professor Michael Kubovy University of Virginia Department of Psychology for mail add: for FedEx or

df1 <- data.frame(a = LETTERS[1:2], b = LETTERS[3:4], c = 1:2) I am looking for an idiom that swaps the elements of df$a and df$b when (e.g.) df$c == 2, resulting in df2 <- data.frame(a = LETTERS[c(1, 4)], b = LETTERS[c(3, 2)], c = 1:2) _____________________________ Professor Michael Kubovy University of Virginia Department of Psychology USPS: P.O.Box 400400 Charlottesville, VA

Dear r-helpers, Does anyone have a straightforward example of putting together three unrelated (expect for a common y-axis) xyplot() figures in what would be in base graphics a par(mfrow = c(1, 3)) arrangement? _____________________________ Professor Michael Kubovy University of Virginia Department of Psychology USPS: P.O.Box 400400 Charlottesville, VA 22904-4400 Parcels: Room 102

library(glmmADMB) #Example for glmm.admb data(epil2) glmm.admb(y~Base*trt+Age +Visit,random=~Visit,group="subject",data=epil2,family="nbinom") Gives: Error in glmm.admb(y ~ Base * trt + Age + Visit, random = ~Visit, group = "subject", : The function maximizer failed ****************** R version 2.4.1 RC (2006-12-14 r40181) powerpc-apple-darwin8.8.0 locale: C