מודול 9: מקרה בוחן II — גרזני ברונזה ורשתות אומנות

שבוע 9 — קווי מתאר ותקופת הברונזה

רקע ארכאולוגי

תקופת הברונזה האירופאית והמתכת שלה

תקופת הברונזה באירופה (כ-2300–800 לפנה”ס) הוגדרה על ידי מטלורגיה. ברונזה — סגסוגת נחושת ובדיל — דרשה רשתות חליפין לטווחים ארוכים, מכיוון שמרבצי נחושת ובדיל כמעט ואינם נמצאים יחד. עצם קיומה של התרבות החומרית של תקופת הברונזה מרמז על קשרים.

גרזני כנפיים הם בין הממצאים המתכתיים הנפוצים ביותר מתקופת הברונזה. הם היו כלים פונקציונליים (עיבוד עץ, לחימה) אך גם חפצי יוקרה, אוצרים במטמונים, ואולי שימשו כחפצי ערך דמויי מטבע.

עיצוב ה”כנפיים” — קצוות גבוהים לאורך הצדדים לסיוע בהצמדה לעץ — היה פתרון טכנולוגי שהתפשט ברחבי אירופה בתקופת הברונזה התיכונה. אבל כיצד התפשטות זו קרתה?

קבוצות G3 ו-G4

מאגר הנתונים שלנו מגיע מגרזני כנפיים ממרכז אירופה שסווגו לשתי קבוצות טיפולוגיות (ייעודים לפי כרונו-טיפולוגיה סטנדרטית של תקופת הברונזה):

G3: קבוצה מוקדמת יותר, כ-1600–1400 לפנה”ס. נמצאת בעיקר במרכז אירופה המערבי.
G4: קבוצה מאוחרת יותר, כ-1400–1200 לפנה”ס. פיזור גיאוגרפי רחב יותר.

הטיפולוגיה המסורתית מסווגת אותם לפי שילוב של פרופורציות (גובה הכנפיים, רוחב הלהב יחסית לאורך הכולל) ותכונות עיטוריות. אך הטיפולוגיה היא תמיד פישוט.

השאלה: האם שונות הצורה בגרזנים אלה מכבדת את גבול G3/G4? האם יש שיפוע מורפולוגי רציף? האם יש וריאנטים אזוריים בתוך כל קבוצה?

מדוע זה חשוב היסטורית

שלושה תוצאות אפשריות:

תוצאה א: הפרדה ברורה

G3 ו-G4 מהווים אשכולות מובחנים ולא חופפים במרחב הצורות. זה תומך ברעיון שכל קבוצה מייצגת מסורת ייצור קוהרנטית — אולי סדנה, אזור או תקופה ספציפיים עם תבנית עיצוב משלהם.

המשמעות: מתכנאי תקופת הברונזה שמרו על תקני צורה קפדניים בתוך כל מסורת, אולי דרך רשתות שוליות עם העברה בין-מסורתית מוגבלת.

תוצאה ב: אשכולות חופפים

G3 ו-G4 חופפים במידה ניכרת במרחב הצורות. שונות הצורה רציפה, לא קטגורית.

המשמעות: גבול הטיפוס G3/G4 עשוי להיות מבנה אנליטי מודרני ולא קטגוריה אמית (פנימית) שקהילות תקופת הברונזה הכירו. הצורה השתנתה בהדרגה, לא בקפיצה בדידה.

תוצאה ג: שיפוע גיאוגרפי חוצה קבוצות טיפולוגיות

צירי PC הלוכדים את רוב השונות מתואמים עם מיקום גיאוגרפי, לא עם חברות בקבוצה G3/G4.

המשמעות: שונות הצורה משקפת מסורות ייצור מקומיות שאינן מתאימות לרצף הטיפולוגי — זהות אומנות אזורית הייתה חזקה יותר משינוי לאורך זמן.

אף אחת מהתוצאות אינה “נכונה” או “שגויה.” כל אחת אומרת לנו משהו חשוב על הארגון החברתי והטכנולוגי של תקופת הברונזה.

ניתוח מלא: הדרכה למטלה 4

פתחו את מחברת שבוע 9

→ מחברת Colab — שבוע 9

מחברת זו מרחיבה את EFA משבוע 8 עם השוואת קבוצות, בדיקה סטטיסטית וויזואליזציה גיאוגרפית.

הבדיקה הסטטיסטית

אנחנו משתמשים ב-MANOVA מבוסס הרשאה (אותו הגיון כמו למטבעות, מוחל על ציוני EFA+PCA):

