similar to: frequencies

A funny thing happened when I wanted a student of mine to install Brugs. Using the InstallPackages in the windows version, firts gives an erro, but trying again works flawlessly. R version is 2.7.0 on WinXP. Any explanation? Bendix Carstensen ______________________________________________ Bendix Carstensen Senior Statistician Steno Diabetes Center Niels Steensens Vej 2-4 DK-2820 Gentofte

I was a bit puzzed by: > formatC(6.65,format="f",digits=1) [1] "6.6" So I experimented and found: > formatC(6.6500000000000001,format="f",digits=1) [1] "6.6" > formatC(6.650000000000001,format="f",digits=1) [1] "6.7" > round(6.6500000000000001,1) [1] 6.7 > round(6.650000000000001,1) [1] 6.7 > version

In version 2.8.1, running Rcmd check on the package foo would leave the file foo-manual.tex in the folder foo.Rcheck. But as of 2.9.0 only foo-manual.pdf and foo-manual.log are there. Is this intentional? Anyway it is inconvenient, because I would occasionally like to include the manual at the end of a set of exercises, and this was a convenient file to \input with a few select %'s added.

Did not seem to reach Saikat DebRoy, this might be the forum. Bendix -----Original Message----- From: BXC (Bendix Carstensen) Sent: 24. november 2008 15:00 To: 'saikat at stat.wisc.edu'; 'stvjc at channing.harvard.edu' Subject: read.ssd It's always annoyed me that that read.ssd crashed on datasets with long variable names, but the other day a collegue of mine pointed out to

Sometimes you want to compute the physical size of a plot based on data. In R itself this is no problem. But is there a way to compute the values of height and width in S-weave, say: <<graph,fig=TRUE,height=xx,width=yy>>= where xx and yy are computed and not physically written in the document? Bendix ______________________________________________ Bendix Carstensen Senior

This is how I get the month names from within R: > mon <- rep(strptime("01/01/1952", format = "%d/%m/%Y"), 12) > mon$mon <- mon$mon + 0:11 > mnam <- months(mon, abbreviate = F) > mnam [1] "januar" "februar" "marts" "april" "maj" "juni" "juli" "august"

I have an odd problem in building a package with only R-code in it. I have a package mainly used by myself which I last build under R 1.9.0. The operation system is Win2000 5.00.2195, Service Pack 3 When I do: c:\stat\r\rw2000\bin\Rcmd install --docs=normal --build --library=c:\stat\R\bxc\library c:\stat\R\bxc\library.sources\xx then after updating help pages I get: preparing package xx for

I have an odd problem in building a package with only R-code in it. I have a package mainly used by myself which I last build under R 1.9.0. The operation system is Win2000 5.00.2195, Service Pack 3 When I do: c:\stat\r\rw2000\bin\Rcmd install --docs=normal --build --library=c:\stat\R\bxc\library c:\stat\R\bxc\library.sources\xx then after updating help pages I get: preparing package xx for

Dear all, I import my data from a csv-file containing one row with date-entries. How can I tell R to treat this data as dates? I've tried to bring it in character-format (as.charachter()) followed by as.POSIXlt() but I get an error message, that the character string is not in a standard unambiguous format, although the date is of format yyyy-mm-dd (what is standard format according to R

I am updating the Epi package. I added functions named pc.points and pc.matpoints. Erroneously I wrote pc.plot and pc.matplot in the NAMESPACE file and of course got an error from Rcmd check. So I corrected the NAMESPACE file, but I still get from r-check: * install options are ' --no-html' Loading required package: utils Error in namespaceExport(ns, exports) : undefined exports:

I run Windows XP. and in a command window I get: > c:\stat\r\R-2.13.1\bin\i386\Rcmd Sweave Lexis.rnw Error in length(arg) : 'arg' is missing Calls: <Anonymous> Execution halted AND: > c:\stat\r\R-2.13.1\bin\i386\Rcmd Sweave --help Error in length(arg) : 'arg' is missing Calls: <Anonymous> Execution halted BUT: > c:\stat\r\R-2.13.0\bin\i386\Rcmd Sweave

I have a samll package that works well and complies nicly on WinXP, using R-2.4.1 Now I want to add a document so i make a folder inst\doc and put the .tex and .pdf in there. But the complation then crashes. Is this because the installin expects some file to be present in inst if an inst folder is there? This is how the thing progress is on my command promp:

In Venables \& Ripley 3rd edition (p. 231) the proportional odds model is described as: logit(p<=k) = zeta_k + eta but polr apparently thinks there is a minus in front of eta, as is apprent below. Is this a bug og a feature I have overlooked? Here is the naked code for reproduction, below the results. ------------------------------------------------------------------------ --- version

Has anyone written a function that will print a difftime in the form: hh:mm:ss or yy-mm-dd hh:mm:ss depending on the actual size. (sloppy notation for months/minutes, but surely you get the point). Bendix ---------------------- Bendix Carstensen Senior Statistician Steno Diabetes Center Niels Steensens Vej 2 DK-2820 Gentofte Denmark tel: +45 44 43 87 38 mob: +45 30 75 87 38 fax: +45 44 43 07

Tha handy thinb about the fig=TRUE option in Sweave is that you do not have to bother about filenames and starting and stpping the device. I want the the resulting LaTeX to look as: \begin{Schunk} \begin{Sinput} > x <- seq(-2 * pi, 2 * pi, 0.1) > plot(x, cos(x), type = "l", lwd = 4) \end{Sinput} \end{Schunk} \begin{figure} \includegraphics[width=0.6\textwidth]{xx-001}

Consider the following two specifications of a model: library( splines ) x <- 1:100 y <- rnorm( 100 ) w <- rep( 1, 100 ) A <- factor( sample( 1:2, 100, replace=T ) ) B <- factor( sample( letters[1:4], 100, replace=T ) ) summary( lm( y ~ ns( x, knots=c(30, 50, 70 ), intercept=T ):A - 1 + B ) ) summary( lm( y ~ ns( x, knots=c(30, 50, 70 ), intercept=T ):A - 1 + B:w ) ) The

It has long been a nuisance to me not being able to form margins on multiway tables in a simple fashion, so i wrote margins(). In my opinion it should go into the base package. The code and the documentation is in: http://www.biostat.ku.dk/~bxc/R/margins/ Please help yourself, and enhance and rename as you see fit. Bendix ---------------------- Bendix Carstensen Senior Statistician Steno

coef() and vcov() have different dimensions if a model contains alised parameters as the following example illustrates. A search on "alised" gave noting as far as I could see. Is this a known bug? Bendix C ---------------------- Bendix Carstensen Senior Statistician Steno Diabetes Center Niels Steensens Vej 2 DK-2820 Gentofte Denmark tel: +45 44 43 87 38 mob: +45 30 75 87 38 fax: +45

Hi Say I have some data, two columns in a table being a binary outcome plus a predictor and I want to plot a graph that shows the percentage positives of the binary outcome within bands of the predictor, e.g. Outcome predictor 0 1 1 2 1 2 0 3 0 3 0

I would like to extract the names, modes [numeric/factor] and levels of variables needed in a data frame supplied as newdata= argument to predict.glm() Here is a small example illustrating my troubles; what I want from (both of) the glm objects is the vector c("x","f","Y") and an indication that f is a factor: library( splines ) dd <- data.frame( D =