Search Archives

I notice a lot of r documentation refer to this reference below. I can't seem to find it anywhere. Does anyone have a link to point to where I can either view it or buy it? *Chambers, J. M...

Run this code tree<-graph_from_literal(1-+2:3,3-+5,1-+4); graph.bfs(tree,root=1, neimode="out",father=TRUE,order=TRUE,unreachable = FALSE) I do not understand why the father values will give NA 1 1...

Hi *I have been looking for a reference to explain how R uses the aov command(at a deeper level)*. More specifically, how R reads the formulae and R computes the sums of squares. I am not interested in understanding...

Below are samples from a kernel density estimated "data" with gaussian kernel. I really like this solution of estimation of a kernel because it is nice and elegant. fit<-density(data) rnorm(N, sample(data, size = N, replace = TRUE), fit...

In https://en.wikipedia.org/wiki/Pareto_distribution, it is clear what the parameters are for the pareto distribution: *xmin *the scale parameter and *a* the shape parameter. I am using rmutil to generate random deviates from a pareto distribution...

Dear R package community, I am uncertain whether this is appropriate for this mailing list. Please let me know. If not, would you be so kind as to point me in a better direction? I am a mathematics major with...

Something does not make sense in R. It has to do with the question of balance and unbalance. *A<-factor(c(1,1,1,1,1,1,2,2,2,2,2,2))* *B<-factor(c(1,1,2,2...

Hi I need some advice. Note: I do not know anything from C apart from my 2 days of research. I am currently trying to make meaning of the modelmatrix function (written in C) and called from R function model...

Hi Peter, Thank you for your good answer. I am sorry for the late reply. *An ortogonalized model matrix generates a decomposition of the model space into orthogonal subspaces corresponding to the terms of the model. Projections onto each of...

Hi again I need help *R-code* debug(model.matrix) model.matrix(~S) *model.matrix code* ans <- .External2(C_modelmatrix, t, data) #t =terms(object) , data="data frame of object" *modelframe C-code* SEXP modelframe(SEXP call, SEXP op, SEXP...

I have something which has been bugging me and I have even asked this on cross validated but I did not get a response. Let's construct a simple example. Below is the code. A<-gl(2,4) #factor of...

Hi, Thanks Richard, That was me playing with too many examples and having too many variables just lying around. Thanks for the tip though. On 22 August 2016 at 23:32, Bert Gunter <bgunter.4567 at gmail.com> wrote: > Thanks...

Hi Michael, Thank you for the reply. I am sorry I forgot to print out the anova table to make my question clear. Df Sum Sq Mean Sq F value Pr(>F) S 2 0.000019 9.630e-06 0...

Below is the covariates for a model ~x1+x2+x3+x4+x5+x6. I noticed that when fitting this model that the coefficient x6 is unestimable.*Is this merely a case that adding more columns to my model matrix will...

Why does the S:x1 column disappear (presumably S:x1 goes into ID but I dont know why)? S is a factor, x1 is a covariate and ID is a factor. rich.side<-aov(y~S*x1+ID) summary(rich...