Part 2: Apologies, some copy/paste gremlins took over my previous post. Attempt 2 at posting my data/code
Apologies for taking up your screens.
Best regards,
Chris
#####################################################
#### IMPORT DATA {Apologies, it's very messy...}
df <-
structure(
list(
Abratenu = c(
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
19.5943297127603,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
19.5943297127603
),
Austmode = c(
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
137.160307989322,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0
),
Capitell = c(
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
58.7829891382809,
0,
78.3773188510412,
0,
0,
0,
19.5943297127603,
0,
0,
0,
0,
58.7829891382808
),
Cerstdrm = c(
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
19.5943297127603
),
Ceraedul = c(
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
19.5943297127603,
901.339166786972,
0,
58.7829891382808,
0,
0,
0,
0,
0,
7994.4865228062,
0,
0,
0,
0,
5231.68603330699
),
Copepoda = c(
0,
0,
19.5943297127603,
0,
0,
0,
19.5943297127603,
0,
0,
19.5943297127603,
0,
0,
0,
0,
0,
0,
0,
607.424221095569,
0,
0,
19.5943297127603,
0,
0,
0,
0,
0,
0
),
Cumogood = c(
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
19.5943297127603,
0,
0,
0,
0,
0,
39.1886594255206,
0,
0,
0,
0,
0,
0,
0,
0,
0
),
Eteolong = c(
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
19.5943297127603,
0,
0,
0,
0,
0,
19.5943297127603,
58.7829891382809
),
Hausaren = c(
0,
0,
19.5943297127603,
0,
0,
0,
0,
0,
0,
19.5943297127603,
0,
0,
0,
0,
78.3773188510412,
0,
39.1886594255206,
0,
0,
0,
0,
0,
0,
0,
0,
19.5943297127603,
0
),
Laniconc = c(
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
195.943297127603,
0,
0,
0,
0,
156.754637702082
),
Limebalt = c(
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
97.9716485638014,
0,
0,
0,
0,
293.914945691404
),
Littlitt = c(
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
19.5943297127603,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0
),
Melipalm = c(
0,
0,
19.5943297127603,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0
),
Mcrphthl = c(
0,
0,
0,
0,
0,
0,
0,
0,
0,
39.1886594255206,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0
),
Nematoda = c(
97.9716485638015,
999.310815350775,
19.5943297127603,
39.1886594255206,
0,
78.3773188510411,
1097.28246391458,
97.9716485638014,
39.1886594255206,
313.509275404164,
156.754637702082,
19.5943297127603,
0,
0,
19.5943297127603,
0,
0,
3879.67728312653,
0,
39.1886594255206,
97.9716485638014,
0,
78.3773188510412,
235.131956553124,
372.292264542445,
0,
0
),
Nemertea = c(
0,
0,
0,
0,
39.1886594255206,
0,
0,
0,
0,
0,
0,
19.5943297127603,
0,
0,
0,
39.1886594255206,
0,
0,
0,
0,
0,
0,
0,
0,
0,
19.5943297127603,
0
),
Parafult = c(
0,
19.5943297127603,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0
),
Parafulg = c(
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
19.5943297127603,
0,
19.5943297127603,
0,
0,
0,
0,
0,
0,
0,
0,
0
),
Periulva = c(
0,
0,
19.5943297127603,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
2664.8288409354,
0,
431.075253680726,
0,
0,
391.886594255207,
39.1886594255206,
19.5943297127603,
0,
0,
58.7829891382809,
0,
19.5943297127603,
0
),
Petrphol = c(
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
19.5943297127603,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0
),
Phylmuco = c(
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
39.1886594255206,
0,
0,
0,
0,
0
),
Protchae = c(
313.509275404165,
0,
0,
0,
0,
0,
0,
0,
0,
39.1886594255206,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
39.1886594255206,
0,
0,
0,
0
),
Pygoeleg = c(
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
97.9716485638015,
0,
0,
0,
0,
470.263913106248
),
Scolsqua = c(
0,
0,
0,
39.1886594255206,
0,
0,
0,
0,
0,
39.1886594255206,
0,
0,
0,
0,
117.565978276562,
39.1886594255206,
19.5943297127603,
0,
176.348967414843,
274.320615978644,
0,
0,
19.5943297127603,
0,
78.3773188510412,
764.178858797652,
0
),
Scolarmi = c(
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
313.509275404165,
0,
0,
0,
0,
333.103605116925
),
Tharyx = c(
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
19.5943297127603
),
Tubipseu = c(
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
176.348967414842
),
Turbellr = c(
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
39.1886594255206,
0,
0,
0,
0
),
Urotpose = c(
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
19.5943297127603
),
Nephtys = c(
0,
0,
78.3773188510412,
19.5943297127603,
0,
0,
0,
0,
0,
0,
0,
39.1886594255206,
176.348967414843,
0,
0,
0,
0,
58.7829891382808,
0,
0,
117.565978276562,
137.160307989322,
19.5943297127603,
0,
0,
0,
19.5943297127603
),
Pontcrts = c(
0,
0,
0,
0,
0,
0,
0,
0,
0,
19.5943297127603,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0
),
Bathypor = c(
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
39.1886594255206,
39.1886594255206,
19.5943297127603,
0,
0,
0,
0,
0,
0,
0,
0,
19.5943297127603,
0
),
Corophim = c(
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
39.1886594255206,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
0,
58.7829891382808
),
Eurydice = c(
19.5943297127603,
0,
0,
19.5943297127603,
19.5943297127603,
0,
0,
0,
0,
39.1886594255206,
19.5943297127603,
0,
0,
0,
39.1886594255206,
19.5943297127603,
333.103605116925,
0,
0,
0,
0,
0,
19.5943297127603,
0,
0,
235.131956553123,
0
)
),
.Names = c(
"Abratenu",
"Austmode",
"Capitell",
"Cerstdrm",
"Ceraedul",
"Copepoda",
"Cumogood",
"Eteolong",
"Hausaren",
"Laniconc",
"Limebalt",
"Littlitt",
"Melipalm",
"Mcrphthl",
"Nematoda",
"Nemertea",
"Parafult",
"Parafulg",
"Periulva",
"Petrphol",
"Phylmuco",
"Protchae",
"Pygoeleg",
"Scolsqua",
"Scolarmi",
"Tharyx",
"Tubipseu",
"Turbellr",
"Urotpose",
"Nephtys",
"Pontcrts",
"Bathypor",
"Corophim",
"Eurydice"
),
class = "data.frame",
row.names = c(NA,
27L)
)
## Load package
require(clustsig)
set.seed(4656)
sim <- simprof(df, method.transform = "squareroot",
method.distance = "braycurtis")
### if I use Euclidean, I don't get an error, so presumably it's a bray-curtis
thisg?
sim
#####################################################
#####################################################
Information in this message may be confidential and may be legally privileged.
If you have received this message by mistake, please notify the sender
immediately, delete it and do not copy it to anyone else.
We have checked this email and its attachments for viruses. But you should
still check any attachment before opening it.
We may have to make this message and any reply to it public if asked to under
the Freedom of Information Act, Data Protection Act or for litigation. Email
messages and attachments sent to or from any Environment Agency address may
also be accessed by someone other than the sender or recipient, for business
purposes.
To report this email as SPAM,
please forward it to [email protected]
_______________________________________________
R-sig-ecology mailing list
[email protected]
https://stat.ethz.ch/mailman/listinfo/r-sig-ecology
