similar to: Problem with greek "beta" symbol in plot using CairoPDF (Windows)

Hi All: I have some basic questions about Cairo graphics engine. I'm trying to use the Cairo package to produce PDF output, mainly because I perceive it to be easy to use with a wide variety of fonts. But right now, I'm stuck trying to figure out what fonts are available to be used with Cairo, specifically the CairoPDF function. I've been able to successfully produce some test PDFs

Hello there, Straight to the point: it seems that CairoPDF from package "Cairo" cannot handle greek letters from expression(). For example, > eta = seq(from=-pi, to=pi, length=100) > f = sin(eta)^2 > pdf(file = "temp_pdf.pdf") > plot(eta, f, type="l", main=expression(f(eta)==sin(eta)^2), xlab=expression(eta), ylab=expression(f(eta))) > dev.off()

Dear list members I'm using CairoPDF to generate PDF (because of its font embedding and support for transparent colours). However, at least on my (Windows) system, the text it outputs seems to have completely wrong kerning. Here's an example: CairoPDF("test.pdf") plot(rnorm(100),xlab="Ovreset") dev.off() The v is (slightly) too far away from the O, it's much

I have a small problem with R graphics output. When I use the lattice package and CairoPDF to generate publication quality graphs I often use the expression to create an axis title that has microlitres or micrometers as a unit. I use something like the following 'expression(paste("Length (", mu,"m )"))' as an argument to the xlabel function. The command works but the mu

Dear R-helpers, This may sound simple to you, but I'm a beginner in this, so please be forgiving. I have a following problem: two analytes were measured in patient's blood on 4 occasions: ProteinA and ProteinB. How to correctly evaluate correlation between ProteinA and ProteinB? I tried: x <- data.frame(Patient.ID=rep(1:10, each=4), Visit=rep(c(1:4),10), ProteinA=rnorm(m=10,

Hello, I have made an R script that prepares a NONMEM dataset and I would like to start the NONMEM run right after the dataset is ready. I am using windows XP, R 2.9.1 and NONMEM 6. I have prepared a run.bat file that looks like this: ---------------------------------------- call K:\nmvi\NMdirectories.bat call K:\nmvi\nmfe6 "path\control.txt" "path\output.txt"

Dear R-helpers, I am working with 'cph' models from 'rms' library. When I build simple survival models, based on 'Surv(time, event)', everything is fine and I can make nice plots using plot(Predict(f, time=3)). However, recently I tried to be more specific and used 'Surv(start, stop, event)' type model. Using this model 'plot(Predict(f))' works OK, but

Dear R users/developers, simple comparison of code execution time of R 2.7.1 and R 2.7.2 shows a dramatic slowdown of the newer version. Rprof() identifies .Call function as a main cause (see the code below). What happened with R 2.7.2? Kind regards Marek Wielgosz Bayes Consulting ######### Probably useful info ############### ### CPU: Core2Duo T 7300, 2 GB RAM ### WIN XP ### both standard

Dear R-helpers I have a data set of roughly 10 million records, 7 columns. It has only about 500MB as a csv, so it fits in the memory. It's painfully slow to do anything with it, but it's possible. I also have another dataset of covariates that I would like to explore - with about 4GB of data... I would like to merge the two datasets and use lmer to build a mixed effects model. Is

Dear R-Helpers, I'm working with immunoassay data and 5PL logistic model. I wanted to experiment with different forms of weighting and parameter selection, which is not possible in instrument software, so I turned to R. I am using R 2.14.2 under Win7 64bit, and the 'nls' library to fit the model - I started with the same model and weighting type (1/y) as in the instrument to see

Dear R-helpers I am running a simple mixed effects model using lme(). The call looks like this: fit <- lme(Analyte~Sample, data=Data, random=~1 | Run) I am particularly interested in the estimated random effects. When I print the 'fit' object, it looks something like example below: (...) Random effects: Formula: ~1 | Run (Intercept) Residual StdDev: 3.483794 3.637523

Dear R-Helpers, I am a novice in survival analysis. I have the following code: for (i in 3:12) print(coxph(Surv(time, status)~a[,i], data=a)) I used it to fit the Cox Proportional Hazard models separately for every available parameter (columns 3:12) in my data set - with intention to compare the Hazard Ratios. However, some of my variables are in range 0.1 to 1.6, others in range 5000 to

Dear All, I encountered the following problem: fun1 <- function(x){ 2*cos(x) } >library(debug) > mtrace(fun1) > fun1(0.1) Error in all.levs[[j]] : subscript out of bounds > sessionInfo() R version 2.5.0 (2007-04-23) i386-pc-mingw32 locale: LC_COLLATE=Polish_Poland.1250;LC_CTYPE=Polish_Poland.1250; LC_MONETARY=Polish_Poland.1250;LC_NUMERIC=C;LC_TIME=Polish_Poland.1250 attached

Dear R-helpers and graphics gurus, I have two problems with plotting confidence bands: 1. First is relatively simple. I am using the Passing-Bablok procedure to obtain "unbiased" regression coefficients. This procedure yields the "a" & "b" coefficient values along with their confidence intervals. I then plot the raw data with the regression line, but I would

Dear R-helpers, I am developing a Mixed-Effects model for a study of immunoassays using 'lme4' library. The study design is as follows: 10 samples were run using 7 different immunoassays, 3 times each, in duplicates. Data attached. I have developed the following model: c.lme <- lmer(Result~SPL + (SPL|Assay/Run) -1, data=data) This model has excellent predictions - the Rsquared of

Hi, I ran into the following peculiarity involving S4 method dispatch and namespaces. I develop a package which have a namespace and which depends on 'pixmap' package, which itself does not have a namespace. Now, in my package I have a class which has a slot for objects from class "pixmap" and I want to have a "plot" method for this new class. Not to clutter the

Dear useRs and wizaRds, I am currently developing a set of functions using S4 classes. On the way I encountered the problem exemplified with the code below. For some reason the 'validity' method does not seem to work, i.e. does not check for errors in the specification of the slots of the defined class. Any hints? My understanding of the whole S4 system was that validity checks are made

Hello everyone, I am new to R, I have been using SAS for a while. Not surprisingly, I find R much better in graphics, which is publication ready right away. Recently, I have been trying to calculate some basic statistics using R. I have a dataset of multiple rows per subject. For example: Subject Date Factor1 Factor2 Factor3 P1 0.5 1 1 3 P1 1 3 2 5 P1 2 3 5 NA ... P2 0.5 1 6 4 P2 1 2 NA 7 P2

I am trying to install the Cairo package on a linux machine, the Cairo package did not install correctly (could not find cairo.h), I am new to R and linux any help on the installation would be great. Below are output from trying to install the Cairo package, thought this might help. > install.packages("Cairo") Warning in install.packages("Cairo") : argument

Dear all, The CRUK MI have multiple positions in and around bioinformatics, biostatistics and computational biology. These range from HPC and mathematical programming, through biostatistics and biomarkers, to computational biology and bioinformatics. Focuses include single cell genomics and personalised medicine in lung and prostate cancer. For full details, please visit our website at: