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 spam at websense.com