similar to: How to generate 'minor' ticks in lattice (qqmath)

Dear group members, I want to compare response variables ("logAUC") of two groups (treatment "Test", "Reference") of a subset ("period == 1") in dataframe "resp" (below): sequence subject period treatment AUC logAUC 1 RT 1 1 Reference 44.1 3.786460 2 RT 1 2 Test 39.1 3.666122 3 RT 2

Hi, I am reading the source code of rpart. I have problems understand the following code and would appreciate for any helps. In rpart.s, there is a line: rpfit <- .C(C_s_to_rp,             n = as.integer(nobs),             nvarx

Hi, Don't have a clue what teh following error message generated by this function call: qqmath( ~val|ind,data = xx ,distribution = function(p){ qt(p,df=20)} ,ylab="Sample Quatinles" ,xlab="Theoretical Quantiles" ,panel=function(x,y) { panel.qqmathline(y , distribution=function(p) qt(p,df=20)

Having solved this problem, I am posting this so that the next time I search for how to do this I will find an answer... Using qqmath(..., groups=num) creates a separate qq distribution for each group (within a panel). Using the 'col' or 'pch' argument does not (usually) work because panel.qqmath sorts the data (but not 'col' or 'pch') before plotting. Sorting

Does the 'groups' option on qqmath just color the points differently in the main distribution or does it actually overlay separate quantile plots for each subset? I would like to be able to do the latter. --Rich Richard Kittler AMD TDG 408-749-4099

In qqmath, how would one go about having 'distribution' change with panel.number? I've tried set.seed(1) mydata <- data.frame(ind = factor(rep(2:4, each = 100))) mydata$val <- rt(300, df=rep(2:4, each = 100)) plot<-qqmath(~ val | ind, layout=c(3,1), data = mydata, prepanel = function(x, distribution, ...) {

David, Duncan, > Hi > > Following on David's rate argument > > try (with modifications of pch and grid) > > rate <- 1/4 > shape = 8 > rate = c(rep(1/4,100),rep(1/3,100)) I don't think the problem is related to the rate argument, which can well be vectorized, as is the case for a number of arguments in distrib functions in R (note that you are redefining it

Hello, I can't seem to make qqmath plot text. v <- rnorm(8) qqmath(~v) # this is ok, but qqmath(~v, panel=function(x, y) { points(x, y) text(0, 0, 'print me', cex=2)}) isn't ok. The axis frame is printed properly (ie, as in qqmath(~v)), but neither 'print me' nor the points show up on the graph. I can, however, see the upper half of the word 'me' in the

Hello, I would like to inlude the Q-Q plot by "qqmath" into a panel with other plots, say, using par(mfrow=c(1,2)). How can this be done given that "qqmath" refreshes the plotting window and there seems to be no series coming out of it? Thanks Serguei [[alternative HTML version deleted]]

The data are not normally distributed when untransformed and I'm trying various transformations to see if any would be appropriate to use. The lattice book (fig. 3.10) shows a 2-sample Q-Q plot with an abline but the code for the figure does not include the line. I'd appreciate a pointer to a reference on how to add an abline to a one-sample qqmath() plot in lattice. Rich

Hello ALL! Can anyone tell me how to specify qqmath() function with distribution qgamma? In help it is bit vague how to pass shape and scale parameters for gamma. Thank you in advance. Best, PM

I want to write a function that get 2 strings y and z and does the following R command. temp<-qq1[qq1$z==y,] for example if it get y="AMI" and z="PrimaryConditionGroup" It should do the following temp<-qq1[qq1$PrimaryConditionGroup=="AMI",] I could do it by the following function that is ugly and I wonder if there is an easier way to do it espacielly when temp

Dear R users, I have a follow-up question on sqlSave(). Since most of the output from the tests I use are lists, I would like to loop to export each element of the list and append it to the sheet. Here is what I do: > library(RODBC) > test <- structure(list(m = structure(c(0.090909090909091, 0.181818181818182, 0.272727272727273, 0.363636363636364, 0.454545454545455,

Dear R list, The code below puts qq-plots for two of three groups on the one plot. However the legend includes all three groups, ie the auto.key ignores the subset instruction. Is there an easy way to get around this, so that only those groups plotted are included in the legend? y1<-rnorm(100); y2<-rnorm(100)+1; y3<-rnorm(100)+2; y<-c(y1,y2,y3)

This is the qqmath example from the lattice package. I added the scales to the example. I would like to switch the axis and not sure how? Meaning I would like the "height" on the x-axis and the probability on the y-axis. Will you show me the correct syntax for this switch thanks. qqmath(~ height | voice.part, aspect = "xy", data = singer, prepanel =

Dear collegues, Is there anyone who ever used qqmath(). Actually, I want to know the distribution of the dependent varible. It's a count number, so it should be poisson or negative binomial distribution. But how could I use R to get the ë and á? Or to demonstrate the distribution to be negative binomial ? I tried my best to find the examples by google, but failed. Thank you for the

Hey R Users I like to control the ticks and labels in a boxplot as described for a xyplot below (thread in maillinglist in may 2003). Does anybody knows how it works? Thanks in advance Thomas Thread from May 2003 (http://tolstoy.newcastle.edu.au/R/help/03a/5604.html) Hello R Users! I'm using lattice to produce some graphs with logaritmic y-scales. I use the command xyplot(hits ~

Hello R Users! I'm using lattice to produce some graphs with logaritmic y-scales. I use the command xyplot(hits ~ c(1:1024), data=eichData, type="S", scales=list(y = list(log=10))) to create the plot. This is fine, except for the automatically choosen tick marks. I'd like to have a major tick at the 10^n location and minor ticks in between which correspond with the native

Full_Name: Bill Wheeler Version: 2.0.1 OS: Windows 2000 Submission from: (NULL) (67.130.36.229) The program to reproduce the error is below. I am calling rpart with a user-defined split function for a binary response variable and one continuous independent variable. The split function works for some datasets but not others. The error is: Error in "$<-.data.frame"(`*tmp*`,

The Pygments highlighter supports many more languages than CodeRay. Pygments is written in Python, but I''ve seen non-Python projects use it (Jekyll, Pygments). There''s also this: https://github.com/nex3/rb-pygments Any idea how much work it would be to add Pygments support to Webgen?