search for: chl

...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 unk Chl Y Nh V1 1 0.056524968 0.04387755 -0.073925372 0.026477146 0.083527021 0.0031950622 0.02164793 2 0.066456284 0.05992579 -0.094609497 0.031772575 0.054881164 0.0022078221 0.01358594 3 0.036383887 0.04601113 -0.061213302 0.021181717 0.026713530 0.00117...

...5 59.5 60.5 62.5 63.5 64.5 ... $ frequency: num 0.5 0.5 1 0.5 1 0.25 0.25 0.5 1 0.25 ... $ gp : Factor w/ 4 levels "fd","fs","md",..: 2 2 2 2 2 2 2 2 2 2 ... $ par : Factor w/ 2 levels "Father","Mother": 1 1 1 1 1 1 1 1 1 1 ... $ chl : Factor w/ 2 levels "Daughter","Son": 2 2 2 2 2 2 2 2 2 2 ... I want to make a 2x2 set of plots of child ~ parent | par+chl, with regression lines and loess smooths, that incorporate weights=frequency. The "frequencies" are not integers, so I can't simply...

...uncton 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 fine with(imp, lm(bmi~hyp+chl)) #Creating a formula object and than passing it as argument does not work: form.obj <- formula(bmi~hyp+chl) with(imp, lm(form.obj)) #The following doesn't work either expr <-"lm(bmi~hyp+chl)" with(imp, expr) Looking at the definition of with.mids reveals that the second argu...

...GpSize Incid SevLo SevHi [1,] "ACC-2688" "1" "3" "HU" "r" "M" "686" "7" "12" "C" "n" "12" "1" "1" "1" [2,] "CHL-0259" "1" "2" "MU" "OT" "F" "1356.58" "4" "10" "L" "n" "1" "4" "0" "0" [3,] "CHL-0259" "1" "2&quot...

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 termination Insuffici...

...; 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 those values from vector mm$chl that are higher than each "treshold value", but how can I compare a vector with 1885 elements with the one with 11? Sorry for this (probably) stupid question... and thanks in advance. Alessandro

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 looking for something like: plot(main=paste(""Raw periodogram of", "...). How can the program understand that the name of the data set is s...

...ssed 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 out: dovecot: 2006-04-12 22:32:39 Error: imap(chl at iaminla.com): open(/var/vmail/iaminla.com/chl/cur/1144905139.V306I31c02bM482275.chldebian.chrislea.com:2,) failed: Permission denied The ls of that directory looks like this: chldebian:/etc/dovecot# ls -al /var/vmail/iaminla.com/chl/cur/ total 12 drwxrwsr-x 2 vmail vmail 4096 2006-04-12 22:12...

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 fin...

...I want to know how to incorporate the Poisson family and log link function in the equation. I would really appreciate if someone could help me with this or direct me towards literature that would help me do this. Here are the coefficients from my model: Call: glm(formula = count ~ slope + sst + chl + dist + ice, family = poisson, data = dfo2003cc) Coefficients: (Intercept) slope sst chl dist ice -6.884e-01 -6.740e-01 5.644e-01 -2.694e+00 -7.278e-05 1.044e-01 Degrees of Freedom: 412 Total (i.e. Null); 407 Residual Null Deviance...

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)] subsets <- c(2:8) ctrl<- rfeControl(functions = lmFuncs, method="cv&...

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 couldn't find the answer. Help a noob out please? I promise I'll write up a detailed howto for my debian +...

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

...2.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 -rw------- 1 atom users 306688 Sep 9 14:31 AppCloserProject.exe -rw------- 1 atom users 7891 Sep 9 14:31 AudioHelp.chl -rw------- 1 atom users 493 Sep 9 14:31 AudioTest.chl -rw------- 1 atom users 4050 Sep 9 14:31 FAQ2Grey2.jpg drwx------ 2 atom users 4096 Sep 9 14:31 Files -rw------- 1 atom users 66560 Sep 9 14:31 KEYPUSH.DLL -rw------- 1 atom users 2498 Sep 9 14:31 KeyMapping.chl -rw------- 1...

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 17 Atlantis 1217 Shallow 0...

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, 1.31406247615814, 1.98046875, 2.61666655540466, NA, 1.44414067268372, 2.34257817268372, 2.89687490463257, 2.15937495231628 ), gov_eff = c(1.76499999562899, 1.85666667421659...

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))* *rownames(coral) <- coral[,1] * *coral <- t(coral[,-1])...

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,...

...ctions 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)"), "Chl","Pro","LP","SP","syn","unk","Y","Nh", "plot.C","read.AC9.data","EstimC","plot.AC9"), namespace=TRUE,force=TRUE) Among others, the package co...

...eFile <- loadGADM(fileName, level) theFile <- changeGADMPrefix(theFile, fileName) theFile } ## this function creates a SpatialPolygonsDataFrame that contains all maps you specify in "fileNames". ## E.g.: ## spdf <- getCountries(c("ARG","BOL","CHL")) ## plot(spdf) # should draw a map with Brasil, Argentina and Chile on it. getCountries <- function (fileNames, level = 0, ...) { polygon <- sapply(fileNames, loadChangePrefix, level) polyMap <- do.call("rbind", polygon) polyMap } spdf <- getCountries(c(...