similar to: "with" and "by" and NA:

Hi folks. If I have the following in my "data" event pH1 pH2 1 4.0 6.0 2 4.3 5.9 3 4.1 6.1 4 4.0 5.9 and on and on..... for about 400 events Is there a way I can get R to plot event vs. pH1 and event vs. pH2 and then do a loess or lowess line for each?? Thanks in advance David [[alternative HTML version deleted]]

Dear r-helpers, I have a question about multiple testing. Here an example that puzzles me: All matrixes and contrast vectors are presented in treatment contrasts. 1. example: library(multcomp) n<-60; sigma<-20 # n = sample size per group # sigma standard deviation of the residuals cov1 <- matrix(c(3/4,-1/2,-1/2,-1/2,1,0,-1/2,0,1), nrow = 3, ncol=3, byrow=TRUE, dimnames =

I'm trying to make a set of contour plots of bivariate kernel density estimates, showing three such plots overlaid, similar to this plot http://euclid.psych.yorku.ca/SCS/Private/Test/ridge-boot2.pdf except that I would like to have the contours *filled* (using transparent colors). To make this reproducible, I've saved the results of KernSmooth::bkde2D() in the following file:

Hello, In the following I tried 3 versions of an example in R help. Only the two first predict command work. After : library(splines) require(stats) 1) fm1 <- lm(weight ~ bs(height, df = 5), data = women) ht1 <- seq(57, 73, len = 200) ph1 <- predict(fm1, data.frame(height=ht1)) # OK plot(women, xlab = "Height (in)", ylab = "Weight (lb)") lines(ht1, ph1) 2)

I have a recurring graphics issue that I've not been able to resolve with R. If I make a series of regression estimates and then plot the estimated function for the regression lines over a scatter plot of the data, e.g., using a sequence of plot( ) and lines ( ) similar to those below

Hi Everyone, I am uncertain that I am writing the contrast statements correctly. Basically, I'm unsure when to use a -1 and a 1 when writing the contrasts. Specifically I am interested in comparing the slopes between different temperature regimes. Temperature is therefore a factor. Time and percent are numerical. Using the gmodels package I made the following model:

Hi, I've just installed libao-1.1.0 on a new Debian install. I also have alsa 1.0.27.1 installed. Debian is running under VirtualBox 4.2.12 under OSX 10.6.8, on a MacBookPro. libao does not work in stereo for the sampling rate of 44100, but it does for 8000, 11025, 22050and 48000. When I send it two channel output at a sampling rate of 44100, the result includes what sounds like skips,

Well, one way to do it is xyplot(y~time|id, data=dat,type='l', panel=function(x,y,subscripts,...){ panel.xyplot(x,y,subscripts,...) panel.abline(v=dat[subscripts,"cm1"]) panel.abline(v=dat[subscripts,"cm2"]) } ) Since I don't know what the dataframe 'mergeData' is I have used 'dat'. Best regards Frede Aakmann

Consider this data frame: dat<-data.frame(id=gl(3,5),time=rep(1:5,3),cm1=rep(c(2,3,4),each=5),cm2=rep(c(2.5,3.5,4.5),each=5),y=rnorm(15)) > dat id time cm1 cm2 y 1 1 1 2 2.5 -1.0824549 2 1 2 2 2.5 -0.7784834 3 1 3 2 2.5 -1.7783560 4 1 4 2 2.5 0.5056637 5 1 5 2 2.5 -0.1967505 6 2 1 3 3.5 1.1643136 7 2 2 3 3.5 -1.4323765 8 2

Dear WizaRds, I tried to construct a two-phase sampling design in Survey just the way I hoped understood in Vienna - I was wrong. I think I am too stupid to create the correct subset for phase 2. Phase1: Sample 1000 parts with 80 defective. Phase2: Sample 100 parts out of these 1000 with 15 defective. 0:ok, 1:defunct. The table below gives the conditional sampling values. Please help me:

Hi I've 3 iax connections to my provider , each of them have own DID , PH1<----| | \/ PH2<-->|-----| <---------------------------> |----|<-- DID1 | A1 | <---------------------------> |ISP |<-- DID2 PH3<-->|-----| <---------------------------> |----|<-- DID3 I had iax phone on each of this connection , but now I want to terminate all

I have a simple cdr configured with the default tables, here is a row of a good cdr report calldate | clid | src | dst | dcontext | channel | ect ..... ect .... 2008-08-29 10:16:49 | "C. BOUTON" <40> | 40 | XXXXXXXXXXX | phonesystems | SIP/40-08776938 | ect ..... ect .... I have replaced the number by

Hola buenas noches, tengo un problema bastante gordo. ¿A alguno le ha dejado de funcionar las funciones pamer.fnc y mcp.fnc con la nueva versión de R? La semana pasada formatee el ordenador y ahora scripts antiguos no funcionan. La cuestión es que me precupa que no funcione el ejemplo de tutorial del autor. Os dejo un script que debería de funcionar y no lo hace

Dear friends: It's possible make a legend in the coplot output? My data are: "dados"<-structure (list ( Qu1 = c(1,2,3,4,5,6,7,8,3,2), Qu2 = c(3,4,5,6,7,8,9,9,6,3), Ku1 = structure (factor ( c(1,1,1,1,2,1,2,2,1,1),levels=1:2),.Label=c("Kmin","Kmax")), Ph1 = structure (factor (

Hi all, I want to build a central call diverter via asterisk (http://www.voip-info.org/tiki-index.php?page=Asterisk%20cmd%20SIPredirect). Calls would come in via SIP from Cisco Callmanager Asterisk would do some searching an diverts the call to an extension, which is also located at Callmanager. When like to use >>sipredirect<< Asterisk complains "No application

Hi, I'm trying to figure out how anova works in R by translating the examples in Sokal And Rohlf's (1995 3rd edition) Biometry. I've hit a snag with planned comparisons, their box 9.4 and section 9.6. It's a basic anova design: treatment <- factor(rep(c("control", "glucose", "fructose", "gluc+fruct",

I have an NT4 server as PDC for a domain: WG. Domain members (on 3 subnets) can either log into the domain, or at least access shares from the network neighboorhood if they don't want to log on. When I boot a samba server as a workgroup member (I don't log it into the domain) everything runs fine for about 24 hours; domain or WG members can access shares on any machine including the

Hi, I want to get the performance data like documents <<Xenoprofoverview&Networking Performance Analysis>> (http://www.xensource.com/files/xs0106_xenoprof.pdf ). But how can I get the modules'' profile in kernel? Now I use command: Opcontrol --no-vmlinux --separate=kernel,cpu --start-daemon ...run benchmark Opcontrol --sutdown Opreport >list The report shows: CPU:

Hola, Lo primero muchas gracias a todos por vuesta habitual ayuda. Llevo meses con las funciones lmer, pero tengo un problema no se como plantear unas comparaciones multiples una vez demostrado que el factor tiene diferencias significativas. Tengo multiples opciones, Solución A: una es hacer un bonferroni a saco (pero eso me chirria ya que tego muchos datos) Solución B: Otra es utilizar glht,

Dear R-Help I am running apply on a data.frame containing factors and numeric columns. It appears to convert are columns into as.character? Does it convert data.frame into matrix? Is this expected? I wish it to recognise numerical columns and round numbers. Can I use another function instead of apply, or should I use a for loop in the case? > summary(xmat) A B