Search Archives

?Dears, Excuse me. I replied to wrong group. Regards, Mario? [[alternative HTML version deleted]]

Dears, I'm having a weird behaviours when setting arguments in functions. fn <- function(x, st="mean", b=NULL, col.range="black", ...){ dots <- list(...) cat("col.range =", col.range, "\n") cat("dots =\n") print(dots) } fn(1, b=2,col...

Hi all, I have a doubt with multiple regression. I have this model: y ~ x + w + z but the response y to z is humped. I check this by plot and run: lm(y ~ z + I(z^2) # better r...

Hi all, I try to save graphic as eps file, but I have probles with quality. My graphics as pdf file is ok. The quality is very good. But when I save as eps the quality is very poor. I...

Hi all, I am doing anova multi factor and I found different Intercept when model has interaction term. I have the follow data: set.seed(42) dt <- data.frame(f1=c(rep("a",5),rep("b",5)), f2=rep(c...

Hi Jesse, excuse me if you tried it yet, but try this: install.packages("AEM", repos = "http://R-Forge.R-project.org") install.packages("PCNM", repos = "http://R-Forge.R-project.org") or download binary from: http://r-forge.r...

Hi all! I am working with species abundance distribution using biomass as abundance measure. For this, I am trying to fit models gamma, weibull and lognormal to my data. I have problems specially with gamma and weibull. When I try...

Hi all, Bert it's all ok! My disappointments are blame of my expectations and not of R. I started with C and Php and learned to be explicit, 'computer do not to guess what you want', but R do...

Hi Duncan and David, Thank you for explanation. I'm really disappointed with this R "resource". I think that partial match, mainly in function args, must be optional and not default. We can have many problems and lost hours find...

Hi Gian, in your code, you selected column and do logic test when what you want is to do test in line In dune data, sites are lines, then: # Select species with abundance > 0 on line 1 spp1 <- names(dune...

Hi all, I trying develop a code to implement a recommendation of Wilson et al (1991). In they text formula: "... LnA' = fitted mean ln abundance; sigma = fitted standar deviation of ln abundance. In ranked-abundance terms: lnAi = lnA' + sigma * fi...

Thank you Eduard! I understood now! One of the methods of optim use derivative method, but with finite deltax. By default deltax=0.001. Using parscale you chance scale of deltax (parscale*deltax) (page 300 of the book). Best regards...

Hi all, my question is more statistic than r-ecology. I am running a protocol to fit niche apportionment species abundance distribution (SAD) proposed by Muillot et al (2003). This protocol is a modification of Bersier and Sugihara (1997) and...

Hi all, excuse me my question. I read that eps have bmp incorporated and this bmp maybe have different resolution. But eps file not have resolution. I try import in Word and OpenOffice but this editors don't import my...

?Dear Thierry, That is the problem. I read that interpretation is the same, but the Intercept value of summary is different: The mean of level "a" of f1 and level "I" of f2 (first level of each factor) is 0...

Hi Oksanen, thank you for help. I had not thought about estimated parameters. I need search how to do this. Converting fitted values to rank abundance from same models is very difficult. I can't found how to converting fitted...

Hi Manuel, when you use cbind(finaldf[, 3:22], cbind function waiting that you gave second data to bind with first. How you use only one data, its result like: DEV <- finaldf[, 3:22] with the same structure of data...

Dear Klaus, Thank you for your answer and good news that ape/phangorn will load tree with singletons! Can I send the tree to your email? My problem with collapse.singles is, for instance, when I have only one genus...

Dear Eliot, Thank you for your answer! My tree is ok, but to plot is needed to collapse singles, as you said. I use ape/phytools to manipulate my tree for date it after with phylocom. Then, all node is...

Hi Paris, excuse-me if I do not understand you. I think that or you standardize you effort to less individual in your samples or you extrapolation your sample to that more abundant. If you have sample with less individual...