similar to: Deviance explined in GAMM, library mgcv

Dear, I would like know the way to calculate the variance explained of a regression tree. I use the function "ctree" from library "party" many thanks [[alternative HTML version deleted]]

<minor rant> All right. I've seen this. I mentioned it along with a few other 2.2.0 bugs I've seen on samba-ntdom along with some possible 2.0.7 bugs related to printing, and even a few possible security issues (not repeated here), one I know I sent directly to a samba developer. I admit that I might be a bit off topic at times, but no one tells me where to go. I seem to

I have the following xts objetct "temp" > str(temp) An ?xts? object from 2010-12-26 to 2011-03-05 containing: Data: num [1:70, 1] 2.95 0.852 -0.139 1.347 2.485 ... - attr(*, "dimnames")=List of 2 ..$ : NULL ..$ : chr "t_n" Indexed by objects of class: [POSIXct,POSIXt] TZ: GMT xts Attributes: NULL > temp t_n 2010-12-26

Hi All, I would like to ask the following. I have an array of data in an objetct, let's say X. I need to use a for loop on the elements of one or more columns of X and I am having a debate with a colleague about the best memory management. I believe that if I do: col1 = X[,1] col2 = X[,2] ... colx = X[,x] and then for(i in whatever){ do something using col1[i], col2[i] ... colx[i] }

Hi, Does anybody knows how to coerce an mcmc objetct to a rv object without loose the information about what parameters. was simulated I tried do it but it looks that this informatio was lost. Thank you very much for any information. Gilberto Matos.

Hi all R-Experts, I'm facing one problem in yahoo data downloading. I'm suing Windows XP, R 2.2.0, and i'm using yahoo.get.hist.quote function to download data. I need 500 companies of S&P index daily 'closing price' data for last ten years. My questions are: 1) I have all the ticker names of S&P 500 companies in a .csv format. I'm reading those names in R and

Hi r-users, I want to use the sn package but I got the following message: > install.packages(repos=NULL,pkgs="c:\\Tinn-R\\sn_0.4-12.zip") Warning: package 'sn' is in use and will not be installed updating HTML package descriptions I did tried to do it a few times but it gives the same message.  ________________________________ From:

Dear all I am trying to install qtiplot in Centos 7 (qtiplot is a graph plotting software similar to Origin). I followed the same step as this mention on this. http://poolinfo.physik.hu-berlin.de/qtiplot-installer/README.html At step 8, I experience the following error. cd fitPlugins && make -f Makefile make[1]: Entering directory `/home/rahman/Downloads/qtiplot-0.9.8.9/fitPlugins'

A maybe trivial and stupid question: In the case of a lm or glm fit, it is quite informative (to me) to have a look to the null deviance and the residual deviance of a model. This is generally provided in the print method or the summary, eg: Null Deviance: 658.8 Residual Deviance: 507.3 and (a bit simpled minded) I like to think that the proportion of deviance 'explained' by the

Dear R users, I read your methods of extracting the variance explained by each predictor in different places. My question is: using the method you suggested, the sum of the deviance explained by all terms is not equal to the deviance explained by the full model. Could you tell me what caused such problem? > set.seed(0) > n<-400 > x1 <- runif(n, 0, 1) > ## to see problem

Un texte encapsul? et encod? dans un jeu de caract?res inconnu a ?t? nettoy?... Nom : non disponible URL : <https://stat.ethz.ch/pipermail/r-help/attachments/20100731/ac3b9e43/attachment.pl>

Hi Everyone, I'm working on building some classification trees, and up to this point I've been using rpart. However, I recently discovered the tree package, and found that it had some useful functions (in particular deviance(), which I would really like to use for my project). I can't seem to find an equivalent function for rpart. I've considered using tree() in place of

I'm running a categorical data analysis with a two-way design of nominal by ordinal structure like the Political Ideology Example (Table 9.5) in Agresti's book Categorical Data Analysis. The nominal variable is Method while the ordinal variable is Quality (Bad, Moderate, Good, Excellent). I rank/quantify Quality with another variable QualityR (1, 2, 3, 4), and run the following:

It seems that the standardized deviance residulas, that one gets on plots of a glm.object x with plot(x) are calculated as r <- residuals(x) s <- sqrt(deviance(x)/df.residual(x)) w <- weights(x) hii <- lm.influence(x)$hat r.w <- if (is.null(w)) r else (sqrt(w) * r) rs <- r.w/(s * sqrt(1 - hii)) This implies that, for example, for binomial B(ni,pi) data the devaince residials

Hi all, I'd like to fit different models using a loop together with the jags package. to do this, I load the following packages "runjags", "R2jags" and "coda" and use the function "jags" to fit my models. The problem is that the following a windowx with the following message appear between each model fitting inside my loop : the following objects in

I have used packages rpart, mvpart and tree for classification and regression trees. I want to calculate fraction of null deviance explained by each node and variable in the tree. For instance, at the first split, this would be (1 - (sum of residual deviance in each of the two leaves)/deviance at the root). In the subsequent splits, this formula is slightly different. There probably is a function

Hi, I am doing a regression tree using the package 'rpart' but could not able to calculate the residual mean deviance. Please help. Narayan [[alternative HTML version deleted]]

Dear all, I would like to get glm() return its results (at least the deviance) for some given parameter values (ie without actually fitting the model). I tried to set `maxit = 0' but this does not work, eg: > glm(y ~ x, start = c(1, 1), maxit = 0) Error in glm.control(...) : maximum number of iterations must be > 0 Any idea? Thanks in advance. Emmanuel Paradis

dear List: glm(a~b+c,family=binomial,data=x)->fit deviance(fit) returns the same as the residual deviance. I don't not know much about logistic regression.Some book tells that: " Deviance (likelihood ratio statistic): Deviance = -2log( likelihoodof the currentmodel /likelihoodof thesaturated model) Note: (1). The current model is the model of interest. (2). The saturated model

Dear all, I'm trying to fit a statistical model to series of measurements. Unfortunately, my knowledge of statistics is rather limited, so I'm a bit at loss of what is going on with the model. First of all, I've prepared a histogram. Then, I've tried to fit a Poisson model to express the relation between the middle points of classes (mids) and the corresponding frequencies