On Tue, 02-Mar-2010 at 02:43PM -0500, Liaw, Andy wrote:
|> In most implementations of boosting, and for that matter, single tree,
|> the first variable wins when there are ties. In randomForest the
That still doesn't explain why with gbm, two identical variables will
"share the glory" (approxima
Dear R-Helpers,
I have a short question related to heatmap.2 function. I wanted to change
the dimension of each cell in the heatmap drawn by this function.
here is a reproducible example from the help page
library(gplots)
data(mtcars)
x <- as.matrix(mtcars)
rc <- rainbow(nrow(x), start=0, end=.
I am
working on bayesian copula but i am new in this area and i do not
have much knowledge about it. Does any one know how to use bayesian copula
using R with
winbugs? and estimating missing data by R?
hena...@yahoo.com
sincerely
[[alternative HTML version deleted]]
__
On 2010-03-05 18:00, Alastair wrote:
Hi,
I'm new to R and I've run into a problem that I'm not really sure how to
express properly in the language. I've got a data table that I've read from
a file containing some simple information about the performance of 4
algorithms. The columns are the name
## dadrivr
## I suggest a lattice plot of all two-way interactions on the
## off-diagonals and all marginal main effects on the main diagonal.
## Please install HH and then ?interaction2wt for some examples.
## install.packages("HH") ## if you don't already have it
library(HH)
?interaction2wt
m
data(warpbreaks) ## given two factors
coplot(breaks ~ 1:54 | wool * tension, data = warpbreaks,
col = "red", bg = "pink", pch = 21,
bar.bg = c(fac = "light blue"))
##
conditionalplot
#
Yes, you are missing something vital.
The log-likelihood is only defined up to an additive constant. extractAIC and
logLik use log-likelihoods with different additive constants, in general. (AIC
my look uniquely defined, but since log-likelihood itself is not uniquely
defined, nor is AIC.)
T
Try this:
> x <- matrix(0,10,5)
> vecToAdd <- 1:5
> # add to rows 1,4,8,9
> # need to transpose because of the way matrix is filled
> # this does not give the right result
> x[c(1,4,8,9),] <- vecToAdd
> x
[,1] [,2] [,3] [,4] [,5]
[1,]15432
[2,]00000
I think both graphs will follow the same technique. This is partly
laid out in how to do graphs of anova output at
http://personality-project.org/r/r.plotregressions.html (which is
part of the short (and out of date) guide to R for psychologists)
http://personality-project.org/r/
At 6:39 P
On 3/6/2010 4:38 PM, casperyc wrote:
> Hi,
>
> I am trying to reproduce a tukey test in R
>
> ==
> x=c(145,40,40,120,180,
> 140,155,90,160,95,
> 195,150,205,110,160,
> 45,40,195,65,145,
> 195,230,115,235,225,
> 120,55,50,80,45
> )
> y2=c
Hi there,
=
x=c(145,40,40,120,180,
140,155,90,160,95,
195,150,205,110,160,
45,40,195,65,145,
195,230,115,235,225,
120,55,50,80,45
)
y2=c(
rep(as.character(1),5),
rep(as.character(2),5),
rep(as.char
Phil, thanks for your suggestion of unstack, which is just what I was
looking for.
Andrew
On Fri, Mar 5, 2010 at 11:52 PM, Phil Spector wrote:
> Andrew-
> Maybe something like this:
>
>> dd = read.table(filename)
>> unstack(dd,V2~V1)
>
> A. B. C.
> 1 1 2 10.0
> 2 34 20 6.7
> 3 2 78 35.0
Hi,
I am trying to reproduce a tukey test in R
==
x=c(145,40,40,120,180,
140,155,90,160,95,
195,150,205,110,160,
45,40,195,65,145,
195,230,115,235,225,
120,55,50,80,45
)
y2=c(
rep(as.character(1),5),
rep(as.c
The usage here is exactly what I am looking for, thanks. However this function
seems to omit some combinations. Continuing with the example below, given an
always true condition in #***, it will only produce 7 combinations (omitting
1,5 1,4 and 2,5).
Am I overlooking something that makes it pr
In browsing the source I see the following construct:
res <- switch(type, working = , response = r, deviance = ,
pearson = if (is.null(object$weights))
r
else r * sqrt(object$weights), partial = r)
I understand that 'switch' will execute the code that is matched
Some code to cut and paste would be helpful. The following may help out.
library(multtest)
# create some p-values
p <- runif(100)
p <- sort(p)
p_adj <- mt.rawp2adjp(p, proc="BH", alpha = 0.05)
> str(p_adj)
List of 4
$ adjp : num [1:100, 1:2] 0.0142 0.0174 0.0254 0.0258 0.0736 ...
..- attr(*
On 02.03.2010 11:01, Sarah Paul wrote:
I am trying to load package fpc in order to use the 'plotcluster' function
however everytime I attempt to do so I get the following warning message:
> library(fpc)
Loading required package: MASS
Error: package 'MASS' could not be l
On 06.03.2010 15:58, Jordi Moya wrote:
Dear friends,
If I use:
model<-lm(y ~ x)
and then extractAIC(model), the value that I obtain does not match (not even
close):
AIC=2*k-2*logLik(model)
However, using AIC from the AICcmodavg(), the value matches exactly the
On 06.03.2010 18:35, Cleber Borges wrote:
Does anyone know if there is SVGAnnotation package for Windows, and if
so, a link to download?
A Google search reveals SVGAnnotation is an Omegahat hosted package.
http://www.omegahat.org tells us to use
install.packages(packageName, repos = "ht
I am using color2D.matplot to plot a matrix about 400 by 200.
The values in the matrix are 0:5 and NA. The resulting plot is not color, but
shaded b/w. I tried to figure out how to add colors, I would like something
like c(blue, green, red, cyan, green)
#example
motifx <- matrix(NA, nrow=10
It is not obvious to me if you can calculate f1, ..., f4, ...
automatically or have to set them manually. Looks like you need to do it
manually. In that case I'd suggest to write something along the lines:
dose <- c(-1.47, -1.1, -0.69, -0.42, 0.0, 0.42)
f1 <- function(x) exp(-x)
f2 <- function
Does anyone know if there is SVGAnnotation package for Windows, and if
so, a link to download?
Thank you!
Cleber N. Borges
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://ww
Hi,
is there a way to take a dataset and extract numeric columns and create
interaction columns from it automatically?
For e.g. there are 5 columns of data: A,B,C,D,E.
CDE are numeric.
Can someone provide code to automatically create more columns such as:
1) C*D, C*E, C*D*E, (C
Hi,
is there a way to take a dataset and extract numeric columns and create
interaction columns from it automatically?
For e.g. there are 5 columns of data: A,B,C,D,E.
CDE are numeric.
Can someone provide code to automatically create more columns such as:
1) C*D, C*E, C*D*E, (C
Hi,
I want to add a vector many time to rows of the matrix without use the loop.
could you please hint me?
thanks
khazei
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-proje
On Mar 6, 2010, at 11:41 AM, Dieter Menne wrote:
David Winsemius wrote:
Dear listeRs;
I am hoping to provide myself the ability to interpolate mortality
estimates within intervals for a range of conditions and ages. I have
four mortality tables derived from the Society of Actuaries 2001 VB
currCombin <- c(1,2) #this is the first possible combination 2 out if 5. Since
the vector has length 2, we are doing 2 out of x here.
selectedSet <- NULL
while (!is.na(currCombin)){
if(currCombin passes test) {
selectedSet <- rbind(selectedSet,currCombin)
}
currCombin <- next.combin(c
Erich-
This approach would be great for my context. However, in the code below I do
not see how to restrict the output to the set of combinations I am looking
for. For example, suppose I am looking for the 10 two element combinations of
1:5. Can you give me some psuedocode that shows how to d
David Winsemius wrote:
>
> Dear listeRs;
>
> I am hoping to provide myself the ability to interpolate mortality
> estimates within intervals for a range of conditions and ages. I have
> four mortality tables derived from the Society of Actuaries 2001 VBT
> tables that contains interval mo
Try this. sprintf formats it with a leading 0 if necessary and then
as.Date processes that:
as.Date(sprintf("%06d", c(280210, 10210)), "%d%m%y")
Note Feb only has 28 days this year so I changed the example.
On Sat, Mar 6, 2010 at 10:17 AM, ManInMoon wrote:
>
> I have text date colum with dates
Dear listeRs;
I am hoping to provide myself the ability to interpolate mortality
estimates within intervals for a range of conditions and ages. I have
four mortality tables derived from the Society of Actuaries 2001 VBT
tables that contains interval mortality rates by starting age and year
There are lots of ways to deal with text. One possibility would be to
check the length with nchar(), and if it is 5, use paste() to add a 0
to the beginning. Granted, not elegant.
Sarah
On Sat, Mar 6, 2010 at 10:17 AM, ManInMoon wrote:
>
> I have text date colum with dates such at 290210 which I
I have text date colum with dates such at 290210 which I can parse with
%d%m%y but I also have date without leading zeo i.e. 10210 (Ist Feb 2010)
which cause any NA.
How can I parse the second one?
--
View this message in context:
http://n4.nabble.com/as-date-parsing-tp1582868p1582868.html
Sent
Dear friends,
If I use:
model<-lm(y ~ x)
and then extractAIC(model), the value that I obtain does not match (not even
close):
AIC=2*k-2*logLik(model)
However, using AIC from the AICcmodavg(), the value matches exactly the
above value.
I read the help of extractAIC() and
Achim Zeileis wrote:
On Sat, 6 Mar 2010, David Winsemius wrote:
On Mar 5, 2010, at 11:54 PM, Patrick Shea wrote:
I'm trying to obtain the robust standard errors for a multinomial
ordered logit model:
mod6 <- lrm(wdlshea ~ initdesch + concap + capasst + qualrat +
terrain,data=full2)
T
On Mar 5, 2010, at 3:45 PM, Ryan Garner wrote:
>
> I just installed R on Redhat Linux at work for the first time and have two
> questions.
>
> 1. I tried to install R to have png and cairo capabilities and was
> unsuccessful. Before running make, I ran ./configure --with-libpng=yes
> --with-x=no
36 matches
Mail list logo
