similar to: Modes again

Hello! Up to now I have been using hist() to display the distributions. Howevere, I noteiced strange numbers on y (vertical) axis, if I used probability = T or freq = F option. I thought it is a bug and launched the R-bug system and found some posts on that matter. Brian Ripley responded to one, that one should look at truehist() for that. Ok I can use truehist() if I want to see the ratios

Hello! I often plot average of y instead of all y values, since one can easily see the trend if there is to many points and/or x might be like 1, 2, 3, 4, 5, ... and you might get a cloud (which can also be informative) and/or columns of points. Anyway, learning with R i get stucked. I have the data in table burnins1, with the following names: names(burnins1) [1] "Model"

Duncan, you are a king! Thanks a lot for this cookie. It really helped me. Thanks for the code as well as detailed explanation at the end. >Hi Gregor. > >Here is a function that will collect all of the nodes in the >XML document whose names are in the vector elementNames > >getElements = >function(elementNames) >{ > els = list() > > startElement = function(node,

I unintentionally missed to cc to r-devel. -------- Original Message -------- Subject: Re: [Rd] levels for list and data.frame Date: Tue, 21 Mar 2006 20:50:21 +0100 From: Gregor Gorjanc <gregor.gorjanc at gmail.com> Reply-To: gregor.gorjanc at gmail.com To: Martin Maechler <maechler at stat.math.ethz.ch> References: <441F2C49.8080703 at gmail.com> <17439.50351.172900.987055

James, thanks for the response. I understand now my puzzle. tmp <- data.frame(y1=1:4, f1=factor(c("A", "B", "C", "D"))) tmp$y2 <- NA tmp[1:2, "y2"] <- 2 Does the job. I see that I should add a full column. In my case adding full column of NAs and then adding values, solves my problems. Thanks to all. james.holtman at convergys.com

Jun hi again! Sorry for bothering. Although I had no problem with building Windows package under Linux, I noticed that man pages are not processed as they should or maybe I have a problem with my setup. My package is called 'GeneticsPed'. Instead of HTML or text files I just got a GeneticsPed.Rd.gz in man directory of a package (this is within GeneticsPed_0.1.zip file!). The file

Hello to everyone! I have built CRAN task view for genetics. For now I have not submit it to CRAN yet and it can be accessible from: http://www.bfro.uni-lj.si/MR/ggorjan/software/R/Genetics.html http://www.bfro.uni-lj.si/MR/ggorjan/software/R/Genetics.ctv I have not submitted it to CRAN, since I would like first some opinion about it. Genetics is really so broad field that I belive one person

G?ran, thanks. Happy easter! -----Original Message----- From: G?ran Brostr?m [mailto:gb at stat.umu.se] Sent: ned 2005-03-27 07:38 To: Gorjanc Gregor Cc: r-help at stat.math.ethz.ch Subject: Re: [Rd] F90 On Sun, Mar 27, 2005 at 03:26:39AM +0200, Gorjanc Gregor wrote: > Hello! > > Has anyone successfully compiled F90 sources in R-package? I found the > same question on r-devel list

And a last post from Paul Gilbert. Thanks to all! This disscusion was really beneficial for me! -----Original Message----- From: Paul Gilbert [mailto:pgilbert at bank-banque-canada.ca] Sent: sre 2005-06-08 21:01 To: Gorjanc Gregor Subject: Re: [R] Random seed problem in MCMC coupling of chains Gorjanc Gregor wrote: > Thanks to Paul and Gabor for additional tips/examples. Actually, I find

On 12 October 2006 at 13:02, Gregor Gorjanc wrote: | Hello! | | I encounter the following problem when I try to load package foreign | | library(foreign) | Error in dyn.load(x, as.logical(local), as.logical(now)) : | unable to load shared library '/usr/lib/R/library/foreign/libs/foreign.so': | /usr/lib/R/library/foreign/libs/foreign.so: undefined symbol: | Rf_allocString | Error:

Hello to Rusers! I am puzzled with R and I really do not know where to look in for my problem. I am moving from SAS and I have difficulties in translating SAS to R world. I hope I will get some hints or pointers so I can study from there on. I would like to do something like this. In SAS I can write a macro as example bellow, which is afcourse a silly one but shows what I don't know how

Hi r-help forum, I have been collecting contact data (with proximity logger collars) between a few different species of animal. All animals wear the collars, and any contact between the animals should be detected and recorded by both collars. However, this isn't always the case and more contacts may be recorded on one collar of the two. This is fine, it depends on battery life and other

A sample from a continuous distribution most likely has no meaningful sample mode. Shown below is a way to estimate the location of mode from a kernel density estimate (using the locfit package): > x <- rnorm(50) > library(locfit) Locfit for R. August 3, 2000. (Updated for R 1.7.0, March 21, 2003) Attaching package 'locfit': The following object(s) are masked from

Hello! Is it possible to set na.rm=TRUE in a global way? I'am constantly forgeting on this when performing analyses. I agree that one should be carefull with this when developing some code, but not necesarilly so in data analysis. Lep pozdrav / With regards, Gregor Gorjanc ---------------------------------------------------------------------- University of Ljubljana PhD student

Hi! I'm running an anova on a lm and I get the message Error: invalid dimnames given for data frame Any idea where it comes from? Many thanks. Here is the fit and the attempted anova: > fit Call: lm(formula = as.formula(paste(chemical, "a~", cmd, sep = "")), singular.ok = TRUE) Coefficients: (Intercept) Odom Displace Weight

I am having trouble fitting correlation structures within nlme. I would like to fit corCAR1, corGaus and corExp correlation structures to my data. I either get the error "step halving reduced below minimum in pnls step" or alternatively R crashes. My dataset is similar to the CO2 example in the nlme package. The one major difference is that in my case the 'conc' steps are

Hello! I am crossposting this to R-help and BioC, since it is relevant to both groups. I wrote a wrapper for Entrez search utility (link for this is provided bellow), which can add some new search functionality to existing code in Bioconductor's package 'annotate'*. http://eutils.ncbi.nlm.nih.gov/entrez/query/static/esearch_help.html Entrez search utuility returns a XML document

Dear all, When analyzing data from a climate change experiment using linear mixed-effects models, I recently came across a situation where: - the summary(model) showed a significant difference between the levels of a two-level factor, - while the anova(model) showed no significance for that factor (see below). My question now is: Is the anova.lme() approach correct for that model? And why does

Hello R-users! I have tried to use Sweave within LyX* and found two ways to accomplish this. I have attached LyX source file for both ways. *<http://www.lyx.org> Lep pozdrav / With regards, Gregor Gorjanc ---------------------------------------------------------------------- University of Ljubljana Biotechnical Faculty URI: http://www.bfro.uni-lj.si/MR/ggorjan Zootechnical

Hello! I am writing some functions and I repeatedly acces internal factor codes. I figured out that internal codes are 1:n where 1 represents 1st level, 2 2nd level etc. This is not documented and I wonder if this is on purpose and subject to change or would be a good idea to add this to documentation for factor? I can prepair a pacth if this is accepted. -- Lep pozdrav / With regards,