similar to: how to make plotmath expression work together with paste

Dear R Helpers: I am trying to annotate a plot. The following code snippet works, but it is kind of a kludge since it adds the partial derivative symbols after creating the plotmath frac(). Is there a more elegant way to write a partial derivative? dev.off() plot(NA, xlim=c(-3,3), ylim=c(0,1.6), xlab="", ylab="", tck=-0.015) text(1.6, 1, expression(paste("slope =

I want to write something like (LaTeX style) b_{norm}=\frac{F\prime(0)}{R\prime(0)} how do I add the "prime" (first derivative) to a R-plot? The help of "plotmath" just talks about "partialdiff". Can you complete this command? text( 30,0.05,labels=expression(b[plain("norm")]==frac(F(0),R(0))) ) Thanks, Thomas

I want to shade the area under the curve of the standard normal density. Specifically color to the left of -2 and on. How might i go about doing this? Thanks -- View this message in context: http://n4.nabble.com/Shade-area-under-curve-tp1586439p1586439.html Sent from the R help mailing list archive at Nabble.com.

Greetings: I would like comments on this example and after fixing it up, I need help from someone who has access to insert this in R's help page for plotmath. I uploaded a drawing http://pj.freefaculty.org/R/Normal-2009.pdf that is created by the following code http://pj.freefaculty.org/R/Normal1_2009_plotmathExample.R This will be a good addition to the plotmath help page/example.

Number of subjects = 25 Mean of Sample = 77 Standard Deviation (s) = 12 sem = 2.4 df = 24 The claim is that population mean is less than 80 * > 80 So our H0 (null hupotheis) is * > 80 > qt(.95,24) [1] 1.710882 > qt(0.05, 24) [1] -1.710882 tstat = -1.25 on t24 falls between 1.711 (.95,24) and *1.711 (.005,24) How Could I sketch t curve for the above data where my * would be at the

Dear list, I got an automatic email complaining than some of my packages didn't pass 'R CMD check' for R-2.10.0 alpha anymore. Both of them make use of the "cases" LaTeX environment. Inspecting the log outputs I got: LaTeX errors when creating PDF version. This typically indicates Rd problems. LaTeX errors found: ! Misplaced alignment tab

Dear R-help, I'm using prop.test() to compute a confidence interval for a proportion under R version 2.4.1, as in: prop.test(x = 340, n = 400)$conf [1] 0.8103309 0.8827749 I have two questions: 1) from the source code my understanding is that the confidence interval is computed according to Wilson, E.B. (1927) Probable inference, the law of succession, and statistical inference. J. Am.

the version 2 parser thinks I have unnecessary braces, but I can't find any. False positive or am I missing something? If a false positive, is there any way to work around the warning? * checking Rd files against version 2 parser ... WARNING Warning: ./man/dbetabinom.Rd:32-34: Unnecessary braces at ?{p(x) = % (C(N,x)*Beta(N-x+theta*(1-p),x+theta*p))/% Beta(theta*(1-p),theta*p)}?

Dear expeRts, How can I vertically adjust an axis tick label so that it is nicely aligned with the other labels? library(lattice) xyplot(0~0, xlim=c(0,3), scales=list(x=list(at=c(1,1.1), labels=c(expression(hat(theta)[italic(n)]),expression(theta))))) ## aim: move the leftmost expression up so that theta is nicely aligned with the second theta. Cheers, Marius

Hello, I have the following question, which I haven't been able to resolve after days of trying. I used to save my plots as jpegs using the savePlot command. However, that seems to result in lost resolution. So now I'm trying to use the jpeg( ) function, but am having trouble because it seems to be incompatible with the windows (width=, height=) command. It's important for me to

Dear R-users, I would like to calculate elasticities and sensitivities of each parameters involved in the following transition matrix: A <- matrix(c( sigma*s0*f1, sigma*s0*f2, s, v ), nrow=2, byrow=TRUE,dimnames=list(stage,stage)) The command "eigen.analysis" avaliable in package "popbio" provides sensibility matrix and elasticity matrix

I spent more time than I should have debugging a script because I wanted x<-seq(0,100)*0.1 but typed x<-seq(O:100)*0.1 seq(0:100) yields 1 to 101, Clearly my own brain to fingers fumble, but possibly one others may want to avoid it. JN

I hope to use the plotmath facility to print titles that mix math and values of R variables. The help for "plotmath" has an example, which after repeated reading, I find baffling. Likewise, I have read the help file for "substitute" (wqhich seems to be needed) without ever understanding what it does, other than being used in some magic incantations. I would like to do

Dear HelpeRs, I would like to do sg similar to: plot(c(0,100), c(0,100), xaxs="i", yaxs="i", type="n") txt <- paste("a =", a, "%\n b =", b, " km2", sep="") text(95, 95, txt, adj=c(1,1)) just with the "km2" formatted with 2 in a superscript. I thought txt <- substitute("a = "*a*" %\n b

hello, I have the following problem setting parameter 'frac' in the sigest function of the kernlab package. ## executing the ?sigest example: library(kernlab) data(spam) srange <- sigest(type~.,data = spam) ## works fine... ## setting 'frac' explicitly ## (in this case even to the default of .25) options(error=recover) srange <- sigest(type~.,data = spam, frac = .25) ##

Hi to all I use a loop to plot 9 different histograms of 9 different transformations of one dataset (x). I want to label the histograms with the mathematical expression of each transformation (e.g. x^3). For that I've prepared a vector with the labeling names for "expression". > trans.expr <- c("x^3", "x^2", "x", "frac(1,x)",

The following line of code seems fairly straight forward, yet it is kicking back an error message: for (i in 1:length(mean.natveg.frac)){month.observed[i]=as.numeric(names(mean.natveg.frac[i]))%%12} Error message: Error in month.observed[i] = as numeric(names(mean.natveg.frac[i]))%%12 : object 'month.observed' not found Seems like its not found because I'm trying to define it. Also

... or if I understand correctly, simply expression(frac(additive ~ HCO[3]^"-", true ~ HCO[3]^"-" ))) Cheers, Bert On Fri, Feb 2, 2024 at 3:06?AM Rui Barradas <ruipbarradas at sapo.pt> wrote: > ?s 10:01 de 02/02/2024, Troels Ring escreveu: > > Hi friends - I'm plotting a ratio of bicarbonates i ggplot2 and > > > >

?s 10:01 de 02/02/2024, Troels Ring escreveu: > Hi friends - I'm plotting a ratio of bicarbonates i ggplot2 and > > ylab(expression(paste(frac("additive BIC","true BIC")))) worked OK - but > now I have been asked to put the chemistry instead - so I wrote > >

I am using both nlminb and optim to get MLEs from a likelihood function I have developed. AFAIK, the model I has not been previously used in this way and so I am struggling a bit to unit test my code since I don't have another data set to compare this kind of estimation to. The likelihood I have is (in tex below) \begin{equation} \label{eqn:marginal} L(\beta) = \prod_{s=1}^N \int