search for: prior

...quot;, ylab = "Y",key = list(space = "top", text = list("B"), points = list(pch=16, col="blue", cex=0.8), text = list("A"), points = list(pch=4, col="red", cex=0.8)) ) -------------- next part -------------- x1,y,x2,x3,x4 1.05,0.550707692,Prior=0.025,m=2,A 1.1,0.719302564,Prior=0.025,m=2,A 1.15,0.846502564,Prior=0.025,m=2,A 1.2,0.927158974,Prior=0.025,m=2,A 1.25,0.958974359,Prior=0.025,m=2,A 1.05,0.508789474,Prior=0.05,m=2,A 1.1,0.674368421,Prior=0.05,m=2,A 1.15,0.846136842,Prior=0.05,m=2,A 1.2,0.9362,Prior=0.05,m=2,A 1.25,0.973957895,Pri...

Hi all, If I run LDA on the same data (2-class classification) with default(no priors specified in the lda function) vs. "prior=c(0.5, 0.5)", the results are different. The (0.5, 0.5) priors give better 1-classify-to-1 rate, and the proportional priors(default, no priors specified in the lda function) give better 0-classify-to-0 rate, for both training and testing data s...

Bugs item #11183, was opened at 2007-05-28 10:28 You can respond by visiting: http://rubyforge.org/tracker/?func=detail&atid=218&aid=11183&group_id=35 Category: Incorrect behavior Group: current Status: Open Resolution: None Priority: 3 Submitted By: Patrick McDonnell (kc9ddi) Assigned to: Kevin Smith (qualitycode) Summary: GLib errors cause crash Initial Comment: Hi - I''m working with some wxruby code, but am getting a crash with glib errors: -------------------------------- pmcdonnell at pmcdonnell10:~/workspace...

...n of the covariates, I have been trying to use the newdata argument to survfit with no success. Here is my model statement, the newdata and the ensuing error. What am I doing wrong? > summary(fit) Call: coxph(formula = Surv(Start, Stop, Event, type = "counting") ~ Week + LagAOO + Prior.f + cluster(interaction(Station, Year)), data = data8, method = "breslow", x = T, y = T) n= 1878 coef exp(coef) se(coef) robust se z p Week 0.00582 1.01 0.0323 0.0239 0.244 8.1e-01 LagAOO 0.71929 2.05 0.1238 0.1215 5.918 3.3e-09 Pri...

Dear R help, What am I doing wrong here? when I don't specify the priors it works just fine but when I specify the priors it breaks.? Does anyone know why and how I can fix it? ---- > N=20000 > ncontrol=ncases=50 > X <- as.matrix(rnorm(N,0,1)) > eta <- -5.3 + X * 1.7 > p <- exp(eta)/(1+exp(eta)) > Y <- rbinom(N,1,p) > controls <- sam...

