Displaying 20 results from an estimated 2000 matches similar to: "Computing plot size in Sweave"
2008 May 18
1
Figure environment and includegraphics options from Sweave
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}
2008 Jul 29
2
FW: Installing BRugs
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
2018 Mar 07
3
Names of variables needed in newdata for predict.glm
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 =
2009 Apr 27
2
The .tex version of the manual in foo.Rcheck
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.
2007 Jun 11
2
Rounding?
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
2018 Mar 08
0
Names of variables needed in newdata for predict.glm
Hi,
Some try:
> names(mi$xlevels)
[1] "f"
> all.vars(mi$formula)
[1] "D" "x" "f" "Y"
> names(mx$xlevels)
[1] "f"
> all.vars(mx$formula)
[1] "D" "x" "f"
When offset is indicated out of the formula, it does not work...
Marc
Le 07/03/2018 ? 06:20, Bendix Carstensen a ?crit?:
> I would like
2004 May 24
2
Month names
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"
2004 Nov 03
1
Building a package under WIN2000 / rw2.0
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
2004 Nov 03
1
Building a package under WIN2000 / rw2.0
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
2008 Nov 24
1
FW: read.ssd
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
2010 Apr 17
1
Fishy error with NAMESPACE when checking package
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:
2011 Jul 16
2
Sweave in 2.13.1
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
2005 Feb 17
1
How to get interction terms first in a model
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
2007 May 11
1
Compilation of source package.
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:
2005 May 21
2
print format for difftime
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
2004 Jan 08
3
Strange parametrization in polr
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
2004 Mar 02
3
Margins on tables
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
2004 Sep 27
1
Funny behaviour of coef() and vcov() if X is singular
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
2018 Mar 31
1
Names of variables needed in newdata for predict.glm
all.vars works fine, EXCEPT, it give a bit too much.
I only want the regression variables, but in the following example I also get "k" the variable holding the chosen knots. Any machinery to find only "real" regression variables?
cheers, Bendix
library( splines )
y <- rnorm(100)
x <- rnorm(100)
k <- -1:1
ml <- lm( y ~ bs(x,knots=k) )
mg <- glm( y ~
2006 Jun 29
2
zero.print in print.table after adding margins
The function addmargins() adds margins to a table, but returns a matrix.
But even after converted to a table the print.zero="." option of
print.table() does not work:
> x <- sample( 1:7, 20, replace=T )
> y <- sample( 1:7, 20, replace=T )
> tt <- table( x, y )
> tx <- as.table( addmargins( table( x, y ) ) )
> print( tt, zero.print="." )
y
x 1 2