# MANOVA בהרשאה על ציוני PCA של EFA
f_stat, p_val = permutation_manova(scores, labels, n_perm=999, n_pcs=8)

print(f"\nMANOVA בהרשאה (G3 לעומת G4):")
print(f"  F = {f_stat:.3f}")
print(f"  p = {p_val:.3f}")
print(f"  {'מובהק' if p_val < 0.05 else 'לא מובהק'} ב-α = 0.05")

# גודל אפקט: פרופורציית שונות מוסברת
r2 = f_stat / (f_stat + (len(labels) - 2))
print(f"  R² ≈ {r2:.3f} ({r2*100:.1f}% מהשונות מוסבר על ידי הקבוצה)")

פרשנות גודל אפקט

ערך R²	פרשנות
< 0.05	אפקט קטן מאוד — קבוצות כמעט ולא שונות
0.05–0.15	אפקט קטן
0.15–0.35	אפקט בינוני
> 0.35	אפקט גדול — קבוצות שונות מאוד

גם אם הבדיקה מובהקת (p < 0.05), R² קטן אומר שהקבוצות שונות באופן אמין אך לא דרמטי.

ויזואליזציה גיאוגרפית

המטה-דאטה הגיאוגרפית לגרזנים מסופקת כקובץ CSV. כאן אנחנו מקשרים ציוני PC למיקומי מציאה:

!pip install python-bidi -q
import pandas as pd
import matplotlib.pyplot as plt
import matplotlib.cm as cm
import matplotlib.colors as mcolors
from bidi.algorithm import get_display
rtl = get_display  # תיקון כיוון טקסט עברי בגרפים של matplotlib

# טעינת מטה-דאטה גיאוגרפית
geo = pd.read_csv("data/flanged_axes/axes_metadata.csv")
# עמודות: name, group, latitude, longitude, country, site_name

# איחוד עם ציוני PCA
df_geo = pd.merge(df, geo, on='name')

# מיפוי ציון PC1 לצבע (כחול = שלילי, אדום = חיובי)
norm = mcolors.TwoSlopeNorm(vcenter=0,
                              vmin=df_geo['PC1'].min(),
                              vmax=df_geo['PC1'].max())

# מפת פיזור פשוטה (מקורבת)
fig, ax = plt.subplots(figsize=(10, 8))
sc = ax.scatter(df_geo['longitude'], df_geo['latitude'],
                c=df_geo['PC1'], cmap='RdBu_r', s=80,
                edgecolor='black', linewidth=0.5, norm=norm)

# תיוג סמנים עם קבוצה
for _, row in df_geo.iterrows():
    ax.annotate(row['group'], (row['longitude'], row['latitude']),
                fontsize=6, ha='center', va='bottom')

plt.colorbar(sc, label=rtl('ציון PC1 (ציר צורה)'))
ax.set_xlabel(rtl("קו אורך"))
ax.set_ylabel(rtl("קו רוחב"))
ax.set_title(rtl("פיזור גיאוגרפי של ציוני PC1 — גרזני כנפיים"))
plt.tight_layout()
plt.show()

הערה: קובץ axes_metadata.csv יסופק. אם נתונים גיאוגרפיים אינם זמינים לכל הפרטים, נעבוד עם תת-הקבוצה הכוללת קואורדינטות.

שאלות לדיון

האם הצורות המשוחזרות של EFA בקצוות PC1 תואמות למה שהייתם מצפים לו מהסתכלות ויזואלית על גרזני G3 לעומת G4?
אם G3 ו-G4 חופפים במרחב הצורות אך בדיקת ההרשאה מובהקת — כיצד תסבירו את שתי התוצאות?
דמיינו שמצאתם גרזן ברונזה חדש מאתר לא ידוע. כיצד תשתמשו בניתוח זה לסיווגו כ-G3 או G4? עם איזה רמת ביטחון?
אילו נתונים נוספים (ניתוח איזוטופי, הקשר מטמון, תאריכים מוחלטים) היו מחזקים את פרשנות דפוסי הצורה שמצאתם?

מטלה 4 מגיעה השבוע

הערה

הגישו קישור למחברת Colab + PDF של 600 מילה + תמונת מפה.

→ מטלה 4: ניתוח קווי מתאר לגרזני ברונזה

→ מודול 10: מקרה בוחן III — כלי חרס ומסורות קדרות

