similar to: lattice: combining panel.xyplot with panel.abline - is thispossible?

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

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:

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 =

Hi, I have a data.frame with many variables for which I am performing the mean by subgroup, for a pair of variables at a time, where one of them for each pair defines the subgroup. The subgroups in the x$cm1 are 0, 1 and 2. x ph1 cm1 0.2345 2 1.2222 1 2.0033 0 0.0000 2 1.0033 1 0.2345 0 1.2222 2 2.0033 0 0.0000 1 1.0033 2 > meanbygroup <- as.vector(with(x, by(x$ph1, x$cm1, mean)))

Hi all, How can one evaluate NAs in a numeric dataframe column? For example, I have a dataframe (demo) with a column of numbers and several NAs. If I write demo.df >= 10, numerals will return TRUE or FALSE, but if the value is "NA", "NA" is returned. But if I write demo.df == "NA", it returns as "NA" also. I know that I can remove NAs, but would like

Hi, I ran two svm models in R e1071 package: the first without cross-validation and the second with 10-fold cross-validation. I used the following syntax: #Model 1: Without cross-validation: > svm.model <- svm(Response ~ ., data=data.df, type="C-classification", > kernel="linear", cost=1) > predict <- fitted(svm.model) > cm <- table(predict,

hi, sorry for the potentially stupid question. I'm a little stuck. when trying to run the stupid update manager i need for my satnav, which - obviously - is windows only, nothing happens. install worked, the icon's on my desktop, it has been made executable, but after a bit of harddrive action nothing happens. terminal says: Code: mzwo at mzwo-laptop:~$ env

Hi, I ran two svm models in R e1071 package: the first without cross-validation and the second with 10-fold cross-validation. I used the following syntax: #Model 1: Without cross-validation: > svm.model <- svm(Response ~ ., data=data.df, type="C-classification", > kernel="linear", cost=1) > predict <- fitted(svm.model) > cm <- table(predict,

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:

Is that possible to create superscript text on the graph legend, for example to put "cm2" (centimeter square) on the legend. Please show me how to do it. Thanks. -- Agus Susanto [[alternative HTML version deleted]]

Hi I've suddenly lost access to some permanent cross-server shares. I think this was following the Nov 5 new kernel, rather than the Centos version upgrade. I don't see any difference to man mount.nfs I have:- [root at hs6 ~]# mount hs8:/data /sysback/hs8-data -v mount: no type was given - I'll assume nfs because of the colon mount: trying 172.26.0.8 prog 100003 vers 3 prot tcp port

I am trying to fit a nonlinear regression to infiltration data in order to determine saturated hydraulic conductivity and matric pressure. The original equation can be found in Bagarello et al. 2004 SSSAJ (green-ampt equation for falling head including gravity). I am also VERY new to R and to nonlinear regressions. I have searched the posts, but am still unable to determine why my data come up

I thought this had been fixed with newer kernels? We have an old Maxtor Snap Server, an 80 GB NAS unit, containing some of our archives. Both of our CentOS servers have it permanently mounted by NFS for convenient access. This weekend, I discovered remotely that, although all systems except the VRTSralus remote agents for Veritas (now Symantec) Backup Exec were working just fine, both servers

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'm working with large data frames and running out of memory. I hope some of you may be able to suggest a more efficient approach. I have grid/lattice data representing a time series of 1 m2 quadrats in a grassland: Each 1 cm2 cell or pixel contains one ecological state (ie grass or bare ground). The goal is to calculate, for each cell, the transition probabilities to all available

Hi I get this on every scan report:- LibClamAV Warning: *********************************************************** LibClamAV Warning: *** This version of the ClamAV engine is outdated. *** LibClamAV Warning: *** DON'T PANIC! Read http://www.clamav.net/support/faq *** LibClamAV Warning: *********************************************************** ----------- SCAN SUMMARY

I'm just wondering if anyone else here is using this? The majority of jobs on our CentOS5 server "fail" by "Communications failure" right at the end, apparently on the very last file which, however, is recoverable, as is the whole archive. In more than half these failures, the Linux agent has stopped and, in about half of those, crashes leaving a backtrace. I had similar

Dear all, I'm facing a problem in estimation of glm model with weibull distribution. I run this : eqn0<-formula(fdh~cup1+cup2+cup3+cup4+fin1+vd1+cm2+cm4+milieu+cpro1+cpro2+cpro3a+cpro3b+schef+log(y)) regWeib0<-vglm(eqn0,family=weibull,subset(br, fdh<1)) I have en estimation but there is a message saying that regularity conditions are violated : Message d'avis : In eval(expr,

An embedded and charset-unspecified text was scrubbed... Name: not available URL: <https://stat.ethz.ch/pipermail/r-help/attachments/20091125/bf6c327e/attachment-0001.pl>