similar to: Form of the equation produced by a GLM with Poisson family and log link function

Dear R experts, I have a Matlab code which I am translating to R in order to examine and enhance it. First of all, I need to reproduce in R the results which were already obtained in Matlab (to make sure that everything is correct). There are some matrix manipulations and '\' operation among them in the code. I have the following data frame > ABS.df Pro syn

In the HistData package, I have a data frame, PearsonLee, containing observations on heights of parent and child, in weighted form: library(HistData) > str(PearsonLee) 'data.frame': 746 obs. of 6 variables: $ child : num 59.5 59.5 59.5 60.5 60.5 61.5 61.5 61.5 61.5 61.5 ... $ parent : num 62.5 63.5 64.5 62.5 66.5 59.5 60.5 62.5 63.5 64.5 ... $ frequency: num 0.5 0.5

Hello I would like to pass a model formula as an argument to the with.mids function from the mice package. The with.mids functon fits models to multiply imputed data sets. Here's a simple example library(mice) #Create multiple imputations on the nhanes data contained in the mice package. imp <- mice(nahnes) #Fitting a linear model with each imputed data set the regular way works

Dear Experts, I have a data object named "data.char". When I use cbind(data.char) in SPlus, I got the following results: > cbind(data.char) data.char data matrix, 2700 excluded list, 3 cluster.var character, 2 strata list, 3 xlog CT link logit gpcorr 1 missing.row numeric, 0 zero.row numeric, 22 infile perc.csv attr(,

Estimados, Les tengo una pequeña consulta. Estoy trabajando con ASReml en Rstudio, me encuentro corriendo un modelo (CHL) para obtener el ANOVA de la interacción de los factores gen x env CHL<- asreml(fixed= MS~geno:env, random = ~rep, data = index) (interacción) Sin embargo, cuando ejecuto el modelo, R indica que siguiente mensaje: Current workspace: 128.000000Mb Warning message: Abnormal

Hi all. I'm new to dovecot. I promise I've looked over the docs and lists for an answer to this, but I didn't find anything. If I missed something I apologize. I'm running dovecot 1.0.beta2 on debian testing. It's almost working, I'm getting permissions errors when the imap server actually goes to read the messages out of the cur Maildir directory. The logs spit these

Hi! OK, I'm trying to select some "useful outliers" from my dataset: I defined 11 "treshold" values (1 for each level of a variable (sampling site) as follows: tresholds<-function(x) { tapply(x,mm$NAME,FUN=mean ,simplify = T, na.rm=T)->med tapply(x,mm$NAME,FUN=sd ,simplify = T, na.rm=T)->standev standev+med } tresholds(mm$chl) Now I'd like to select

Dear R users, I am having a problem with the main title of a plot. I have written a function returning a plot and I would like that the main title of this plot shows actually the name of the dataset I have used. So, for instance: if my dataset is called Chl, I would like that my plot to be entitled : "Raw periodogram of Chl". I know the function "paste". I suppose I would be

I have the following data frame: > algaedata = year DIV cellsperml 2001 BAC 72.808 2001 CHL 3.273 2002 BAC 14.002 2002 CYA 220.896 2002 UNI 464.699 2003 BAP 0 2003 BAC 1.782 2004 CYA 315.799 2005 UNI 222.532 2005 BAP 0.2 2005 CYA 163.627 2005 BAC 324.949 2006 CHL 1.636 2006 BAC 199.145 2007 CHL 19.635 2007 CYA 134.174 2007 BAC 485.405 2007 CHL 11.454 2007 CYA 104.721 ...which makes a fine

I'm trying to run "rfe" for variable selection in the caret package, and am getting an error. My data frame includes a dummy variable with 3 levels. x <- chlDescr y <- chl #crate dummy variable levels(x$State) <- c("AL","GA","FL") dummy <- model.matrix(~State,x) z <- cbind(dummy, x) #remove State category variable w <- z[,c(-4)]

Hi guys, the SUPER codec by Erightsoft http://www.erightsoft.net/SUPER.html contains lots of GPL and LGPL code: mplayer, ffmpeg, x264, musepack, theora, which they admit and give credit for. Still, their product is proprietary, and they insist on it. I tried to get the source through their forum, but they of course won't give it: http://www.erightsoft.net/Supforum.html I'll forward

Hi all. I got my previous problem fixed. Effective uid and gid fix in dovecot.conf file made it work. However, I'm still seeing this in my /var/log/dovecot.conf dovecot: 2006-04-13 21:09:00 Info: imap(chl at iaminla.com): Effective uid=1000, gid=1000 The uid and gid *should* both be 1009, and I can't for the life of me figure out how to configure that. I looked around the list and

Hello I have troubles running TextAloundMp3 under wine. TextAloudMP3 version is 1.459 Wine version is 20050725 My TextAloudMp3 directory looks like: $ ls TextAloudMp3 -rw------- 1 atom users 842 Sep 9 14:31 27bluL2.gif -rw------- 1 atom users 842 Sep 9 14:31 27bluR2.gif -rw------- 1 atom users 842 Sep 9 14:31 27blulft.gif -rw------- 1 atom users 845 Sep 9 14:31 27blurgt.gif

Hi All, I am trying to reproduce this using R instead. [image: Full-size image (38 K)] I tried using the following code *SChla <- read.csv("SM_Chla_data.csv")* *Atlantis <- SChla[16:66,]* *head(Atlantis)* * * Seamount Station Depth Pico Nano Micro Total_Ch dbar Latitude Longitud 16 Atlantis 1217 Surface 0.0639 0.1560 0.0398 0.2597 2.082 -32.71450 57.29733

I am trying to draw maps for the world using: library(rworldmap) library(maptools) library(RColorBrewer) tmp2<- dput(head(pece,10)) structure(list(iso3 = c("AUS", "AUT", "BEL", "CAN", "CHE", "CHL", "CZE", "DEU", "DNK", "ESP"), eps_score = c(0.877343773841858, 2.68984365463257,

Hi All, I'm using the following code to produce some stacked bar graphs. *setwd("C:\\Users\\Tinus\\Documents\\NMMU\\R\\Seamounts")* *SChla <- read.csv("SM_Chla_data.csv")* * * *#Extract mean values from data file* * * *Coral <- SChla[185:223,] #Reduce SChla to Coral only* *coral <- with(Coral , aggregate(cbind(Pico, Nano, Micro), list(Depth),FUN=mean))*

Hi, I am looking at the effects of two explanatory variables on chlorophyll. The data are an annual time-series (so are autocorrelated) and the relationships are non-linear. I want to account for autocorrelation in my model. The model I am trying to use is this: Library(mgcv) gam1 <-gam(Chl~s(wintersecchi)+s(SST),family=gaussian, na.action=na.omit, correlation=corAR1(form =~

I am writing a package, which contains several functions and variables. The variables are for internal use by functions. Functions are ment to be callable by a user. However, the function cannot access these variables. The package contains only R code, and was created using package.skeleton package.skeleton(name="ac9", list=c(ls(pattern="(AB2)|(ac9)"),

Dear All, Please see the short script at the end of the email, which I assembled looking for bits and pieces on the web. It essentially does what I need: it plots several countries as a color-coded map. I just would like to fine-tune a bit the final image, in particular (1) Select my own colors for "painting" the countries (i.e. associate manually a color to every level) (2) Be

Hi Miluji, Try this: arrows(-100,-140,100,-140,code=3) Jim On Fri, Apr 8, 2016 at 10:24 PM, Miluji Sb <milujisb at gmail.com> wrote: > I am trying to draw maps for the world using: > > library(rworldmap) > library(maptools) > library(RColorBrewer) > > > tmp2<- dput(head(pece,10)) > structure(list(iso3 = c("AUS", "AUT", "BEL",