--- title: "מודול 9: מקרה בוחן II — גרזני ברונזה ורשתות אומנות" subtitle: "שבוע 9 — קווי מתאר ותקופת הברונזה" lang: he --- ## רקע ארכאולוגי ### תקופת הברונזה האירופאית והמתכת שלה תקופת הברונזה באירופה (כ-2300–800 לפנה"ס) הוגדרה על ידי מטלורגיה. ברונזה — סגסוגת נחושת ובדיל — דרשה רשתות חליפין לטווחים ארוכים, מכיוון שמרבצי נחושת ובדיל כמעט ואינם נמצאים יחד. עצם קיומה של התרבות החומרית של תקופת הברונזה מרמז על קשרים. **גרזני כנפיים** הם בין הממצאים המתכתיים הנפוצים ביותר מתקופת הברונזה. הם היו כלים פונקציונליים (עיבוד עץ, לחימה) אך גם חפצי יוקרה, אוצרים במטמונים, ואולי שימשו כחפצי ערך דמויי מטבע. עיצוב ה"כנפיים" — קצוות גבוהים לאורך הצדדים לסיוע בהצמדה לעץ — היה פתרון טכנולוגי שהתפשט ברחבי אירופה בתקופת הברונזה התיכונה. אבל כיצד התפשטות זו קרתה? ### קבוצות G3 ו-G4 מאגר הנתונים שלנו מגיע מגרזני כנפיים ממרכז אירופה שסווגו לשתי קבוצות טיפולוגיות (ייעודים לפי כרונו-טיפולוגיה סטנדרטית של תקופת הברונזה): - **G3**: קבוצה מוקדמת יותר, כ-1600–1400 לפנה"ס. נמצאת בעיקר במרכז אירופה המערבי. - **G4**: קבוצה מאוחרת יותר, כ-1400–1200 לפנה"ס. פיזור גיאוגרפי רחב יותר. הטיפולוגיה המסורתית מסווגת אותם לפי שילוב של פרופורציות (גובה הכנפיים, רוחב הלהב יחסית לאורך הכולל) ותכונות עיטוריות. אך הטיפולוגיה היא תמיד פישוט. השאלה: **האם שונות הצורה בגרזנים אלה מכבדת את גבול G3/G4? האם יש שיפוע מורפולוגי רציף? האם יש וריאנטים אזוריים בתוך כל קבוצה?** --- ## מדוע זה חשוב היסטורית שלושה תוצאות אפשריות: ### תוצאה א: הפרדה ברורה G3 ו-G4 מהווים אשכולות מובחנים ולא חופפים במרחב הצורות. זה תומך ברעיון שכל קבוצה מייצגת מסורת ייצור קוהרנטית — אולי סדנה, אזור או תקופה ספציפיים עם תבנית עיצוב משלהם. **המשמעות:** מתכנאי תקופת הברונזה שמרו על תקני צורה קפדניים בתוך כל מסורת, אולי דרך רשתות שוליות עם העברה בין-מסורתית מוגבלת. ### תוצאה ב: אשכולות חופפים G3 ו-G4 חופפים במידה ניכרת במרחב הצורות. שונות הצורה רציפה, לא קטגורית. **המשמעות:** גבול הטיפוס G3/G4 עשוי להיות מבנה אנליטי מודרני ולא קטגוריה אמית (פנימית) שקהילות תקופת הברונזה הכירו. הצורה השתנתה בהדרגה, לא בקפיצה בדידה. ### תוצאה ג: שיפוע גיאוגרפי חוצה קבוצות טיפולוגיות צירי PC הלוכדים את רוב השונות מתואמים עם מיקום גיאוגרפי, לא עם חברות בקבוצה G3/G4. **המשמעות:** שונות הצורה משקפת מסורות ייצור מקומיות שאינן מתאימות לרצף הטיפולוגי — זהות אומנות אזורית הייתה חזקה יותר משינוי לאורך זמן. **אף אחת מהתוצאות אינה "נכונה" או "שגויה."** כל אחת אומרת לנו משהו חשוב על הארגון החברתי והטכנולוגי של תקופת הברונזה. --- ## ניתוח מלא: הדרכה למטלה 4 ::: {.callout-important} ## פתחו את מחברת שבוע 9 [→ מחברת Colab — שבוע 9](https://colab.research.google.com/github/shaigordin/comparch/blob/2026/morphometrics/notebooks/week-09.ipynb) ::: מחברת זו מרחיבה את EFA משבוע 8 עם השוואת קבוצות, בדיקה סטטיסטית וויזואליזציה גיאוגרפית. ### הבדיקה הסטטיסטית אנחנו משתמשים ב-MANOVA מבוסס הרשאה (אותו הגיון כמו למטבעות, מוחל על ציוני EFA+PCA): ```python # MANOVA בהרשאה על ציוני PCA של EFA f_stat, p_val = permutation_manova(scores, labels, n_perm=999, n_pcs=8) print(f"\nMANOVA בהרשאה (G3 לעומת G4):") print(f" F = {f_stat:.3f}") print(f" p = {p_val:.3f}") print(f" {'מובהק' if p_val < 0.05 else 'לא מובהק'} ב-α = 0.05") # גודל אפקט: פרופורציית שונות מוסברת r2 = f_stat / (f_stat + (len(labels) - 2)) print(f" R² ≈ {r2:.3f} ({r2*100:.1f}% מהשונות מוסבר על ידי הקבוצה)") ``` ### פרשנות גודל אפקט | ערך R² | פרשנות | |--------|---------| | < 0.05 | אפקט קטן מאוד — קבוצות כמעט ולא שונות | | 0.05–0.15 | אפקט קטן | | 0.15–0.35 | אפקט בינוני | | > 0.35 | אפקט גדול — קבוצות שונות מאוד | גם אם הבדיקה מובהקת (p < 0.05), R² קטן אומר שהקבוצות שונות *באופן אמין* אך לא *דרמטי*. --- ## ויזואליזציה גיאוגרפית המטה-דאטה הגיאוגרפית לגרזנים מסופקת כקובץ CSV. כאן אנחנו מקשרים ציוני PC למיקומי מציאה: ```python !pip install python-bidi -q import pandas as pd import matplotlib.pyplot as plt import matplotlib.cm as cm import matplotlib.colors as mcolors from bidi.algorithm import get_display rtl = get_display # תיקון כיוון טקסט עברי בגרפים של matplotlib # טעינת מטה-דאטה גיאוגרפית geo = pd.read_csv("data/flanged_axes/axes_metadata.csv") # עמודות: name, group, latitude, longitude, country, site_name # איחוד עם ציוני PCA df_geo = pd.merge(df, geo, on='name') # מיפוי ציון PC1 לצבע (כחול = שלילי, אדום = חיובי) norm = mcolors.TwoSlopeNorm(vcenter=0, vmin=df_geo['PC1'].min(), vmax=df_geo['PC1'].max()) # מפת פיזור פשוטה (מקורבת) fig, ax = plt.subplots(figsize=(10, 8)) sc = ax.scatter(df_geo['longitude'], df_geo['latitude'], c=df_geo['PC1'], cmap='RdBu_r', s=80, edgecolor='black', linewidth=0.5, norm=norm) # תיוג סמנים עם קבוצה for _, row in df_geo.iterrows(): ax.annotate(row['group'], (row['longitude'], row['latitude']), fontsize=6, ha='center', va='bottom') plt.colorbar(sc, label=rtl('ציון PC1 (ציר צורה)')) ax.set_xlabel(rtl("קו אורך")) ax.set_ylabel(rtl("קו רוחב")) ax.set_title(rtl("פיזור גיאוגרפי של ציוני PC1 — גרזני כנפיים")) plt.tight_layout() plt.show() ``` **הערה:** קובץ axes_metadata.csv יסופק. אם נתונים גיאוגרפיים אינם זמינים לכל הפרטים, נעבוד עם תת-הקבוצה הכוללת קואורדינטות. --- ## שאלות לדיון 1. האם הצורות המשוחזרות של EFA בקצוות PC1 תואמות למה שהייתם מצפים לו מהסתכלות ויזואלית על גרזני G3 לעומת G4? 2. אם G3 ו-G4 חופפים במרחב הצורות אך בדיקת ההרשאה מובהקת — כיצד תסבירו את שתי התוצאות? 3. דמיינו שמצאתם גרזן ברונזה חדש מאתר לא ידוע. כיצד תשתמשו בניתוח זה לסיווגו כ-G3 או G4? עם איזה רמת ביטחון? 4. אילו נתונים נוספים (ניתוח איזוטופי, הקשר מטמון, תאריכים מוחלטים) היו מחזקים את פרשנות דפוסי הצורה שמצאתם? --- ## מטלה 4 מגיעה השבוע ::: {.callout-note} הגישו קישור למחברת Colab + PDF של 600 מילה + תמונת מפה. [→ מטלה 4: ניתוח קווי מתאר לגרזני ברונזה](../assignments/04-axes.qmd) ::: [→ מודול 10: מקרה בוחן III — כלי חרס ומסורות קדרות](10-pottery.qmd)