2.6 Organizing groups

Build lists of named group factors for easiness of reference.

First, create a list aiming to define workshops productions, so no shipwrecks:

factor_list <-
  list(
    Site = factor(cleanAmphorae$Site_Name[!isShipwreck]),
    FabricGroup = factor(cleanAmphorae$FabricGroup[!isShipwreck]),
    ChemGroup = factor(cleanAmphorae$ChemReferenceGroup[!isShipwreck]),
    ProvGroup = factor(ProvenanceGroup[!isShipwreck])
  )

	Site	FabricGroup	ChemGroup	ProvGroup
ACM001	Ca L’Arnau-Can Pau Ferrer	ACM-2	ACM-C	ACM
ACM097	Ca L’Arnau-Can Pau Ferrer	ACM-1	ACM-A	ACM
BIF006	BDN-Pompeu Fabra	BIF-1	BIF-2	BIF
BIF048	BDN-Pompeu Fabra	BIF-1	BIF-3	BIF
CRC002	Cal Ros de les Cabres	CRC-1	CRC-1	CRC
ELV002	El Vilarenc	ELV-1	ELV-1	ELV
ELV051	El Vilarenc	ELV-1	ELV-2	ELV
FEU009	Can Feu	FEU-IND1	FEU-IND1	FEU-IND1
LLA010	Llafranc	LLA-1	LLA-A2	LLA
MOR018	El More	MOR-2	MOR-3	MOR
SAL025	La Salut	SAL-2	SAL-1	SAL
SBL045	Sant Boi Historic Centre	SBL-2	SBL-2	SBL

Then, create a second list, aiming to assign shipwreck observations to workshop productions, so with shipwreck samples but no true outliers:

factor_list_Shipwreck <-
  list(
    Site = factor(cleanAmphorae$Site_Name[!isTrueIND]),
    FabricGroup = factor(cleanAmphorae$FabricGroup[!isTrueIND]),
    ChemGroup = factor(cleanAmphorae$ChemReferenceGroup[!isTrueIND]),
    ProvGroup = factor(ProvenanceGroup[!isTrueIND])
  )