similar to: Error: could not find function "invlogit" and "bayesglm"

Dear R users: I have run a logistic regression, used Gelman et al.'s car package to simulate the parameter estimates of that model, and have plotted the probability (using Gelman et al.'s invlogit() function) of the dependent variable being 1 given the value of a particular independent variable is at its mean. The plot has probabilities on the y-axis and the number (1-1000) of the

-----BEGIN PGP SIGNED MESSAGE----- Hash: SHA1 While I'm not an R expert, I have used R on Windows XP. Now I've moved to Ubuntu (Intrepid), and I'm trying to configure R to work with the Gelman and Hill _Data Analysis Using Regression and Multilevel/Hierarchical Models_. So far, it's not working. I start by following the instructions for installing arm and BRugs at

Dear list, I have troubles with using Bayesian logistic model with count data in bayesglm. If I consider the following artificial data set, with a response y and a covariate x, in a form of row data (a) and count (b) : a<-data.frame(y=c(rep(1,10),rep(0,6)),x=c(rep(5,6),rep(4,4),rep(5,1),rep(4,5))) a$un<-1 b<-aggregate(a$un,list(y=a$y,x=a$x),sum) names(b)[3]<-"w" A

Hi R bayesians, I need an advise how to resolve the two different estimates applying a traditional glm (TG) and a bayes glm (BG), and different results depending on the data formats of response data and the prior specs using bayesglm in R. I'm not familiar with bayes estimate and my colleague asked me to look into this because the EPA from France reported a quite different estimates for

Hi folks, Warning: I don't know if the result I am getting makes sense, so this may be a statistics question. The fitted values from my binomial lmer mixed model seem to consistently overestimate the cell means, and I don't know why. I assume I am doing something stupid. Below I include code, and a binary image of the data is available at this link:

[This email is either empty or too large to be displayed at this time]

Hello, I am a bit over my head on this issue. My colleagues and I are running R off of our server. We all have admin rights and prior to yesterday we all had our own libraries. Our main system administrator advised us that we should have a shared library. So, I am trying to do this. I have downloaded the latest version of R and installed it on the main drive of our server in the "Program

Dear everyone, I'm encountering an unknown error message when using gamm function: > fitoutput <- gamm(cvd~as.factor(dow)+pm10+s(time,bs="cr",k=15,fx=TRUE)+s(tmean,bs="cr",k=7,fx=TRUE) + ,correlation=corAR1(form=~1|city),family=poisson,random=list(city=~pm10),data=mimp) Maximum number of PQL iterations: 20 iteration 1 iteration 2 iteration 3 iteration 4

Dear List, I am making a list of all the files in a folder and I want to exclude one file called "proj.current". This is the 31st file in a folder. Currently I use the command: tdirs <- list.files(pattern="proj.")[-31] However I would like to exclude it based on its name and not the position. Any advice would be useful. Thanks, Daisy Daisy Englert Duursma Room

Hello help, I have changed around some graphing code and made it into a function. Previously they y label of the axis was inserted as text in its own layout box. text(1,1, expression(~degree~C),cex=1) This worked great and resulted in the symbol for degree. In the function, I have changed it so: text(1,1,paste(b_unit),cex=1) and b_unit<-expression(~degree~C) This now inserts ~degree~C

Dear all, I am encountering some issues with my data and need some help. I am trying to run glm analysis with a presence/absence variable as response variable and several explanatory variable (time, location, presence/absence data, abundance data). First I tried to use the glm() function, however I was having 2 warnings concerning glm.fit () : # 1: glm.fit: algorithm did not converge # 2:

Hola, Deseo simular un modelo logit con interacción, estimar sus coeficientes y comprobar si son o no parecidos al modelo teórico. Con este ejemplo obtengo que los coeficientes estimados no se asemejan mucho a los originales. ¿Se le ocurre a alguien cuál es el motivo de esta discrepancia? ¿y cómo solucionarlo? Muchas gracias Emilio logisticsimulation <- function(n){ dat <-

Hi, winetools 2.1.1 (installed from wine.sf.net deb recently) has a few missing double quotes on line 2568. Here's the diff: tjzeeman@carrington:~$ diff /chroot/ia32-sid/root/winetools.orig /chroot/ia32-sid/usr/bin/winetools 2568c2568 < if [ $WINEVER -lt 20040914 ] || [ $WINEVER -gt 20050111 ]; then --- > if [ "$WINEVER" -lt "20040914" ] || [ "$WINEVER"

Hello, We are comparing some features of our product against a competitor. Since the product is produced in lots and we have proof that ambient temperature is relevant for its functioning we used a mixed effect model (MEM) structured this way (NB: R lme4 notation): `out ~ Brand * Temperature + (Temperature | Lot)` with out being the various outcomes, a random intercept on Lot and a random slope

Hi All, I've been having a little trouble with creating a loop that will run a a series of t.tests for inspection, Below is the code i've tried, and some checks i've looked at. I've used the get(paste()) idea as i was told previously that the use of the eval should try and be avoided. I've run a single syntax to check that my systax is correct and works without any problems

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: