Thanks a lot Duncan, Jim and Marc!
2017-05-11 22:59 GMT-03:00 Marc Schwartz :
>
> Bingo.
>
> The 'inset' argument is what I was missing. That allows this to be done
> with one step, rather than the two that I had.
>
> Thanks Duncan.
>
> Mar
[[alternative HTML version deleted]]
_
Thanks for you attention Mark. Your approach is really very interesting. I
will try to apply it in my code.
All the best
Antonio
2017-05-11 16:19 GMT-03:00 Marc Schwartz :
>
> > On May 11, 2017, at 1:36 PM, Antonio Silva
> wrote:
> >
> > Hello r-users
> >
&g
Hello r-users
I want to plot some barplots inside a looping with the legend placed
outside the plotting area.
No matter the number of bars I want the legend to be placed centered on the
x-axis.
See the example below
for(i in 1:10) {
var_1 <- sample(1000:10,sample(3:8,1))
ymax <- max(var_1)
*3.5))
>
>
> Bill Dunlap
> TIBCO Software
> wdunlap tibco.com
>
> On Sat, May 6, 2017 at 2:50 PM, Antonio Silva
> wrote:
>
>> Hello
>>
>> I want to make a plot with two "y" axes. The labels at axis 4 should be
>> the
>> values fro
Hello
I want to make a plot with two "y" axes. The labels at axis 4 should be the
values from axis 2 multiplied by a scale factor (as 3.5).
In the example below I draw axis 4 exactly as axis 2. But I could not find
a way to multiply its values by 3.5 (e.g.).
plot(rnorm(100,30,5))
axis(4) # I'd l
Reading the help page with attention I saw that this function fits a
regression model for "circular dependent and linear independent".
So my question now is, is there a way to fit a model been with the circular
as the dependent variable?
Thanks
Antonio
2015-10-24 20:19 GMT-02:
Dear R users
I'm trying to reproduce the results from Lowry et al. 2007 Lunar landings -
Relationship between lunar phase and catch rates for an Australian
gamefish-tournament fisheryFisheries Research 88: 15–23
Basically we have two columns: Lunar days and CPUE (catch per unit effort).
The aim
Thanks Giorgio, David and Rui
With the suggestions my problem was solved in different ways.
Best regards
Antonio
2015-09-28 15:52 GMT-03:00 Giorgio Garziano :
> Try this:
>
> X<-c("A","B","C","D","E")
> Y<-c(0,1,2,3,4)
>
> for (i in 0:3) {
> Y<-Y+i
> data<-data.frame(X,Y)
> fe.flag <- f
Dear R users
I want to write a file that contains several data frames generated in a loop
ing.
I also want the column names be written to file only when it is created in
first loop.
In the example below, when I run each line separately without "for (i in
...) { }" it works, but when I run the lo
Hello,
I have some doubts on the usage of some gdal tools.
After converting a hdf file to tif I want to reproject to SIRGAS2000 and
clip between lats of 22 to 29 S and lats of 40 to 50 W.
HDF file can be downloaded at
https://app.box.com/s/16cf7qv6af6gsz1v66staori2mtneu0r
Basically I'm followin
nice week, best regards
Antonio Olinto
2015-04-04 22:22 GMT-03:00 Jim Lemon :
> Hi Antonio,
> Is it possible to use "add=TRUE" and display the plot in two passes?
>
> Jim
>
>
> On Sat, Apr 4, 2015 at 11:44 PM, Antonio Silva
> wrote:
>
>> Hello everybod
Hello everybody
The problem is that species names are shown in the ordination diagram when
the data set has a maximum of 80 rows, and mine has 81 (!!!).
So now the question is: it is possible to go through this limitation?
Thanks again,
Antônio Olinto
2015-04-03 20:08 GMT-03:00 Antonio Silva
Dear R users
I'm trying to do a RDA analysis based on Borcard et al. 2011 Numerical
Ecology with R examples.
What I cannot understand is why when I run the script (RDA.R) using the
dataset from book (dataset1.rar) RDA triplot shows species names and when I
use my dataset (dataset2.rar) species na
our approach:
> > try:
> >
> > boxplot(mydata[,names(sort(apply(mydata,2,median)))])
> >
> > B.
> >
> >
> > On Mar 21, 2015, at 6:49 PM, Antonio Silva
> wrote:
> >
> >> Thanks Bill and David
> >>
> >> Here goes
factor levels according to a function of the data in each group. E.g.,
> compare the following two plots:
>
>d <- data.frame(Response=cos(1:15), Group=rep(c("A","B","C"),c(6,5,4)))
>par(mfrow=c(1,2))
>boxplot(Response ~ Group, data=d)
>b
Hello
I'm using a dataframe (mydata) where row names are sampling points and
column names are species in a multivariate analysis.
If I write boxplot(mydata) I'll have boxplots for each species abundance in
alphabetical order.
How to get the boxes orderer by the median?
Usually for this I write
Thanks Henrik and Tom, job done. Nice solutions.
All the best
Antonio Olinto
2015-02-02 17:31 GMT-02:00 Henrik Bengtsson :
> On Mon, Feb 2, 2015 at 9:56 AM, Antonio Silva
> wrote:
> > Hi R users
> >
> > I want to save a plot after using the command identify.
> &g
Hi R users
I want to save a plot after using the command identify.
I use identify to place labels manually near the points in order to avoid
overlapping lines and numbers.
x<-c(1,2,3,4,5,6)
y<-c(20,30,15,7,25,40)
plot(x,y,type="b",ylim=c(0,45))
identify(x,y,labels=y)
But when I add
png(filenam
OK, I found the answer:
2*pt(t,summary(fit1)$df[2])
Thanks anyway ...
2013/2/7 Antonio Silva
> 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,
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-sta
t; birthstimeseriesJanFeb1<-ts(birthstimeseriesJanFeb,frequency=2,start=c(1946,1))
> plot.ts(birthstimeseriesJanFeb1)
> birthstimeseriesJanFebHW<-HoltWinters(birthstimeseriesJanFeb1)
> plot(birthstimeseriesJanFebHW)
> A.K.
>
>
>
>
> - Original Message -
>
:wind], start = start(x), frequency = f), seasonal) : time
series has no or less than 2 periods
What could be going wrong?
Best regards
Antonio
2012/12/4 PIKAL Petr
> Hi
>
> > -Original Message-
> > From: r-help-boun...@r-project.org [mailto:r-help-bounces@r-
> &
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.co
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 Rc
ext(cumsum(sapply(r,length)),
rep(mean(tail(unique(cluc$hei),2)), length(r)),
paste(unique(grp[cluc$ord])))
Best regards,
Antonio
2012/3/2 Antonio Silva
> Hi,
>
> Is it possible to place the group number in a cluster?
>
> The best example I found was at
> https://www.
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
>
> require(car)
> ?recode
>
> for an alternative.
>
> HTH,
> Jorge
>
>
> On Wed, Aug 31, 2011 at 1:27 PM, Antonio Silva <> wrote:
>
>> Hello
>>
>> I'm translating a Visual Basic routine to R.
>>
>> Which command in R is si
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 ele
28 matches
Mail list logo
