search for: pubhealth

...unced message when I tried to submit a bug from http://bugs.r-project.org/cgi-bin/R (B) The bug report itself is at the bottom. Thanks, Chris (A) -------- Original Message -------- Subject: Undelivered Mail Returned to Sender Date: Wed, 6 Aug 2008 17:08:02 +0200 (CEST) From: MAILER-DAEMON at pubhealth.ku.dk (Mail Delivery System) To: candrews at buffalo.edu This is the Postfix program at host mail.pubhealth.ku.dk. I'm sorry to have to inform you that your message could not be delivered to one or more recipients. It's attached below. For further assistance, please send mail to <p...

Hi, I am new to R and I'm looking to perform some method comparison analysis and would be grateful for some help regarding package installation. There is a package 'MethComp' that I want to use (http://staff.pubhealth.ku.dk/~bxc/MethComp/Courses/Sthm.2007/.) Unfortunately I am unable to install it. I have tried numerous approaches but to no avail. this is the first time I have tried to install a package 'outside' of R. My operating system is Windows XP. Here's one attempt that does not work: where=...

...y to use is to split follow-up time so finely that you can safely assume that rates are constant in each interval, and then just stuff it all into a Poisson model. This allows you to use any kind of time-dependent variables as well as accommodating competing risks. In the Epi package (http://staff.pubhealth.ku.dk/~bxc/Epi/) on CRAN there is a machinery for representation and manipulation of multistate data (of which competing risks are a special case), see the function Lexis(). It also has a time-splitting function there, splitLexis(). If you take a look at practical 14 from this year's SPE cours...

Hi, The following function works, but is there a neater way to write it? f = function(x,...) { # return a character vector of the arguments passed in after 'x' gsub(" ","",unlist(strsplit(deparse(substitute(list(...))),"[(,)]")))[-1] } > f(x,a,b,c*d) [1] "a" "b" "c*d" > Thanks. [[alternative HTML

...than a matrix because I have columns of different types in real life. > rownames(d) = NULL Error in "dimnames<-.data.frame"(`*tmp*`, value = list(NULL, c("a", "b" : invalid 'dimnames' given for data frame -----Original Message----- From: pd at pubhealth.ku.dk [mailto:pd at pubhealth.ku.dk] On Behalf Of Peter Dalgaard Sent: 08 December 2005 18:57 To: Matthew Dowle Cc: 'r-help at stat.math.ethz.ch' Subject: Re: [R] data.frame() size Matthew Dowle <mdowle at concordiafunds.com> writes: > Hi, > > In the example below why is...

Hi folks, I'm having a spot of bother with tempdir(). If I set the environment variable TMPDIR in my shell prior to invoking R it works as expected: > Sys.getenv("TMPDIR") TMPDIR "/eurotempest/tmp/R" > tempdir() [1] "/eurotempest/tmp/R/Rtmp0xY4XD" However if I don't set it prior to entry, but attempt to set it from within

There is now an Epi-package on CRAN. It is intended for epidemiological analysis in R. It has its own homepage, http://www.pubhealth.ku.dk/~bxc/Epi The package has been used at the course "Statistical practise in Epidemiology with R", see http://www.pubhealth.ku.dk/~bxc/SPE. A mailing list for people using R for epidemiology, R-sig-Epi, has been set up, see https://stat.ethz.ch/mailman/listinfo/r-sig-epi or through...

Hi all. Has anyone ever done a Deming Regression in R? I'm wondering if there's a simple way to do it. Thanks for all your help! -Ed -- View this message in context: http://www.nabble.com/Deming-Regression-tp17949318p17949318.html Sent from the R help mailing list archive at Nabble.com.

...lem in the most appropriate way? convolve2=function2(x,y) { x<- c(rep.int(0,exp.length-1),x) n <- length(y<-c(y, rep.int(0, irf.length -1))) x <- fft(fft(x) * Conj(fft(y)), inverse=TRUE) return(Re(x)/n) } -- Best wishes/bedste hilsner Gunnar Hellmund cand. scient., PhD http://staff.pubhealth.ku.dk/~guhe/ --- ?If everyone is thinking alike, then somebody isn't thinking.?

There is now an Epi-package on CRAN. It is intended for epidemiological analysis in R. It has its own homepage, http://www.pubhealth.ku.dk/~bxc/Epi The package has been used at the course "Statistical practise in Epidemiology with R", see http://www.pubhealth.ku.dk/~bxc/SPE. A mailing list for people using R for epidemiology, R-sig-Epi, has been set up, see https://stat.ethz.ch/mailman/listinfo/r-sig-epi or through...

Hi, under the new version R 2.7.1 I can't use WinEdt. It seems that there is a bug in that version. To use older versions of the package RWinEdt doens't bring success. The Use of R 2.7.0 doesn't cause that Problem. Regards, Hans-J?rgen Auinger -- ---------------------------------------------- Hans-J?rgen Auinger Lehrstuhl f?r Pflanzenz?chtung Wissenschaftszentrum Weihenstephan

Hello, I have a dataframe which I read from a file: df <- readtable(myFile); The dataframe has 4 columns: "model", "count", "value" and "date" where "model" and "date" are made of charactors and "count" and "value" are numbers. An example is like: model count value date A 4 20.8

I have a dataframe with 700,000 rows and 2 vectors (columns): ?group? and ?score?. I wish to calculate a third vector of length 700000: the average score by group. Even though the avarge value will repeat, I wish to return the average for that particular group for each row. (I know I can do this by calculating each group?s average and then using the merge command, but as my calculations get

I am producing 2X2 correlation matrices by a class variable. I need to extract a vector of correlation coefficients only. I am doing that in a loop (see below) but I am sure there would be a simpler way. Please help! > by(d1[,c(2,3)],d1[,1],cor) d1[, 1]: 1 c e c 1.0000000 0.1972309 e 0.1972309 1.0000000

Hi, In the example below why is d 10 times bigger than m, according to object.size ? It also takes around 10 times as long to create, which fits with object.size() being truthful. gcinfo(TRUE) also indicates a great deal more garbage collector activity caused by data.frame() than matrix(). $ R --vanilla .... > nr = 1000000 > system.time(m<<-matrix(integer(1), nrow=nr, ncol=2)) [1]

Dear List: When I click on the link to download a reference manual for a package on cran, I get an error message that the file is damaged and could not be repaired. I randomly chose various packages and the same error message appears. Are the links actually broken? I have also restarted my machine and closed and re-opened acrobat. I am using Windows XP, Acrobat Professional 6.0.0.5, and

All, When I take a subset of a factor the reduced factor still maintains all the original levels of the factor when say forming the key in a plot. The data is correct, but the variable still "remembers" the original levels. See below for reproducible code. Does anyone know how to fix this? cheers, dave fact = as.factor(c(rep("A", 3),rep("B", 3), rep("C",

Dear R-helpers Apologies for the basic question, but I just got stuck: I would like to write values from a vector into array cells with the same names > count[1:10] 10010 10014 10015 10017 10030 10080 10100 10230 10250 10280 0 0 0 0 0 1 1 0 2 0 >data[1:10,,1] [,1] [,2] [,3] [,4] [,5] 10010 NA NA NA NA NA 10014 NA NA NA

Hi Folks! I'm reading in some data from a .csv file that has a date column. How do I use colClasses to get read.csv to recognize the date column? The documentation on this seems to be nil - And yes, I've read help and R Data Import/Export and can't figure out what the colClasses syntax is. Thanks, john [[alternative HTML version deleted]]

Dear list, sometimes my function generates too much data and shows them on screen, i cannot view first several lines until program ends and I have to scroll my mouse up to get them. Is there any re-direction function in R to pipeline outputs to "more" or "less" type functions? Thanks mike [[alternative HTML version deleted]]