similar to: Add data count to plot?

I created a 2 D plan: > plot(range(-2.5,0.95),range(0.00,1.00),type="n",axes=TRUE) I made a projection of points with their coordonates (X,Y) in that plan > fullpoints=read.csv2("fullpoints.csv",h=T) > plot(fullpoints) The points are listed in that .csv file it is organized this way: 0,48875 0,142857143 0,409 0,142857143 0,45611 0,25 0,49833 0,222222222 0,61158

Hello dear R-help group. My task looks simple, but I can't seem to find a "smart" (e.g: non loop) solution to it. Task: I wish to randomize a data.frame by one column, while keeping the inner-order in the second column as is. So for example, let's say I have the following data.frame: xx <-data.frame(a= c(1,2,2,3,3,3,4,4,4,4) , b =

Hello dear R members. I have been learning the Anova syntax in order to perform an SS type III Anova with repeated measures designs (thank you Prof. John Fox!) And another question came up: where/what are the (between/within) residuals for my model? ############ Play code: phase <- factor(rep(c("pretest", "posttest", "followup"), c(5, 5, 5)),

Hi everybody, I would like to count how many times names in list L, nombreL, apear in list C, nombreC. Can I improve the next program? cuenta <- 0 topL <- length(nombreL) topC <- length(nombreC) for (i in 1:topL) { for (j in 1:topC) { k <- grep(noquote(nombreL[i]),nombreC[j])

Hello R help members. I came across today with an article on Barnard's exact test (http://www.cytel.com/Papers/twobinomials.pdf), that is supposed to give a more powerful fisher.test - Because it doesn't assume that we know the row and column totals are in advance. Any pointers to such a function ? Thanks, Tal -- ---------------------------------------------- My contact information:

I would like to know if you have any suggestions how to visualize the results from a paired t-test (see the example data below). I tried to produce plots that show the mean and CI's from the original data and the estimate of the difference between means and the confidence intervals (see below) from the t-test. I really don't know what would be the best way to graphically display the

Hi All, I am new comer here and for Statistics R. I want to know how to use R to calculate one way repeated measure ANOVA, for example, on three groups of data. I appreciate if someone can help me solve this problem. John

*(note*: This is an R community question, not a statistical nor coding question. Since this is my first time writing such a post, I hope no one will take offence of it.) Hello all, I will be attending useR 2009 next week, and was wondering if there are any of you who are *bloggers *intending to participate and report on useR 2009? If so - I would love to know your blogs URL so as to follow you.

Hello Martin Maechler and All, A simple question (I hope): How can I compute the "sum of the dissimilarities" that appears in the pam command (from the cluster package) ? Is it the "manhattan" distance (such as the one implemented by "dist") ? I am asking since I am running clustering on a dataset. I found 7 medoids with the pam command, and from it I have the

(I use R 2.9.1 with win XP) If I run this code: seq(-0.1,.9, by = .05)[seq(-0.1,.9, by = .05) <= 0.5] I get this output: [1] -0.10 -0.05 0.00 0.05 0.10 0.15 0.20 0.25 0.30 0.35 0.40 0.45 Why is 0.50 not in the results ? (It seems that it gives a slightly bigger number then 0.5 but I don't understand why it does that) Where as if I try: seq(-0.1,.9, by = .05)[seq(-0.1,.9, by = .05) <=

Hello all. I wish to run k-means with "manhattan" distance. Since this is not supported by the function "kmeans", I turned to the "pam" function in the "fpc" package. Yet, when I tried to have the algorithm run with different starting points, I found that pam ignores and keep on starting the algorithm from the same starting-points (medoids). For my

Hello dear R mailing list members. I have recently became curious of the possibility applying model selection algorithms (even as simple as AIC) to regressions of large datasets. I searched as best as I could, but couldn't find any reference or wrapper for using step or stepAIC to packages such as biglm. Any ideas or directions of how to implement such a concept ? Best, Tal --

Hello people. I wish to reorder a simple vector of numbers by another vector of the order (and then do the same, but with a data frame rows) I try this (which doesn't work) : > aa <- c(3, 1 ,2 ) > aa[aa] [1] 2 3 1 The same won't work if I try to order a data frame: > data.frame(matrix(c(3,1,2), 3,2))[c(3,1,2),] X1 X2 3 2 2 1 3 3 2 1 1 How should I do this? p.s:

(Have searched for this already) Hi, How do you find the strength of correlation between two variables using an ANOVA table? "Pr(>F)" gives the statistical significance of the association, but not the strength of the correlation. See data (from R) below Readable: "Df" "Sum Sq" "Mean Sq" "F value"

Hello dear R-help group (and David Smith from REvolution), I would like to perform parallel computing using R with Condor (hopefully using foreach or other recommended solutions, if available) for some "Embarrassingly parallel" problem. I will start by listing what I found so far, and then go on asking for help. So far I found the a manual by Xianhong Xie from Rnews_2005-2 (see page

Hi, I'm trying to use the function "error.bars", but the program don't find it, and I dont't found any package with this function. Is there some another functin to draw barplots with error bars? Sueli Rodrigues Eng. Agr?noma - UNESP Mestranda - USP/ESALQ PPG-Solos e Nutri??o de Plantas Fones (19)93442981 (19)33719762

Hello all, The following line of code includes a right-to-left language text, yet the R graphics engine displays it from left to right. One problem this causes is when there are parenthesis in the test, here is a basic example? plot(1:10, main = "שלום (טקסט)") Is there a way to make sure the text is displayed from right to left? Many thanks for any suggestions, Tal

Hello dear R help group. My question is statistical and not R specific, yet I hope some of you might be willing to help. *Experiment settings*: We have a list of subjects. each of them went through two tests with the answer to each can be either 0 or 1. *Goal*: We want to know if the two experiments yielded different results in the subjects answers. *Statistical test (and why it won't

Hello all. I wish to read a large data set into R. My current issue is in getting the data so that R would be able to access it. Using read.table won't work since the data is over 1GB in size (and I am using windows XP), so my plan was to read the file chunk by chunk and each time move it into bigmemory (I'll play with that when the time will come, maybe ff is better ?!). I encountered

I have the following codes but can anyone make it shorter i.e making these FOR loop into one loop ... thanks... par(mfrow=c(2,4)) for(i in 16:23){ hist(data[,i],main=paste(colnames(data)[i],sep=""),ylab="Frequency",xlim=c(1,5),xlab="Score",ylim=c(0,100)) } png("histogram.png") dev.off() par(mfrow=c(2,4)) for(i in