search for: plankton

...is also difficult to alter models if I want to change a term. For example, below are the set of models I would like to run. Is there a way to specify the full model and have R generate the rest? I.e. specify m1234567<-glm.convert(glm.nb(mantas~site*year+cosmonth+sinmonth+coslunar+sinlunar+plankton, data=mydata)) and have R run all the other models. library(MASS) #Intercept only m0<-glm.convert(glm.nb(mantas~1,data=mydata)) #One term - 7 models #Manta abundance is greater at one of the two sites m1<-glm.convert(glm.nb(mantas~site,data=mydata)) #Manta abundance increases ea...

...message mean? > Boot2.Pillai <- function(x, ind) { + x <- as.matrix(x[,2:ncol(x)]) + boot.x <- as.factor(x[ind, 1]) + boot.man <- manova(x ~ boot.x) + summary(manova(boot.man))[[4]][[3]] + } > > man.res <- manova(as.matrix(pl.nosite) ~ + as.factor(plankton.new[,1]))$residuals > boot2.plank <- cbind(plankton.new[, 1], man.res) > boot.sep <- boot(boot2.plank, Boot2.Pillai, R = 1000, + strata = plankton.new[, 1]) Error in summary.manova(manova(boot.man)) : residuals have rank 5 < 6 Execution halted A sample of...

Dear R-users -- -- Heather Anne ------------------------------------------- Heather A. Wright, PhD candidate Ecology and Evolution of Plankton Stazione Zoologica Anton Dohrn Villa Comunale 80121 - Napoli, Italy website<http://www.szn.it/SZNWeb/cmd/ShowArchiveItem?TYPE_ID=HUMANRESOURCE&ITEM_ID=21942&LANGUAGE_ID=2> Lab: +39 081 583 3201 Cell: +39 366 484 2196 Fax: +39 081 764 1355 * * [[alternative HTML version deleted]]

...s I obtain in R and excel? Adding the mean to this line instead of abs resulted in obtaining the exact same output dates. Happy to provide more code but this is long enough. Thanks for your input! ------------------------------------------- Heather A. Wright, PhD candidate Ecology and Evolution of Plankton Stazione Zoologica Anton Dohrn Villa Comunale 80121 - Napoli, Italy Lab: +39 081 583 3201 Fax: +39 081 764 1355

...are binary, categorical, and continuous variables.  I ran glm.nb to fit my full model containing all first order terms.  The model I ran was: glm.nb.full<glm.nb(mantas~site+year+cosday+sinday+daylength+lunarpercent+sstmean+sststd+ sshmean+sshstd+cosdir+sindir+spd+temp+alt+tideht+high+falling+low+plankton)   However, when I use stepAIC on the model I get the message:   > stepAIC(glm.nb.full) Start:  AIC=19240.46 mantas ~ site + year + cosday + sinday + daylength + lunarpercent +     sstmean + sststd + sshmean + sshstd + cosdir + sindir + spd +     temp + alt + tideht + high + falling + low + pl...

Hi, Suppose I have: > summary(manova(plank.man)) Df Pillai approx F num Df den Df Pr(>F) plankton.new[, 1] 1 0.5267 9.8316 6 53 2.849e-07 *** Residuals 58 --- Signif. codes: 0 `***' 0.001 `**' 0.01 `*' 0.05 `.' 0.1 ` ' 1 My understanding is the MANOVA summary returns a list. However I am not sure how to ext...

Hello, I am currently trying to show the abundance of two species of zooplankton within the North Sea as pie chart bubble plots. I followed Werner Wernersen's advice in R help (http://finzi.psych.upenn.edu/R/Rhelp02a/archive/48644.html) and used Paul Murrell's paper "Integrating Grid Graphics Output with Base Graphics Output" (in R News) to try and do this (us...

Hello eveybody, I am studying the vertical distribution of plankton and want to study its variations relatively to several factors (time of day, species, water column structure etc.). So my data is special in that, at each sampling site (each observation), I don't have *one* number, I have *several* numbers (abundance of organisms in each depth bin, I s...

The Sieve plugin page of Wiki documentation (http://wiki.dovecot.org/LDA/Sieve?highlight=%28filter%29) says: "It's also not a good idea to build the plugin against self-compiled Dovecot sources, but then actually use a prebuilt binary package of Dovecot." Why isn't a good idea? Ciao, luigi -- / +--[Luigi Rosa]-- \ Every time you've heard someone say they saw a

Hi, I was wondering whether there is an implementation of the Kolmogorov-Smirnov goodness of fit test for discrete, ordinal scale data in R - I've only managed to find the test for continuous data. Thanks! Gila

Dear all I would like to use R and make some maps. I want to have strict control, over the details of the produced map, like remove borders, city names, add markers, add labels. Is there any package apart Rgooglemaps that can do something like that? B.R Alex [[alternative HTML version deleted]]

I know citation() gives the R citation to be used in publications. Has anyone put this into endnote nicely? I'm not very experienced with endnote, and the way I have it at the momeny the 'R Development Core Team' becomes R. D. C. T. etc. Cheers.

Hi all, I have never worked with this kind of data before, so Please help me out with it. I have the following data set, in a csv file, looks like the following: Jan 27, 2010 16:01:24,000 125 - - - Jan 27, 2010 16:06:24,000 125 - - - Jan 27, 2010 16:11:24,000 176 - - - Jan 27, 2010 16:16:25,000 159 - - - Jan 27, 2010 16:21:25,000 142 - - - Jan 27, 2010 16:26:24,000 142 - - - Jan 27, 2010

Dear R-helpers, I have 4-dimensional atmospheric data (x,y,z,t), which I want to analyse on spatio-temporal diversities. As far as I understand there only exists the possibility to construct time series as two-dimensional matrices (mts). For the moment, I hold it in different objects: 1. a four-dimensional array for the spatial related analyses 2. a two-dimensional mts timeserie, which was

...December 2002: an additional month for PASTECS, but not allocated yet between the development of functions, or of the user interface. + January to March 2003: about 30-40% of my time will be dedicated to further develop the GUI in the framework of another project based on R (Zooscan, for automatic plankton recognition). I will thus be able to go further in the development of SciViews Insighter/Xchanger until end of March 2003, but it will certainly not be finished at that time. Since my work is driven by short-term contracts, I have still no plan for the future after March 2003. It leaves some oppor...

Hello, i have a xls file. I will read it in r, what library-command i use for this?? any ideas?? Thanks Michele [[alternative HTML version deleted]]

...F2D-0C0C29C1CE26}' '\\beth' Remoting NetServerEnum to \\beth on transport \Device\NetBT_Tcpip_{D3F23406-7D0D-475B-8F2D-0C0C29C1CE26} with flags ffffffff 3 entries returned. 3 total. 891 milliseconds \\BETH NT 00.00 (W,S,TS,PQ,XN,NT,SS,PBR,MBR,00080000) \\PLANKTON NT 00.00 (W,S,NT,PBR) \\VAV NT 00.00 (W,S,NT,PBR) NOTE: I don't know what the "00080000" is doing! Also, BETH is not listed with the PDC flag, despite the fact that the "getpdc" operation lists BETH. I assume that the reason is that the abo...

...on of Missing Values * wq (0.2-2) Alan Jassby http://crantastic.org/packages/wq Functions to assist in the processing and exploration of data from monitoring programs for aquatic ecosystems. The focus is on time series data for physical and chemical properties of water, as well as the plankton. The package is intended for programs that sample approximately monthly at discrete stations. Updated packages ---------------- actuar (1.1-0), agricolae (1.0-9), AlgDesign (1.1-2), analogue (0.6-23), animation (1.1-0), aroma.affymetrix (1.5.0), aroma.core (1.5.0), arulesSequences (0.1-9), bac...