Hi all, I'm hoping I might be able to get some help with some issues specifying priors for MCMCglmm. I'm trying to fit a gaussian glmm using MCMCglmm to a data set with two (correlated) response variables. The response variables are both logit-transformed proportions (there are a few reasons why I've chosen these with gaussian error over binomal glmm, which I won't go...

I am fairly new at using R/programming in general so I apologize if I am leaving crucial parts of the puzzle out, but here goes. First and foremost this is the error I am receiving: Error in muPriors[priors[, 1:2]] <- priors[, 3] : NAs are not allowed in subscripted assignments This occurs only when I am using observational priors and some number of nominal variables, it does not occur when I change to specifying the variables as ordinal, or if I do not use priors. Using a doctored ve...

Hi Im testing wxRuby. I have a problem when ruby try to load the load wx. Im newbie to ruby, maybe a make a simple mistake : Best regards. Install wxRuby-------------------------------------------------------------- pedro@la-vaca-azul:~$ sudo gem install wxruby2-preview Need to update 2 gems from http://gems.rubyforge.org .. complete Select which gem to install for your platform (i486-linux)

Bugs item #6633, was opened at 2006-11-12 16:46 You can respond by visiting: http://rubyforge.org/tracker/?func=detail&atid=218&aid=6633&group_id=35 Category: Incorrect behavior Group: None Status: Open Resolution: None Priority: 2 Submitted By: Alex Fenton (brokentoy) Assigned to: Kevin Smith (qualitycode) Summary: Serious (but harmless?) GTK warnings on Linux Initial Comment: Whenever any wxruby app is run, GTK throws a lot of warnings which sound serious but appear harmless. This may be an issue with Ubuntu 6.10 GTK...

Dear all, I am trying to use the logistic regression with MCMClogit (package: MCMCpack/Coda) and I want to put a beta prior on the parameters, but it's giving me error message (please see output below) no matter what shape 1 or 2 I use. It works perfect with the cauchy or normal priors. Do you know if there is a catch there somewhere? Thanks logpriorfun <- function(beta,shape1,shape2){ sum(dbeta(beta,shape1,sha...

Hi, everyone, I am a new user of R. Recently, I tried the following commands, but couldn't make them work. If any one of you has some ideas about it, please help me. The commands are >std<-1000 >mean<-8000 >prior<-function(n){1/(sqrt(2*pi)*std)*exp(-1.0*(n-mean)^2/(2*std^2))} >plot(prior,1,15000) >post<-function(n){ + if(n < 9543) + 0 + else + prior(n)/n + } >plot(post,1,15000) # This command didn't work. The error message is " x and y lengths differ " I was really confu...

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.

Hi all, I previously emailed about a multinomial model, and after seeking some additional help, realized that since my response/outcome variables are not mutually exclusive, I need to use a multi-response model that is *not* multinomial. I'm now trying to figure out how to specify the priors on the multi-response model. Any help would be much appreciated. My data look like this: X other focal village present r teaching Opp_teacher Dir_teacher Enh_teacher SocTol_teacher Eval_teacher Total_teacher f_Age f_Ed Age Ed1 61 10202 10213 0 15 0.250000000 2...

I have a problem in where i generate m independent draws from a binomial distribution, say draw1 = rbinom( m , size.a, prob.a ) then I need to use each draw to generate a beta distribution. So, like using a beta prior, binomial likelihood, and obtain beta posterior, m many times. I have not found out a way to vectorize draws from a beta distribution, so I have an explicit for loop within my code for( i in 1: m ) { beta.post = rbeta( 10000, draw1[i] + prior.constant , prior.constant + size.a - dra...

I'm wondering whether anyone has any insight into why the 'simulate' methods for the built-in glm() families (binomial, Poisson, Gamma ...) extract the prior weights using object$prior.weights rather than weights(object,"prior") ? At first I thought this was so that things work correctly when e.g. subset= and na.action=na.exclude are used. However, the current versions of the functions don't seem to work particularly well when these opt...

...t; I previously emailed about a multinomial model, and after seeking some > additional help, realized that since my response/outcome variables are not > mutually exclusive, I need to use a multi-response model that is *not* > multinomial. I'm now trying to figure out how to specify the priors on the > multi-response model. Any help would be much appreciated. > > My data look like this: > > X other focal village present r teaching Opp_teacher > Dir_teacher Enh_teacher SocTol_teacher Eval_teacher Total_teacher > f_Age f_Ed Age Ed1 61 10202 10213...

Hello I wonder whether someone can tell me what I am doing wrong. Here is the code (from Bayesian Computation with R - Chapter 2.3.R) that I am trying to run #################################### # Section 2.3 Using a Discrete Prior #################################### graphics.off() # Close all graphics rm(list=ls()) # Clear all variables library(LearnBayes) library(lattice) p = seq(0.05, 0.95, by = 0.1) prior = c(1, 5.2, 8, 7.2, 4.6, 2.1, 0.7, 0.1, 0, 0) prior = prior/sum(prior) windows() plot(p,...

...g a trial to move from S3 to S4... I have created some classes of interest and they work acceptably well for the purpose. I am now wondering how to make them operate in a package. In clear when a package is loaded (eg library(mypackage)) where should I put the class descriptions: setClass("Prior",representation(Distrib="character",Params="list")) setClass("SamplePrior",representation("Prior",Sample="list")) so that they are created and then usable for functions after a simple call to library(mypackage). It is probably something triv...

Hi, Using the example in ?VSEM in the package BayesianTools I'm attempting to iteratively update the prior but find the plotTimeSeriesResults produces the following errors when I extend the VSEM example in BayesianTools. With the Code below (the errors) I get: " Error in quantile.default(x, probs = quantiles[i]) : missing values and NaN's not allowed if 'na.rm' is FALSE" The e...

...1:9,warn=TRUE) There were 11 warnings (use warnings() to see them) > mc best model: unequal variance with 3 components > mc$mu NULL > sqrt(mc$sigmasq) Error in sqrt(mc$sigmasq) : Non-numeric argument to mathematical function > warnings() Warning messages: 1: In meV(data = data, z = z, prior = prior, control = control, ... : sigma-squared falls below threshold 2: In mclustBIC(data = c(1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, ... : there are missing groups 3: In meV(data = data, z = z, prior = prior, control = control, ... : sigma-squared falls below threshold 4: In mclustBIC(data = c(1...