similar to: Specifying Prior Weights in a GLM

With ggplot2, I can plot the glm stat_smooth for binomial data when the response is binary or a two-level factor as follows: data("Donner", package="vcdExtra") ggplot(Donner, aes(age, survived)) + geom_point(position = position_jitter(height = 0.02, width = 0)) + stat_smooth(method = "glm", family = binomial, formula = y ~ x, alpha = 0.2, size=2) But how can I

Hi all, Sorry to bother you with the potentially not too bright question, but since I cannot get it to work, I'll have to ask someone.. I've got this vector: > token [1] sal skal skal bak b_r d_ gal kal l_r n_r pak p_r [13] sal sl_r sn_r spak sp_r st_r s_ bak d_ gal l_r pak [25] p_r sak sl_r sp_r t_ kal n_r s_ st_r sak kal

Dear R Users I am maximizing a user defined log likelihood function. It includes variance parameter (sigma). I used R function optim with BFGS maximization method. However, it stops before the solution saying ?sqrt(sigma): NaNs produced? Could anybody know a proper transformation for sigma which can be passed in the function? For the correlation parameter I used Fishers? transformation so it

X.Org Security Advisory: October 25, 2023 Issues in X.Org X server prior to 21.1.9 and Xwayland prior to 23.2.2 ===================================================================== Multiple issues have been found in the X.Org X server implementation published by X.Org for which we are releasing security fixes for in xorg-server-21.1.9 and xwayland-23.2.2. The first issue (CVE-2023-5367) can be

Hi all- I would like to enable search in my application via Ferret and a number of my searchable classes contain dates (of type Time). When these classes are indexed, the dates are written per the standard to_s output (e.g. "Wed Dec 07 08:56:04 CDT 2006"), which is not suitable for date range queries. Is it possible to format the dates before they are indexed (e.g.

I''m working my way through the AWD book and have just started the testing chapter. The first two unit tests run as expected. When I add the test for the ''destroy'' function to product_test.rb, however, I''m getting a failure I don''t understand. The error I''m getting says that in the setup method for test_update, it "couldn''t

Hi, In Winbugs, one can assign a distribution called "dflat()", which corresponds to an improper uniform prior on the whole real line. What would be the equivalent in R? I ask because I'm trying to run Winbugs through R, and the model I've constructed in Winbugs includes a "dflat()" distribution. For all of the normal priors in Winbugs, I've been using

Hello, I'm trying to use Mclust to fit a Gaussian Mixture Model to a mulitdimensional data set. Because of the specific source of my data, I know that all components have the same variance and that the covariance between dimensions is zero (modelname=VII). Furthermore, I have a reliable estimate of the variance of the components. I want to to use this estimate as a prior in mclust, hoping

I'm sorry, this is probably a very simple question. I am relatively new to R and am having difficulty understanding how to use the mclust options to change the prior probabilities for a mixture model. I am dealing with a situation where I have 2 groups of data and it is known that the sample sizes are different in the population. I want to incorporate the probability associated with

Hey, R-listers, I am now going to use Bayesian mathod to estimate a matrix parameter C. It is assumed that C is an orthogonal matrix already. We know, if C is an arbitrary column vector, we may use multivariate Gaussian prior on it. However, now it is a matrix, so what can I do to define a proper prior probability density function for C? Thanks for your point. Have a nice day! Fred

Hi R folks, How can I generate the discriment function from lda? I have an unbalanced data set. one class has about 25 entries and another class has about 200 entries. I used lda for classification > z<- lda(V3 ~ V1+V2, data) > z Prior probabilities of groups: 0 1 0.91111111 0.08888889 Group means: V1 V2 0 0.4445161 0.04723951 1 0.4058900

Give this one a try which will require the user/hosts permissions to be setup prior to starting postgres --- installer/modules/ovirt/manifests/postgres.pp | 2 +- 1 files changed, 1 insertions(+), 1 deletions(-) diff --git a/installer/modules/ovirt/manifests/postgres.pp b/installer/modules/ovirt/manifests/postgres.pp index c46b360..381d67c 100644 ---

Hola, si tengo estas dos funciones de distribucion y quiero en una grafica verla, como le hago en R? [image: Imágenes integradas 1] Cual es el comando en R para graficar estas funciones? Saludos, Tania ------------ próxima parte ------------ Se ha borrado un adjunto en formato HTML... URL: <https://stat.ethz.ch/pipermail/r-help-es/attachments/20121215/7d185b91/attachment-0001.html>

Hi there, I have a question on the LLVM internals. Is it possible to dump an InstructionList (i.e. a (possibly) naively scheduled assembly) prior register allocation? Does LLVM use infinite (virtual) registers similar to MachSUIF? This is, of course, meant for a given target in contrast to MachSUIF that features the SUIFvm ISA as low-level IR and such a dump is possible at this point. Plus:

As requested by Chris. Neil. -------------- next part -------------- Index: include/llvm/ADT/APInt.h =================================================================== --- include/llvm/ADT/APInt.h (revision 41148) +++ include/llvm/ADT/APInt.h (working copy) @@ -19,9 +19,7 @@ #include <cassert> #include <string> -#define HOST_CHAR_BIT 8 -#define compileTimeAssert(cond) extern int

Hello all, As I'm sure is pretty common, I have some extensions that dial mobile numbers after a local timeout. I would like to prompt the caller to record their name after the local timeout and have the recipient be able to hear the name prior to accepting the call. Recording the message is easy enough, so I thought about doing something like dumping them into MeetMe after they record

In another thread I talked about running CI testing on a packge before release. I thought I would explain that a bit more AND ask for tests, if they are needed. Prior to pushing newly released packages into the main CentOS os, updates, extras, cr, or fasttrack repositories, we run our t_functional CI test suite on the entire repo set. So if we are releasing an update, it would be against the os

On Wed, 3 Apr 2019 12:33:55 +0100 Stephen via samba <samba at lists.samba.org> wrote: > Hi everyone, i've just been following the instructions about setting > up a Samba domain member as a file shares. I am using Samba version > 4.5.16-Debian (yes its old, but i'm stuck with it for now ;) ) and I > have been following the official Samba documentation found here: >

hello, I am using LDA and QDA function of MASS library. I understand Fisher LDA is a method non-probabilistic in nature, so I wonder what happens when I try to predict my test set examples as in: > fit <- lda(labels~., data=train.table, prior=c(.5,.5)) > pred <- predict(fit, data=test.table, prior=c(.5,.5)) Specifically I ask this because in my problem there are 700 examples

Hi! Could you please explain the difference between "prior" and "weight" in rpart? It seems to be the same. But in this case why including a weight option in the latest versions? For an unbalanced sampling what is the best to use : weight, prior or the both together? Thanks a lot. Aur?lie Davranche.