search for: peces

Estimados miembros de la lista. Estoy usando 2 funciones del paquete vcd y obtengo resultados que parecen discrepar entre ambas funciones. > dput(peces) c(26L, 63L, 1L, 21L, 20L, 50L, 0L, 104L, 19L, 0L, 8L, 32L, 0L, 24L, 5L, 2L, 0L, 36L, 0L, 5L, 27L, 1L, 13L, 24L, 15L, 25L, 23L, 5L, 37L, 26L, 29L, 14L, 15L) La función Ord_plot me dice que la mejor distribución para los datos es binomial negativa. Ord_plot(peces) La función distplot parece indic...

Hola a tod en s, Trabajo con datos de abundancia (número de individuos) por tallas (medidas al medio centímetro inferior) de diferentes especies de peces y mis datos base son el número medido de peces por talla. Necesitaría vuestra ayuda para transformar las distribuciones de frecuencias de talla a los datos brutos que las generaían. Es decir, pasar de algo como esto 8           1 8.5        3 9          5 10        7 10.5     4 11        2 11...

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,

...amlss para este caso es "beta inflated" (Stasinopulos and Rigby 2007. Journal of Statistical Software 23(7)). El modelo que intento correr es: m1<-gamlss(Teleosteos ~ cs(LT_max)+random(Familia/Genero,lambda=T),family=BEINF, data = datos) donde LT_max es la talla máxima de especies de peces y la variable respuesta es la proporción del consumo de peces teleósteos de cada especie. El efecto aleatorio es anidado de géneros y familias taxonómicas de especies (género dentro de familia). El problema que tengo es que los gamlss no toman el efecto aleatorio como random(Familia/Genero) o (1|Fa...

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

Hola, buen día. Una consulta, no sé conozco en extenso todas las bondades de R, pero quería saber exactamente si hay algún paquete que permita realizar análisis ecológicos específicamente para comunidades marinas. Existen muchos índices, indicadores, estimadores en general en muchos programas que están a nivel de usuario, pero quería saber si todos estos pueden estar ya disponibles en R.

Forgot to copy the list Dear Jim, Thank you for your reply. I must be doing something wrong, If this is my command to plot a map of Europe: eps_europe <- mapCountryData(n, nameColumnToPlot="eps_score", mapTitle="EPS Score - Europe",colourPalette=colourPalette, catMethod="fixedWidth", missingCountryCol = "white",

Dear David, Thank you for your answer. Sorry for the embarrassing mistake. However, even with when I generate a map for the whole world using: eps <- mapCountryData(n, nameColumnToPlot="eps_score", mapTitle="EPS Score",colourPalette=colourPalette, catMethod="fixedWidth", missingCountryCol = "white", addLegend=FALSE) And then use:

> On Apr 9, 2016, at 8:13 AM, Miluji Sb <milujisb at gmail.com> wrote: > > Forgot to copy the list > > Dear Jim, > > Thank you for your reply. I must be doing something wrong, If this is my > command to plot a map of Europe: > > eps_europe <- mapCountryData(n, nameColumnToPlot="eps_score", mapTitle="EPS > Score -

Amigos de la erre. He creado mi primer bucle con for para entrenar unos modelos con GAM. La respuesta es quasipoisson porque estoy trabajando con densidades de peces. Sin embargo, tengo un problema, no se muy bien como añadir el offset a la formula siguiente cuando creo el bucle. GAM.A1 <-gam ((DYO)~s(DMA,k=4)+ s(WOD,k=4)+s(CIN,k=4)+s(DRA,k=4)+s(DBR,k=4) +offset(log(ARE)),family=quasipoisson(link = "log"),data=as.data.frame(cbind(DYO,Conjunto)))...

> On Apr 9, 2016, at 11:18 AM, David Winsemius <dwinsemius at comcast.net> wrote: > > >> On Apr 9, 2016, at 10:46 AM, Miluji Sb <milujisb at gmail.com> wrote: >> >> Dear David, >> >> Thank you for your answer. Sorry for the embarrassing mistake. However, even with when I generate a map for the whole world using: >> >> eps <-

...n la forma. Se me está pasando algo por alto y no doy con ello ¿Alguna idea? Un saludo y Gracias Juan Carlos *************************************** Si las cuerdas del trapecio las corté fue para subir y atarlas y ver a la luna otra vez. Y volverlas a cortar una y mil veces y boquear como los peces cuando les toca perder. [[alternative HTML version deleted]]

Hello David, This is exactly what I want but I still can't get the arrows. R and R studio is updated. Thanks again! Sincerely, Milu On Sat, Apr 9, 2016 at 10:29 PM, David Winsemius <dwinsemius at comcast.net> wrote: > > > On Apr 9, 2016, at 1:27 PM, David Winsemius <dwinsemius at comcast.net> > wrote: > > > > > >> On Apr 9, 2016, at 11:18 AM,

...evels) of a bird presence/absence (pres.f) on 14 artificial ponds (charca.f). The ponds are integrated in 4 greater spatial units (zepa). I have 11 explanatory variables, some of them are continuous (prof1,prof2,prof4,capac,perim_veget,div_veg,div_herp) and some are factors (ungul.f,ganado.f,turb.f,peces.f). Not all artificial ponds have been monitored in all seasons, so I have some missing values in the dependent variable. I have thought that the correct model should be, ignoring the spatial pseudorreplication: m1<-lmer(pres.f~ganado.f+ungul.f+turb.f+prof1+prof2+prof4+capac+perim_veget+pec...

> On Apr 9, 2016, at 10:46 AM, Miluji Sb <milujisb at gmail.com> wrote: > > Dear David, > > Thank you for your answer. Sorry for the embarrassing mistake. However, even with when I generate a map for the whole world using: > > eps <- mapCountryData(n, nameColumnToPlot="eps_score", mapTitle="EPS Score",colourPalette=colourPalette, >

Dear David, The device was the issue. The quartz() device works fine but pdf() does not. Now I just need to figure out the limits for map for Europe. Thanks for all your help and patience. Sincerely, Milu On Sun, Apr 10, 2016 at 7:10 PM, David Winsemius <dwinsemius at comcast.net> wrote: > > > On Apr 10, 2016, at 4:12 AM, Miluji Sb <milujisb at gmail.com> wrote: > >

> On Apr 10, 2016, at 4:12 AM, Miluji Sb <milujisb at gmail.com> wrote: > > Hello David, > > This is exactly what I want but I still can't get the arrows. R and R studio is updated. Thanks again! I didn't try it in Rstudio until just now (and I don't remember that you ever mentioned RStudio as a possible issue.) The plotting I see in the default graphics Rstudio

> On Apr 9, 2016, at 1:27 PM, David Winsemius <dwinsemius at comcast.net> wrote: > > >> On Apr 9, 2016, at 11:18 AM, David Winsemius <dwinsemius at comcast.net> wrote: >> >> >>> On Apr 9, 2016, at 10:46 AM, Miluji Sb <milujisb at gmail.com> wrote: >>> >>> Dear David, >>> >>> Thank you for your answer.

> On Apr 10, 2016, at 1:45 PM, Miluji Sb <milujisb at gmail.com> wrote: > > Dear David, > > The device was the issue. The quartz() device works fine but pdf() does not. Now I just need to figure out the limits for map for Europe. Thanks for all your help and patience. After plotting a map of Europe with base graphics the coordinates of the lower-left and upper-right

Dear David, Thank you very much for your replies! I didn't know about par('usr'). I get different coordinates though: [1] -19.75966 54.75966 33.60000 71.40000 But the arrow is not at the bottom of the map. I will keep playing with this. Thanks again! Sincerely, Milu On Mon, Apr 11, 2016 at 12:00 AM, David Winsemius <dwinsemius at comcast.net> wrote: > > > >