search for: olinto

...rlay (like add=TRUE) the plots 1 and 2. With plots "at the same screen" I meant overlaid and not one beside the other, i.e., 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): > > > q...

...ite in x axis label "fitted value of lambda" (lambda in greek letter). xlab=expression(lambda) gives the "lambda", I tryed things like xlab=paste ("fitted value of ", expression(lambda)) but I didn't get the greek letter. Thanks in advance for any hint. Antonio Olinto ------------------------------------------------- WebMail Bignet - O seu provedor do litoral www.bignet.com.br

Hello, I want to know how do R calculates the number of intervals between tick-marks in the y axis in a plot. I'm making a three y-axes plot and this information would help me a lot. Thanks in advance. Antonio Olinto ---------------------------------------------------------------- Webmail - iBCMG Internet http://www.ibcmg.com.br

Hello R-Users, I have I simple question that I could not solve: How to extract the elevation and the slope values from a linear model (lm) separately? coef(model.lm) gives both of them. Thanks, Antonio Olinto -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the "body", not the subject !) To: r-help-request at st...

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 litoral www.bignet.com.br

...if. codes: 0 `***' 0.001 `**' 0.01 `*' 0.05 `.' 0.1 ` ' 1 Residual standard error: 10.32 on 97 degrees of freedom Correlation of Parameter Estimates: Linf K K -0.9631 t0 -0.6472 0.7827 Well, I would be grateful to receive any help. Best regards, Antonio Olinto Fisheries Institute Sao Paulo - BRASIL ------------- > growth.dat fish age lt 1 1 1 65.3 2 1 2 144.1 3 1 3 208.6 4 1 4 261.4 5 1 5 304.7 6 1 6 340.1 7 1 7 369.1 8 1 8 392.8 9 1 9 412.2 10 1 10 428.1 11 2 1 68.5...

...rths, frequency=12, start=c(1946,1)) birthstimeseries plot.ts(birthstimeseries) birthstimeseriesHW <- HoltWinters(birthstimeseries) plot(birthstimeseriesHW) How to proceed the plotting and HoltWinters smoothing considereing only Januarys, Februarys, etc. separately. Thanks in advance. Antonio Olinto [[alternative HTML version deleted]]

...Robert H. Shumway Introductory Time Series with R (Use R) by Paul S.P. Cowpertwait and Andrew V. Metcalfe I would like to receive any suggestion of which is most appropriate for a non-statistician (I am a biologist). Reading only the index I could not evaluate it. Thanks in advance, Antonio Olinto Marine Fisheries Center Sao Paulo Fisheries Institute Brazil

...1996)was used to estimate the von Bertalanffy growth function parameters assuming that data from the same fish are repeated-measures with autoregressive correlation structure and fitting with a nonlinear repeated-measures model. Is there any similar procedure that I can use in R? Thanks, Antonio Olinto -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the "body", not the subject !) To: r-help-request at stat....

....975,summary(fit1)$df[2]) # Two-sided p-value should be estimated with dt(t,summary(fit1)$df[2]) # isn't it? But this value is different from summary(fit1)$coef[8] My question is: how to reach to the same p-value indicated in Pr(>|t|) or summary(fit1)$coef[8]? Thanks in advance, Antonio Olinto [[alternative HTML version deleted]]

...Nevertheless, when I ask for the summary (summary(data2)) the others species names still apearing. Also I tried data2<-data1(data1$species %in% c("sp1","sp2")),] But I got the same result. How can I get a "clean" subset? Thanks for any help. Best regards. Antonio Olinto ------------------------------------------------- WebMail Bignet - O seu provedor do litoral www.bignet.com.br

...ist (Year = Y[S=="A" & F=="trawl"]), fun=sum), but it didn't worked. To overcome this problem I did a subset as a new data frame and then I used aggregate. But I'm sure there's a way to apply a filter in data frame. Thanks for any help. Best regards. Antonio Olinto

...elect 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 Olinto [[alternative HTML version deleted]]

...),pos=3,cex=0.5) But the labels are sometimes overlapped. Is there any way to use identify, or a similar command, to plot the row names in a given position? Identify would be perfect but it indicates the row number and I'd like to have the row name given in nmds$points. Many thanks, Antonio Olinto Sao Paulo Fisheries Institute Brasil ------------------------------------------------- BCMG Internet Webmail www.bcmg.com.br

...above each bar. See the example below. I could not place the numbers on the middle of each bar because I could not calculate the right position of each x-axis tick. axTicks(1) indicated a unitary step, but it does not seem work. I appreciate any help or suggestions. Best regards, Antonio Olinto ====================== CAT VAR1 VAR2 Category 01 17.5 9 Category 02 15.2 20 Category 03 10.3 500 Category 04 8.4 150 Category 05 20.3 5000 # Coping data from a spreadsheet dat.data <- read.delim("clipboard",header=T) summary(dat.data) CAT VAR1 VAR2...

...t(x,y) axis(1,at=pretty(x),labels=chartr(".", ",", as.character(pretty(x)))) axis(2,at=pretty(y),labels=chartr(".", ",", as.character(pretty(y)))) I got an "2" in the middle of the "2.0". How could I avoid it? Thanks very much, Antonio Olinto -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the "body", not the subject !) To: r-help-request at sta...

...f GLM sun of squares: in "type I" the factors are added in sequence and "type III" this sequence doesn't matter. In R I always get the message: "Terms added sequentially (first to last)". Is there any way to make the so called Type III procedure? Thanks! Antonio Olinto Fisheries Institute Sao Paulo - Brazil -------------- next part -------------- An HTML attachment was scrubbed... URL: https://stat.ethz.ch/pipermail/r-help/attachments/20001105/3eaf1b64/attachment.html

..."Antonio Olinto" To: "R-help" <r-help at stat.math.ethz.ch> <aolinto at bignet cc: (bcc: Andrzej P. Jaworski/US-Corporate/3M/US)...

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 http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the "body", not the subject !) To: r-help-request at sta...

...-value. so try 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 coeffcient...