3
B
B04_R1
0.974706
4
B
B04_R2
0.966461
4
B
B04_R3
0.949697
4
B
B05_R1
0.958339
5
B
B05_R2
0.952168
5
B
B05_R3
0.96301
5
B
Thank you very much.
Regards.
---
MSc. <mailto:r.alui...@gmail.com> R
Hi, a correction in the previous script. Anyone who tries to use it must
ignore (remove) the [,-7] after the Grupos object.
Sorry.
De: Rodrigo Aluizio [mailto:r.alui...@gmail.com]
Enviada em: terça-feira, 11 de janeiro de 2011 11:01
Para: R Help
Assunto: Postscript function Bug at R x64
rsion nodename machine
login
"Windows" "7 x64" "build 7600" "RODRIGO-PC" "x86-64"
"Rodrigo"
user
"Rodrigo"
Hope you can help me find out what is going on!
Regards
-
Well, at least I could help (even if unintentionally), as a
compensation for all the help I got from the list.
Cheers.
Em sexta-feira, 10 de dezembro de 2010, Jim Lemon escreveu:
> On 12/10/2010 08:48 PM, Rodrigo Aluizio wrote:
>
> OK that's it. Working nicely. I sent the final
gt;> Point RainSalt
>> Fev/03 365.6 13
>> Mar/03 235 18
>> Abr/03 115.1 18
>> Mai/03 47.4 18.75
>> Jun/03 112 15
>> Jul/03 156.8 17
>> Ago/03 66.115
>> Set/03 149.8 14
>> Out/03 167,1 11.5
>> Nov/03 269.
-feira, 9 de dezembro de 2010 17:57
Para: Rodrigo Aluizio; 'R Help'
Assunto: RE: [R] Barplot with "Independent" Lines Y axis
Without seeing a reproducible example we cannot be sure, but my guess is
that you are letting twoord.plot set the limits and the function does not
remove mi
this time
it's not a choice of mine (unfortunately). Thank you anyway for the
highlight, maybe it will help me arguing.
Regards
Rodrigo.
-Mensagem original-
De: Greg Snow [mailto:greg.s...@imail.org]
Enviada em: quinta-feira, 9 de dezembro de 2010 15:56
Para: Rodrigo Aluizio; R Help
p=c(0,370,10))
lines(mp,Dados$Salt,type='b',pch=19,lty='dotted')
Any ideas?
Thank you for the attention.
Regards
-----------
MSc. Rodrigo Aluizio
Centro de Estudos do Mar/UFPR
___
on and sorry for the silly question.
---
MSc. <mailto:r.alui...@gmail.com> Rodrigo Aluizio
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posti
my system components (drivers)?
Thank you very much for the attention.
Rodrigo Aluizio
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http:
I'm not sure if this is what you want. But the function coordinates() in sp
package gives you the coordinates of SpatialObjects.
Regards.
Rodrigo.
2010/6/9 Nikhil Kaza
> You need to execute gpclibPermit() to enable gpclib.
>
> library(maptools) should have issued a warning to that effect.
>
>
Well, it's may not be the best way to do so, but you may use this:
colnames(foo[n])
where n is the column number, or the column name inside '' ''.
The $ can't be used because it turn the column content into a vector,
loosing its "properties" (e.g.: column and row names). But when you using
[], t
Hi Martin, here as example code that may help you to get what you want (once
you didn't specified were in the plot you want the symbol).
plot(1:10)
title(main= "\u20ac", font = 5)
Ps.: This coding may differ between operational system, this one worked in a
Windows Vista 32 bits.
This works for al
Or maybe another output image format, such as tiff or png?
2010/2/10 Uwe Ligges
>
>
> On 10.02.2010 13:16, Jason Rupert wrote:
>
>> Hello R-Help,
>>
>> I've got a bit of an issue with WMF's. I am working on WindowsXP and
>> outputting WMF format images. I then take the WMF format images and
insert the coordinates that will close the BP[1] line into
this object without changing the object properties. Any Ideas?
Thank you in advance.
-
MSc. Rodrigo Aluizio <mailto:r.alui...@gmail.com>
Centro de Estudos do Ma
R2HTML_1.59-1 Hmisc_3.7-0
[5] survival_2.35-8
loaded via a namespace (and not attached):
[1] cluster_1.12.1 grid_2.10.1 lattice_0.17-26 svMisc_0.9-56
[5] tools_2.10.1
-
MSc. <mailto:r.alui...@gmail.com> Rodri
s probably the Cochran Q Test.
So, any ideas?
Thank you for the pacience.
-
MSc. <mailto:r.alui...@gmail.com> Rodrigo Aluizio
Centro de Estudos do Mar/UFPR
Laboratório de Micropaleontologia
[[alternative HTML v
,lwd=2,ci.lty=0,ci.type=
'poly',ci.col='lightgray',ylab='Riqueza',ylim=c(0,250),yaxp=c(0,250,10),bty=
'n',main='Fauna Viva - Bacia de Campos - Março 2009')
Thank you very much.
Regards.
-----------
,sep=''),c('Days','Hum','Temp')),byrow=F))
Data$Hum<-as.numeric(Data$Hum)
Data$Temp<-as.numeric(Data$Temp)
Then you can plot it.
-----Mensagem original-
De: joris meys [mailto:jorism...@gmail.com]
Enviada em: terça-feira, 6 de outubro de 2009 12:08
This may work for you.
Then you can custom your graphics with ?par.
Data<-as.data.frame(matrix(c("Monday", "Tuesday", "Wed", "Thu", "Fri",
"Sat","Sun",56,57,60,75,62,67,70,76,77,81,95,82,77,83),ncol=3,nrow=7,list(pa
ste('Row',1:7,sep=''),c('Days','Hum','Temp')),byrow=F))
plot(Data$Days,Data$Temp)
Well, on Linux => Emacs+ Ess
On Windows => Tinn-R
-
MSc. Rodrigo Aluizio
Centro de Estudos do Mar/UFPR
Laboratório de Micropaleontologia
__
R-help@r-project.org mailing list
-------
MSc. <mailto:r.alui...@gmail.com> Rodrigo Aluizio
Centro de Estudos do Mar/UFPR
Laboratório de Micropaleontologia
Avenida Beira Mar s/n - CEP 83255-000
Pontal do Paraná - PR - Brasil
[[alternative HTML version deleted]]
_
I'm not shure but I guess that you miss a " and putted it in the wrong
place!
Try this: postscript("/pathto/filename.eps",horizontal=FALSE,onefile=FALSE)
Hope It helps.
-Mensagem original-
De: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] Em
nome de e-letter
Enviada e
Hi, take a look at the gdata package (any plataform) function read.xls(),
and at the package RODBC (if you are on a Windows Plataform) at the
functions odbcConnectExcel(), odbcConnectExcel2007() and sqlFetch().
Hope this helps
Rodrigo.
2009/8/11 rajclinasia
>
> Hi Every one,
> I have a problem
You can save plots as specified below. For details see ?png, ?postscrip,
?tiff, ?jpeg, ?bmp.
postscript(file = ifelse(onefile, "Rplots.ps", "Rplot%03d.ps"),
onefile, family, title, fonts, encoding, bg, fg,
width, height, horizontal, pointsize,
paper, pagecentre, pr
Hi list,
By the entire last 2 weeks I was looking for a way to directly import xlsx
files to R in a Linux OS (Ubuntu 9.04). I already read the R Import/Export
guide, and I know how to use gdata to import xls files and read.table to
import .csv. My problem is that all data that I receive is in the x
names=F)
mo
Hope it helps.
---------
MSc. Rodrigo Aluizio
Centro de Estudos do Mar/UFPR
Laboratório de Micropaleontologia
Avenida Beira Mar s/n - CEP 83255-000
Pontal do Paraná - PR - Brasil
-Mensagem original-
De: r-help-boun...@r-project.org [mailto:r-help-boun...@r-pr
Thank you so much Gilles, that is exactly what Im looking for.
Rodrigo.
De: LE PAPE Gilles [mailto:lepape.gil...@neuf.fr]
Enviada em: sexta-feira, 19 de junho de 2009 10:14
Para: r.alui...@gmail.com
Assunto: Post-hoc test
Hi Rodrigo,
you can use the nparcomp function, library(nparcom
to define a way to it.
Do You have any suggestion? Other equivalent test maybe!
Thanks in advance.
-
MSc. <mailto:r.alui...@gmail.com> Rodrigo Aluizio
Centro de Estudos do Mar/UFPR
Laboratório de Micropaleontologia
A
Aluizio
Assunto: Re: [R] Arc/Info Ungenerate - shp - Maps
Hi Aluizio
have you tryed to rename .dat to .e00?
sometimes this can help the arcgis understand the format.
I never done this.
good luck
milton
On Wed, Jun 3, 2009 at 8:17 AM, Rodrigo Aluizio wrote:
Hi list.
I would like to
.
-
MSc. <mailto:r.alui...@gmail.com> Rodrigo Aluizio
Centro de Estudos do Mar/UFPR
Laboratório de Micropaleontologia
Avenida Beira Mar s/n - CEP 83255-000
Pontal do Paraná - PR - Brasil
[[alternative HTML version d
Well, on windows an excellent and free software for R scrip editing is Tinn-R.
Works on XP and Vista and it is constantly reviewed and updated.
Rodrigo.
-Mensagem original-
De: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] Em nome
de Peter Flom
Enviada em: terça-fei
200],Pontos$lat[Spp10ml[,i]>=51 &
Spp10ml[,i]<=200],pch=23,bg='yellow',col='black',cex=1.5)
points(Pontos$long[Spp10ml[,i]>=201 &
Spp10ml[,i]<=400],Pontos$lat[Spp10ml[,i]>=201 &
Spp10ml[,i]<=400],pch=22,bg='orange',col='black',cex=1.
nsagem original-
De: Greg Snow [mailto:greg.s...@imail.org]
Enviada em: sexta-feira, 24 de abril de 2009 18:23
Para: Rodrigo Aluizio; R Help
Assunto: RE: [R] Text Contrast in a Plot
Try this function as one approach:
shadowtext <- function(x, y=NULL, labels, col='white', bg=
5,-4.56,-36.111,-4.565,density=NULL,angle=45,col='white',border='
black',lwd=0.1)
text(c(-36.20,-36.155,-36.111,-36.215),c(-4.573,-4.573,-4.573,-4.562),c('0',
'5','10','km'),col='white',cex=1)
Norte(-36.155,-4.60,0.02)
box(which=
Hi Sebastien, take a look at the par(pty='s') argument. Maybe its can solve
your issue.
Best wishes.
Rodrigo.
-Mensagem original-
De: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] Em
nome de Sebastien Bihorel
Enviada em: segunda-feira, 20 de abril de 2009 12:26
Para:
ilto:palsp...@hortresearch.co.nz]
Enviada em: segunda-feira, 23 de março de 2009 19:55
Para: Rodrigo Aluizio; R Help
Assunto: RE: [R] Plot Means Line with Standard Deviation as "Whiskers"
Tena koe Rodrigo
You could always make up some data and then show us what you have tried to
do. I would guess y
lui...@gmail.com> Rodrigo Aluizio
Centro de Estudos do Mar/UFPR
Laboratório de Micropaleontologia
Avenida Beira Mar s/n - CEP 83255-000
Pontal do Paraná - PR - Brasil
Fone: (0**41) 3455-3620 ramal 217
Fax: (0**41) 3455-3623
[[alternative HTML version d
Luis, take a look at the ?legend function, it can produce almost any kind of
legend you may need!
Hope It helps.
Rodrigo.
-Mensagem original-
De: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] Em
nome de Luis Ridao Cruz
Enviada em: terça-feira, 17 de março de 2009 11:
Take a look on it at the CRAN web site.
Hope it helps.
Best wishes.
__
Rodrigo Aluizio
Em 15/03/2009, às 00:27, "Mathew, Abraham T" escreveu:
I'm an R newbie and had a question about the read.xls function. I've
heard that this is often not a reliable functi
Thanks a lot Yihui, It's perfect, in fact exactly identical the I have as a
png file in my computer.
Thank you so much.
Regards
Rodrigo
-Mensagem original-
De: Yihui Xie [mailto:xieyi...@gmail.com]
Enviada em: quarta-feira, 11 de março de 2009 02:48
Para: Rodrigo Aluizio
Cc: R
='2'],col='yellow',cex=0.5,font=2)
text(Coord$Longitude[Coord$Réplicas=='3'],Coord$Latitude[Coord$Réplicas=='3'
],rownames(Coord)[Coord$Réplicas=='3'],col='blue',cex=0.5,font=2)
points(Coord$Longitude,Coord$Latitude-0.045,pch=20,cex=0.7)
ample')
odbcClose(Dados)
library(xlsReadWrite)
setwd('C:/Documents and Settings/Administrador/Meus
documentos/Artigos/Disserta/')
Abio<-read.xls('AF.xls',sheet=1,rowNames=T)
Bio<-read.xls('AF.xls',sheet=2,rowNames=T)
SAbio<-decostand(Abio,method='stan
riginal-
De: Prof Brian Ripley [mailto:rip...@stats.ox.ac.uk]
Enviada em: domingo, 4 de janeiro de 2009 11:15
Para: Rodrigo Aluizio
Cc: R Help
Assunto: Re: [R] Customized LDA (MASS) object plot
Hmm: what do you want to do that a customized panel function cannot do?
plot.lda does not accept typ
question is: How could I reach these axis coordinates inside my
object (supposing they are there somewhere, once the plot.lda and pairs.lda
are able to plot them)?
Thanks in advance.
___
MSc. <mailto:r.alui...@gmail.com> Rodrigo Aluizio
Centro de Estudos do Ma
Also, you can take a look at ?matplot function, depending on the type of
data and graphic output you want, it may b useful.
Best Wishes
___
MSc. Rodrigo Aluizio
Centro de Estudos do Mar/UFPR
Laboratório de Micropaleontologia
Avenida Beira Mar s/n - CEP 83255-000
rgb(229,229,229,max=255),add=T)
axis(1,xaxp=c(,,))
axis(2,yaxp=c(,,))
points(Pontos$long,Pontos$lat,pch=20)
text(Pontos$long,Pontos$lat,rownames(Pontos),pos=1,cex=0.5)
text(Prof$Long,Prof$Lat,rownames(Prof),col='darkgray',cex=0.8,srt=-24)
Ps.: Sorry for hiding the maps coordin
Claudia, Thomas and William thank you all so much.
It's exactly what we needed.
Works perfectly.
Rodrigo.
-Mensagem original-
De: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] Em nome de Claudia Beleites
Enviada em: quinta-feira, 27 de novembro de 2008 14:01
Para: r-help@r-project.org
Assun
e
He told me that Statistica 7.0 has this Index, but he doesnt own it. So
Is
there a way on R to do this correctly?
Thanks for the attention.
___
MSc. <mailto:[EMAIL PROTECTED]> Rodrigo Aluizio
Centro de Estudos do Mar/UFPR
Laboratório de Micropaleontol
Take a look at the links you will found on this previous post of the list.
http://tolstoy.newcastle.edu.au/R/help/06/01/18259.html
I myself don't know anything about this subject.
Sorry, but you probably will found what you need there.
Best Wishes
Rodrigo.
-Mensagem original-
De: [EMA
Take a look at the las() function o ?par. It permits you to choice
labels directions.
__
Rodrigo Aluizio
Em 21/11/2008, às 17:14, "Richardson, Patrick" <[EMAIL PROTECTED]
g> escreveu:
!#x000a
Is there any way to change the orientation of the labels on
Sorry, the code is incomplete.
You get a better result this way...
postscript('Circle.eps',paper='special',width=4,height=4)
par(mar=c(0,0,0,0))
plot.new()
points(0.5,0.5,pch=21,cex=50,bg='gray')
dev.off()
-----Mensagem original-
De: Rodrigo Aluizio [mailto
A fast and simple way to do that would be something like this (the example
is a gray circle)
postscript('Circle.eps')
par(mar=c(0,0,0,0))
plot.new()
points(0.5,0.5,pch=21,cex=50,bg='gray')
dev.off()
You just have to chance the symbol (pch) and color (bg), the output will be
an EPS file, open it o
m original-
De: David Winsemius [mailto:[EMAIL PROTECTED]
Enviada em: sábado, 15 de novembro de 2008 15:57
Para: Rodrigo Aluizio
Cc: R Help
Assunto: Re: [R] PostScript File Dimensions
You may be new to R but that does not mean you should stay "new" to
the help pages. It looks as t
with other
plots Ive done.
Thank you for the attention.
___________
MSc. <mailto:[EMAIL PROTECTED]> Rodrigo Aluizio
Centro de Estudos do Mar/UFPR
Laboratório de Micropaleontologia
Avenida Beira Mar s/n - CEP 83255-000
Pontal do Paraná - PR - BRASIL
[[alternative HTML version
Philipp you actually solved my problem when you mentioned matplot function.
I hadn't read anything about it before, when you mentioned I searched and
I'm getting what I need, now it's just a question of minor adjustments.
Below is a simple solution that I will now improve.
matplot(Perfil$3.W,Perf
AIL PROTECTED]> Rodrigo Aluizio
Centro de Estudos do Mar/UFPR
Laboratório de Micropaleontologia
Avenida Beira Mar s/n - CEP 83255-000
Pontal do Paraná - PR - BRASIL
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
ch + Line, random = ~ 1 |
Block/Strata)
T.lmemod <- lme(response ~ (Season + Beach + Line)^2, random = ~ 1 |
Block/Strata)
anova(T.lmemod)
?anova.lme
## Use package multcomp for doing post-hoc analysis.
Regards, Mark.
Rodrigo Aluizio wrote:
>
> Hi list, I was searching how to properly write
, and be able to run post hoc
tests.
Thanks for your attention and patience.
___
MSc. <mailto:[EMAIL PROTECTED]> Rodrigo Aluizio
Centro de Estudos do Mar/UFPR
Laboratório de Micropaleontologia
Avenida Beira Mar s/n - CEP 83255-000
Pontal do Paraná - PR -
Sorry for don't make the data available. My bad.
But anyway the Richard's fourth solution works quite well.
Thank you very much!
Rodrigo.
-Mensagem original-
De: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED]
Enviada em: quarta-feira, 5 de novembro de 2008 10:59
Para: Rodrigo
='solid',horiz=T,beside=T,las=1,col=c('lightgray','bl
ack'),main='Main Speceis')
legend('topright',c('Living Fauna','Dead
Fauna'),fill=c('lightgray','black'),bty='n')
Thanks for your attention
g and searching.
Rodrigo.
--
From: "Jim Lemon" <[EMAIL PROTECTED]>
Sent: Saturday, October 25, 2008 9:16 AM
To: "Rodrigo Aluizio" <[EMAIL PROTECTED]>
Cc: "R Help"
Subject: Re: [R] Barplot Labels Problem
Rodrig
,las=1,
legend.text=T,width=5,cex.names=0.8)
Thank you for your attenction.
___
MSc. Rodrigo Aluizio
Centro de Estudos do Mar/UFPR
Laboratório de Micropaleontologia
Avenida Beira Mar s/n - CEP 83255-000
Pontal do Paraná - PR - BRASIL
[[alternative HTML version deleted]]
___
x27;,cex=0.6,lwd=0.5,lty='dotted')
Well it's data manipulation, I won't use that...
So... Any one knows how to move the labels, keeping the vectors position on a
cca-like biplot?
Thanks for your attention!
.5,lty='dotted')
Well it's data manipulation, I won't use that...
So... Any one knows how to move the labels, keeping the vectors position on a
cca-like biplot?
Thanks for your attention!
___
MSc. Rodrigo Aluizio
Centro de Estudos do Mar/UFPR
Labo
2] * -1
plot(vare.cap)
for(i in c(2:8)) vare.cap$CCA[[i]][, 1:2]<-vare.cap$CCA[[i]][, 1:2] * -1
plot(vare.cap)
Cheers, Mark.
Rodrigo Aluizio wrote:
Well this time I have to assume, I'm not understanding what is wrong now.
And I have to say: 'Thank You for your patience', cause
]][,2]<-CAPpotiFTI,"
## Using your example to switch axis 1. Make sure that CAPpotiFT exists
and
is your
## original object. Copy and paste the lines below into your R console
window and run
CAPpotiFTI <- CAPpotiFT
for(i in c(2:8,15)) CAPpotiFTI$CCA[[i]][, 1]<-CAPpotiFTI$CCA[[i]][, 1]
$CCA[[i]][, 2] <- mynew,cca$CCA[[i]][, 2] * -1
plot(mynew.cca)
## This does axis 1
for (i in 2:8) mynew.cca$CCA[[i]][, 1] <- mynew,cca$CCA[[i]][, 1] * -1
Regards, Mark.
Rodrigo Aluizio wrote:
Ok Mark, it worked for the species.
I still get an error if I try with biplot (cn,sp), but it
Just as a curiosity. Is there on R any function able to run a SimPer
(Similarity Percentage Analysis)? It's an analysis of the Primer software. I
guess, but I'm not sure!
I was looking for it on CRAN, but I didn't found anything.
_______
MSc. Rodrigo Al
Temp.data$sites[, 2] <- Temp.data$sites[, 2] * -1
Temp.data$species[, 2] <- Temp.data$species[, 2] * -1
Temp.data$centroid[, 2] <- Temp.data$centroid[, 2] * -1
Do this for all the things you are plotting whose axes you want reversed.
Check this against your original plot to ensure that the mir
object and yax = obj[, 2]
(xax
being obj[, 1]). Obviously, copying isn't necessary.
mynew.cca <- my.cca
mynew.cca$scores$species[, 2] <- my.cca$scores$species[, 2] * -1
mynew.cca$scores$sites[, 2] <- my.cca$scores$sites[, 2] * -1
mynew.cca$scores$centroids[, 2] <- my.cca$scores$centroids
,],pch=20)
points(ScoresCAPFT$sites[FatoresRep$BiofaciesBC=='3',],pch=22,bg='gray')
points(ScoresCAPFT$sites[FatoresRep$BiofaciesBC=='4',],pch=24)
points(ScoresCAPFT$sites[FatoresRep$BiofaciesBC=='5',],pch=24,bg='black')
points(ScoresCAPFT$sites[Fato
,],pch=20)
points(ScoresCAPFT$sites[FatoresRep$BiofaciesBC=='3',],pch=22,bg='gray')
points(ScoresCAPFT$sites[FatoresRep$BiofaciesBC=='4',],pch=24)
points(ScoresCAPFT$sites[FatoresRep$BiofaciesBC=='5',],pch=24,bg='black')
points(ScoresCAPFT$sites[Fato
Nevermind people, I just got it easily
just realized that what I wanted is inside the object.
So...
write.xls(TukeySS$Biofacies,'TukeySS.xls')
Solved the problem.
Thanks anyway.
From: Rodrigo Aluizio
Sent: Thursday, September 18, 2008 6:06 PM
To: R Help
Subject: Saving a TukeyHSD
s not in table format, but there isn't anyway to
create a table with this object?
I wasn't able to find anything about it in the R documentation.
Thanks in advice.
___
MSc. Rodrigo Aluizio
Centro de Estudos do Mar/UFPR
Laboratório de Micropaleontologia
Ave
If the students have the permission, maybe run R as administrator solves the
question too.
Rodrigo.
--
From: "Duncan Murdoch" <[EMAIL PROTECTED]>
Sent: Thursday, September 18, 2008 10:48 AM
To: "Green, Paul" <[EMAIL PROTECTED]>
Cc:
Subject: Re: [
ng about these in order to
use proj4 correctly.
James
On Tue, Sep 16, 2008 at 12:35 PM, Rodrigo Aluizio <[EMAIL PROTECTED]>
wrote:
Hi everybody,
I'm looking for some package or function of R that can convert
Coordinates of Longitude and Latitude in different formats, such like
matics to this one,
but I'll need to convert other formats in the future and it won't be that easy.
So, any suggestion well be appreciated.
Thanks in advice.
___
MSc. Rodrigo Aluizio
Centro de Estudos do Mar/UFPR
Laboratório de Microp
Hi Everybody,
Now I have a doubt that is more statistical than R's technical. Im working
with ecology of recent Foraminifera.
At the lab we used to perform cluster analysis using 1-Pearsons R and Wards
method (we already saw it in bibliography of the area) which renders good
results with our
Hi Everybody,
Now I have a doubt that is more statistical than R's technical. Im working
with ecology of recent Foraminifera.
At the lab we used to perform cluster analysis using 1-Pearsons R and Wards
method (we already saw it in bibliography of the area) which renders good
results with our
through vegan carefully to see the
full
offering. Nevertheless, it's pretty easy to do it from the ground up once
you get the knack.
You would do something like the following. Say you have a factor-column
called Site, the rows of which match the rows of the data frame (=mydf)
you used for the
an always get your scores [via: scores(ord.object): see
?scores] from your saved ordination object and add them as per above.
Regards, Mark.
Rodrigo Aluizio wrote:
Hi, I'm a beginner with R, but I'm getting excellent results with it.
Well I've got an capscale object (vegan package)
",add=T)
dbRDA
#Preparing to generate and save the graphic
SymbolsRep<-read.xls('Fatores.xls',sheet=2)
tiff('dbRDAPontos.tif',width=1250,height=1250,res=150)
plot.cca(dbRDA,type='none',display=c('bp','sites'))
text.cca(dbRDA,dis='cn'
;bray",add=T)
dbRDA
#Preparing to generate and save the graphic
SymbolsRep<-read.xls('Fatores.xls',sheet=2)
tiff('dbRDAPontos.tif',width=1250,height=1250,res=150)
plot.cca(dbRDA,type='none',display=c('bp','sites'))
text.cca(dbRDA,dis=
",add=T)
dbRDA
#Preparing to generate and save the graphic
SymbolsRep<-read.xls('Fatores.xls',sheet=2)
tiff('dbRDAPontos.tif',width=1250,height=1250,res=150)
plot.cca(dbRDA,type='none',display=c('bp','sites'))
text.cca(dbRDA,dis='cn'
85 matches
Mail list logo