similar to: [OT] 'gv' and fractional points

Hi, I want to draw a scatter plot with 1M and more points and save it as pdf. This makes the pdf file large. So i tried to save the file first as png and than convert it to pdf. This looks OK if printed but if viewed e.g. with acrobat as document figure the quality is bad. Anyone knows a way to reduce the size but keep the quality? /E -- Dipl. bio-chem. Witold Eryk Wolski MPI-Moleculare

Hello, i am searching now the R-Documentation many hours, but i cannot find a solution to my problem. I hope, i can solve my issue with R. i have data-tables, 15*40 data-tables i want to visualise this. this example-code i copied makes a 3x3-grid: x <- 0:12 y <- sin(pi/5 * x) op <- par(mfrow = c(3,3), mar = .1+ c(2,2,3,1)) for (tp in c("p","l","b",

I have two CSV files (exported from Excel), say file1 and file2. The have the same number of rows, and each has several columns, with names on the first line; and some of the columns in file1 are repeated in file2. Using the "foreign" package, I can read these in separately to dataframes say d1 and d2 with > d1<-read.csv("file1") >

The formula attribute of the builtin CO2 dataset seems a bit strange: > formula(CO2) Plant ~ Type + Treatment + conc + uptake What is one supposed to do with that? Certainly its not suitable for input to lm and none of the examples in ?CO2 use the above.

Hi Folks, I don't think I'm up to implementing a function for drawing French Curves, as requested by "dream". However, the discussion about it took me back so vividly to the old days that, paranormally, I felt once again the urge to ascertain magnitudes as it used to be, and should be, done. I have therefore implemented, and hereby donate to the R community, a new function

The following can't be right, first rw2010: > 1 %% 0.001 [1] 0.001 Then rw2001: > 1 %% 0.001 [1] -2.081668e-17 > and the last seems about right. -- Kjetil Halvorsen. Peace is the most effective weapon of mass construction. -- Mahdi Elmandjra -- No virus found in this outgoing message. Checked by AVG Anti-Virus.

Hi again folks, I seem to have got the trellis (lattice + grid) business working basically OK, in that I can get my 5x5 array of histograms laid out as they should be. However, there is a behaviour I can't see my way round. Even when displayed in the R Graphics Window, there is some clipping (the extreme top, bottom, left and right edges are not there, so that labels, and parts of numerical

Can someone clarify what I cannot make out from the documentation? The function 'clogit' in the 'survival' package is described as performing a "conditional logistic regression". Its return value is stated to be "an object of class clogit which is a wrapper for a coxph object." This suggests that its usefulness is confined to the sort of data which arise in

Hi Folks, I'm interested in looking (in a beginner's amateurish way) at prosodic/phonetic analysis of recorded speech. In particular I would want to use R to formulate and evaluate specific models. So I would like to ask R people for their recommendations for a program which would a) Take as input a sound file in one of the common formats (".wav", ".au") b)

Hi Folks, I note that, while the "chisq" functions dchisq(x, df, ncp=0, log = FALSE) pchisq(q, df, ncp=0, lower.tail = TRUE, log.p = FALSE) qchisq(p, df, ncp=0, lower.tail = TRUE, log.p = FALSE) rchisq(n, df, ncp=0) all have a slot for the non-centrality parameter "ncp", of the functions for the t and F distributions: dt(x, df, log = FALSE)

Greetings Folks, When R code (as entered or read from a courced file) is executed, is it interpreted from the input form every time having once been read in, or do subsequent invocations use an "intermediate" (pre-interpreted) form? Or, putting it another way, is the execution of R code faster second time time round (and later) because the pre-interpretation has already been done once

Hi My simulation with the followin R code works perfectly: sim <- replicate(999, sum(exp(rgamma(rpois(1,2000), scale = 0.5, shape = 12)))) But now I do not want to have values in object "sim" exceeding 5'000'000, that means that I am just using the beginning of densitiy function gamma x < 15.4. Is there a possibility to modify my code in an easy way? Thanks for any help!

Hi Folks, I'm seeking advice about how to resume an outer loop following failure of a function within the loop (which issues an error message). Essentially, I'm repeating (simulating) a process which involves random sampling, EM and MCMC. I'm walking on very edge of rather thin ice -- data rather thinly spread over many parameters! Occasionally some component of the loop fails, with

Hi folks, I have a vaiable Y and an associated factor Z at several (13) levels. boxplot(Y~Z) produces a nice array of boxplots, one for each level of Z, and each duly labaelled with its level of Z. I would like to superpose on each boxplot the actual data points which it represents, i.e. do something conceptually (though not in real R) expressed as points(Y~Z) or points(Z,Y) It can

Hi Folks, While I was browsing in the R-help archives yesterday, I got curious about the time series of the sizes of the monthly archives in MB. This turned out to have an unexpected feature or two, which I leave to readers to explore for themselves. I'm now wondering at what point in time we might expect to be receiving 1000MB/month (30+MB/day). It's not that far away, it seems, but

Sorry Folks, I'm sure I could suss out the answer myself but I need it soon ... ! 1. Given a set of 4 variables X,Y,Z,W in a dataframe DF, I make a scatter-plot matrix using splom(DF). 2. I do all regressions of U on V using lm(U~V), where U and V are all 12 different ordered pairs from X,Y,Z,W. 3. Now I would like to superpose the regression lines from (2) onto the corresponding

Hi folks, Suppose I have a series of cases each with categorical factors A, B. What is the best way to "glue" A and B together into a single factor? For example, given A0 B1 ... A1 B1 ... A0 B2 ... A1 B0 ... A0 B0 ... A1 B2 ... then I'd like to end up with a single factor with levels A0B0, A0B1, A0B2, A1B0, A1B1, A1B2 according to all the combinations which actually occur in

Hi Folks, What is the best way to avoid a function being read in anew (and masking an exiting function) when a definition of it has already been established in R? Reason: Fernando Tusell and I are working up Schafer's 'CAT' for R (basically done now, just needs some cosmetic tidying up). This uses a function 'slice.index', present in S but not in the versions of R we were

Doese anyone know a R function to find the median of a gamma distribution? [[alternative HTML version deleted]]

I am looking for a way to produce a "distribution graph" as in the example: (http://cecsweb.dartmouth.edu/release1.1/datatools/dgraph.php?year=2003&geotype=STD_HRR&event=A01_DIS&eventtype=UTIL Anybody who can help? Christian von Plessen Department of Pulmonary Medicine Haukeland university hospital Bergen Norway