search for: aolinto

Dear R-List members, I want to check if a set of measurements follows better a gamma or a lognormal distribution (see data below). Using shapiro.test I can test for normality (shapiro.test(log (Lt)). To test for gamma (and normal) distribution I would use ks.test but I need to specify its shape and scale. How should I calculate these values in R? I tried > Lt.fit <- glm(Lt ~ 1,

Hi, I'd like to find which values of x will give me a y. In other words, in the example of a gaussian curve, I want to find the values of x that will give me a density, let's say, of 0.02. curve(((1/(sqrt(2*pi)*10))*exp(-((x-50)^2)/(2*10^2))),xlim=c(0,100)) Thanks for any help, Antonio Olinto ------------------------------------------------- WebMail Bignet - O seu provedor do

...., only one plot with all the points (1 to 10) in blue (or any other color), the points 5 to 9 in red and a regression line from points 5 to 9. Best regards, Antonio Olinto ----- Original Message ----- From: "Rashid Nassar" <rnassar at duke.edu> To: "Antonio Olinto" <aolinto at bignet.com.br> Sent: Sunday, September 23, 2001 2:20 PM Subject: Re: [R] plot and lm > One way to do this (if I understand you correctly): > > > qu > x y > 1 1 1.0 > 2 2 2.0 > 3 3 3.0 > 4 4 4.0 > 5 5 4.4 > 6 6 4.1 > 7 7 3.2 > 8 8...

Dear R-list members, Some months ago I wrote a message on the usage of gnls (nlme library) and here I come again. Let me give an example: I have a 10 year length-at-age data set of 10 fishes (see growth.dat at the end of this message) and I want to fit a von Bertalanffy growth model, Li= Linf*(1-exp(-k*(ti-t0))) where Li = length at age i, Linf= asymptotic length, k= curvature parameter, ti=

Dear list members I want to analyze separately the months of a time series. In other words, I want to plot and fit models for each month separately. Taking the example of http://a-little-book-of-r-for-time-series.readthedocs.org/en/latest/src/timeseries.html births <- scan("http://robjhyndman.com/tsdldata/data/nybirths.dat") birthstimeseries <- ts(births, frequency=12,

Dear list members I have a doubt on how p-values for t-statistics are calculated in the summary of Linear Models. Here goes an example: x <- rnorm(100,50,10) y <- rnorm(100,0,5) fit1<-lm(y~x) summary(fit1) summary(fit1)$coef[2] # b summary(fit1)$coef[4] # Std. Error summary(fit1)$coef[6] # t-statistic summary(fit1)$coef[8] # Pr(>|t| summary(fit1)$df [2] # degrees of freedom #

Dear R-list members Here goes a question on data frames: I want to create a new data frame excluding some records (rows) from an existing one. Considering the data frame "fish.dat" at the end of this message, I can make a boxplot(Ring~Radius) and identify the points at rows 11 and 25 as outliers. With the command fish.dat[c(11, 25),] I can identify which fishes are related to

...Olinto" To: "R-help" <r-help at stat.math.ethz.ch> <aolinto at bignet cc: (bcc: Andrzej P. Jaworski/US-Corporate/3M/US) .com.br> Subject: [R] coefficient of determination on a nls regression...

...summary(fit.lm)$coefficients[,1:2] to get the coefficients and their SEs. For nls, it takes a bit more digging, as the documentation is a bit sparse. Something like: summary(nlsfit)$parameters[,1:2] will give what you want. Cheers, Andy -----Original Message----- From: Antonio Olinto [mailto:aolinto at bignet.com.br] Sent: Wednesday, October 30, 2002 11:05 AM To: R-Help Subject: [R] extracting Std. Error value from lm/nls Hi, How to extract Std. Error values from a lm or a nls? With coef(model), coef(model)[1], coef(model)[2] and df.residual(model) I can get the coeffcients and the degree...

Hello I'm translating a Visual Basic routine to R. Which command in R is similar to "select case" in VB? See the example: select case int(AG) case 0, 5, 8, 10 VAR = 37 case 1, 4 VAR = 70 case 3, 6, 9 VAR = 90 case 2, 7 VAR = 112 end select A chain of "if else" would not be very elegant. Is there another option? Thanks in advance. Best regards. Antonio

Hi, Is it possible to place the group number in a cluster? The best example I found was at https://www.crops.org/images/publications/cs/42/5/1584f4.jpeg (see I, II, IIIa, IIIb) It would help to identify the group given by cutree since branches can turn around and they necessarily are not plotted in order (group 2 is not necessarily the second from the left). Thanks for any help. All the

Hello R-list members, I’m trying to make a min-avg-max plot (like a boxplot, without box and outliers, i.e., for each class of "x" a vertical line indicating the range (max-min) and an horizontal line or point indicating the mean) with data like: Year Avg Min Max 98 10 7 13 99 12 10 15 01 11 6 14 Which command should be used? Many thanks, Antonio Olinto

Hi, Some time ago I asked about how to use comma as decimal separator in plots and Mr. Paul Murrell wrote: ------- You could try something with axis() and chartr(), like ... par(mfrow=c(2,1)) plot(1:10/11, rep(1, 10), main="Standard X-Axis") plot(1:10/11, rep(1, 10), main="Customised X-Axis", axes=F) axis(1, at=pretty(1:10/11),

As to call files of another softwares(as, Excel, minitab) to be manipulated inside of R-plus? Respectfully Regiane -------------- next part -------------- An HTML attachment was scrubbed... URL: https://stat.ethz.ch/pipermail/r-help/attachments/20011112/23b21666/attachment.html

This works fine: fit42<-nls(Vfs~SSlogis(Months,Asym.Int+Asym.Group*Groupdum,xmid,scal), data=df, start=c(Asym.Int=22,Asym.Group=5,xmid=2,scal=6), na.action=na.omit) But this, identical except using gnls, doesn't converge: fit43<-gnls(Vfs~SSlogis(Months,Asym.Int+Asym.Group*Groupdum,xmid,scal), data=df, start=c(Asym.Int=22,Asym.Group=5,xmid=2,scal=6), na.action=na.omit) Error in gnls(Vfs

Hi to all, So I'm also a new user. I downloaded the program last week and I think it's great. Thanks to those who have developed R. I have a special interest in GLM as a tool to describe fisheries and its variables and I'm just begging to study it. As I could understand there's two types of GLM sun of squares: in "type I" the factors are added in sequence and

Dear R-list members, Cynthia M. Jones wrote a paper (Fitting growth curves to retrospective size-at-age data, Fisheries Research 46(2000):123-129; abstract at http://www.elsevier.nl/gej-ng/10/19/44/70/24/37/abstract.html)where the SAS procedure MIXED, Macro NLINMIX (Littell et. al., 1996)was used to estimate the von Bertalanffy growth function parameters assuming that data from the same fish are

Dear R-Users, Is it possible to omit the automatic x labeling in an interaction.plot? xaxt="n" doesn't seem to work. Also I tried to set col.axis="white" but it changes only y annotation color. Thanks, Ant?nio Olinto -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read

Hello I have a large dataset to make a cluster analysis and I'm trying to use Rclusterpp to do so. Nevertheless in a sample matrix analysis (see http://dl.dropbox.com/u/755659/gillnet.txt) I'm getting different results from Rclusterpp and hclust. Following the example from "An Introduction to Rclusterpp" (

Hi, Here goes two plot questions I could not solve by myself: 1) I have a data frame with X (factors) and Y data, ordered by Y values. When I make plot(X,Y) x-axis is rearranged. How to keep the original order? To keep the original order (and to have points plotted, instead of horizontal lines) I tried, for X= names of 8 species and Y= median length: > par(lab = c(8, 5, 7)) >