similar to: Bug?: summary() fails after use of na.action="na.exclude" in lm()

Hello! I have already reported on this behaviour at [1] and now I have checked the same example on r-devel and the problem is still there. Thomas Lumley inidicated problem [2]. Example: n <- 50 x <- runif(n=n) y1 <- 2 * x + rnorm(n=n) y2 <- 5 * x + rnorm(n=n) y2[sample(1:n, size=5)] <- NA y <- cbind(y1, y2) ## Goes ok here fit <- lm(y1 ~ 1,

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

Hello! I posted on saturday mail with the same subject on r-help seeking for help in my work, but now I realized that this list is more appropriate for this. I think I found I bug. Bellow are comments and reproducible examples: # Create a data frame (tmp <- data.frame(y1=1:4, f1=factor(c("A", "B", "C", "D")))) y1 f1 1 1 A 2 2 B 3 3 C 4 4 D #

Hello! I posted on saturday mail with the same subject on r-help seeking for help in my work, but now I realized that this list is more appropriate for this. I think I found I bug. Bellow are comments and reproducible examples: # Create a data frame (tmp <- data.frame(y1=1:4, f1=factor(c("A", "B", "C", "D")))) y1 f1 1 1 A 2 2 B 3 3 C 4 4 D #

Hello R users! I have written one function (look at the end), which will ease my work with analysis of data in another programme, for which I need sometimes a special data structure. However I encountered several problems with a created data frame. --------------------------------------------------------------- The data frame (produced from the example at the end) looks like the way I want

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

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

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,

Hello! na.encode does not have any effect on format of NA values of Date and POSIXct (POSIXlt?) "atomic" classes in a data.frame. Here is the example (the same in R 2.3.1 and 2.5.0 (2006-09-19 r39409)): testData <- data.frame(num=c(NA, 2.6), int=c(1, NA), fac=factor(c(NA, "abc")), cha=c("a",

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

Hi! If I use format() on numeric vector, na.encode argument does not have any e= ffect. This was reported before: - https://stat.ethz.ch/pipermail/r-help/2007-October/143881.html - http://tolstoy.newcastle.edu.au/R/e2/devel/06/09/0360.html It works for other (say character) classes! > format(c("a", NA), na.encode=3DTRUE) [1] "a " "NA" >

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

Hello R users! I just found out that I have 49 invitations for Gmail (gmail.google.com). I have been using it now for a while and is really nice. Don't forget 1 GB for free. I will invite those who respond to this mail by FIFO. -- Lep pozdrav / With regards, Gregor Gorjanc ------------------------------------------------------------------------ University of Ljubljana Biotechnical

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 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

Hello! I have a continous distribution and would like to get mode (the most frequent value in distribution). I easily found mean, median and other basic thing but not mode function. Can anyone help? I know there my might be problems with multiple modes, but still I think that there should be a mode function in R. Please send mail to R-help list and me, so I can get response faster. Thank

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

Hello! Has anyone successfully compiled F90 sources in R-package? I found the same question on r-devel list from 2002 and I wonder if there is any progress. I heard that g95 and gfortran can be usable. And one more thing. I have a program written in F90 and lets suppose that compiler is not a problem. I've read manual 'Writing R Extensions' and successfully loaded one F77 subroutine

Hello! Has anyone successfully compiled F90 sources in R-package? I found the same question on r-devel list from 2002 and I wonder if there is any progress. I heard that g95 and gfortran can be usable. And one more thing. I have a program written in F90 and lets suppose that compiler is not a problem. I've read manual 'Writing R Extensions' and successfully loaded one F77 subroutine

Hello! Can someone provide me an estimate or relevant source of information on how long is command winepath available or from wich version of wine? I would need this as I work on some R[1] package R2WinBUGS[2] that calls a Windows program WinBUGS[3] and that program needs windows paths. Currently winepath or parsing of .wine/config is used and I wonder if I can remove fiddling with .wine/.config.