similar to: Conditional plot length in R

Hello again R users, I have a devilishly hard problem, which should be very simple. I hope someone out there will have the answer to this on the tip of their tongue. Please consider the following toy example: x <- read.table(textConnection("y x1 x2 indv.1 bagels 4 6 indv.2 donuts 5 1 indv.3 donuts 1 10 indv.4 donuts 10 9 indv.5 bagels 0 2 indv.6 bagels 2 9 indv.7 bagels 8 5 indv.8

This is a general statistics question so I'm sorry if its outside the field of r help. Anyway, I have a suite of female and male traits and I have made a matrix of correlation coefficients using rcorr(). This results in a 6 by 6 matrix like this.. [1] 0.11287990 0.20441361 0.23837442 0.04713234 0.04331637 0.01461611 [7] 0.22627981 0.11720108 0.14252307 0.19531625 0.29989953 0.09989502

Dear R-list members, I am calculating the linear extrapolation for a data set, using the function found in Hmisc. x=c(0.0265,-0.0003,0.0142,0.0263,0.0634,0.1145,0.2504) y=c(58,107,152,239,362,512,724) x1=0.0393216 approxExtrap(x,y,x1, method="linear") approx(x,y,x1) #to see what is happening: plot(x, y, typ="o") abline(v=x1, col=8) Which gives x=0.03 and y=163, instead of

Hi, I have a 2D chart that is divided into four quadrants, I, II, III, IV: plot(1:10,ylim=c(0,10),xlim=c(0,10),type="n") abline(v=5,h=5) text(x=c(7.5,7.5,2.5,2.5),y=c(2.5,7.5,7.5,2.5),labels=c("I","II","III","IV")) I would like to fill each quadrant with a background color unique to the quadrant. Does anyone know how to do this in R? Thanks, --

Hello All, Let's say I have points on a x-y plane. x ranges from 0-2 and y from 0-2. There are points in quadrants x[0:1]---y[0:1] and in x[1:2]----y[1:2]. I would like to get the mean and std of the points in the x[0:1]----y[0:1] quadrant alone. Is there a straight forward way to do it? I asked a similar question a few days ago regarding plotting a subset of data using conditions. The

I just upgraded our optimizer to LLVM 3.0 from 2.8 and noticed that the scalar replacement of aggregates pass takes a lot longer for some code. Has there been a performance regression in this pass, or does it do more work? LLVM 3.0: Total Execution Time: 1.0600 seconds (1.0526 wall clock) ---User Time--- --System Time-- --User+System-- ---Wall Time--- --- Name --- 0.5100

Hello, My name is Craig and I need help plotting a "line" for a multiple linear regression in R. Here is my sample data (filename: convis.txt) Output of convis.txt is (vis and density being predictors of either avoidance or entrance): vis den avoid entrance 1 10 1 0.0000 0.0000 2 10 3 0.8750 0.0000 3 8 3 0.8180 0.0300 4 8 3 0.6670 0.0667 5 8 1

Hi, I'm reading part of a table from postgres, so I'm getting a data frame. how can I extract the numerica values so I can operate on them. > res <- dbGetResult(mydata) > str(res) `data.frame': 5 obs. of 8 variables: $ cyx.1: num 0.149 -0.278 0.114 0.060 0.109 $ cyx.2: num 0.158 -0.070 0.063 0.149 0.150 $ cyx.3: num 0.052 -0.350 0.114 0.126 0.238

Hello - I have a Premisys Slimline Channel Bank connected to a Digium TE110P. I am not able to call the FXS extensions or get dialtone on them. The channel bank is connected via a T1 crossover to the cable and lights show green. I really need to get this functioning by end of day. If anyone can help me out I would be greatly appreciative. Thanks, Dave zaptel.conf loadzone =

Hi everyone, I have a PRI conection on an * system running Asterisk 1.0.9, libpri 1.0.9 and zaptel 1.0.9.2 connected to an AXE 10 (APZ 21220 System 64) in the network side. I know the system and the wildcard I´m using are ok because I´ve used them before with other PRI connections (to a Siemens EWSD) without any problem. First the PRI didn´t work (I got the TE110P alarmed with red). After

Warning: I am a complete newbie to R. I have read ISwR, but I am still finding myself completely stuck on some simple concepts. I have tried everything I can think of to solve this one, and finally decided that enough was enough and I need a pointer to a solution. I have the following summary from lm(): ---- > summary(lm(nu1~nu4)) Call: lm(formula = nu1 ~ nu4) Residuals: Min 1Q

Hi my friends, I'm trying to make a scatterplot like this. 1) I have a 3-variable dataset. They are ID, x, and y. 2) "x" is for the X-axis, "y" for the Y-axis, and "ID" is used to label all the cases in the scatterplot. 3) After creating the scatterplot, I need to add both a X-axis reference line and a Y-axis reference line. The X-axis reference line is a

I am new to R.I want to draw grid from a csv file which contains latitude minimum ,latitude maximum ,longitude minimum ,longitude maximum.The grid should be divided into exactly 4 quadrants. The map is of NY state of USA. I want to know how can I do it. Help would be appreciated. Thanks Jaimin [[alternative HTML version deleted]]

Hi everyone, I am Evelina and i have just joined the forum and hope to be able to support soembody as well as finding some help for my stats problem. I am trying to plot Frequencies against Densities of an ecological community using a Olmstead-Tukey diagram (O-T diagram). it is based on the Corner test for Association. I seem to undestand that the test can be run in R using the cor.test()

Dear all, I am running a MCMC on my data, and I want to plot the results of the simulation. I want to dedicate a page in the PS file for each element in my sample, and in each page I need to plot 8 quadrants. In attachment you find my first attemp (just the first page). Here are my troubles (I'm sorry if they are stupid, but I am a newby with R): -) I need to substitute the greek letters

Hi all. This is perhaps more a statistics question but I'm hoping someone can help me. I have a group of patients for whom I'm looking at beat to beat RR interval changes. I have plotted the difference between one beat length and the next against the difference between the previous beat length and the current one. This gives me a plot with four quadrants: the bottom left corresponding to

Hello all: I am confused about the output from a lm() model with an incomplete design/missing level. I have two categorical predictors and a continuous covariate (day) that I am using to model larval mass (l.mass): leaf.species has three levels - map, syc, and oak cond.time has two levels - 30 and 150. There are no response values for Map-150, so that entire, two-way, level is missing.

Hi all. For a while, gcc has supported a stack protection mechanism (-fstack-protector and friends, available in gcc 4.1.2 and up). Can anyone think of a good reason not to enable it if the compiler supports it? A quick test here shows minimal difference in runtime over a full regress pass (~10sec over 8.5 minutes, and since the machine is not entirely idle that could be experimental error).

Hello, While testing recent snapshots (20131023 and 20131024) I encountered a problem with ProxyCommand. The regression tests all passed, but the use of ProxyCommand's in my ~/.ssh/config resulted in name resolution errors; even if CanonicalizeHostname was explicitly set to "no." The patch included inline below fixed the issue: Index: ssh.c

Hello. I'm a new user of R, and I have a question regarding the use of aov and lm-functions. I'm doing a fractional factorial experiment at our production site, and I need to familiarize myself with the analysis before I conduct the experiment. I've been working my way through the examples provided at http://www.itl.nist.gov/div898/handbook/pri/section4/pri472.htm