• Nie Znaleziono Wyników

ANALIZA WRAŻLIWOŚCI KRZYWEJ PHILLIPSA Z ENDOGENICZNĄ CZĘSTOTLIWOŚCIĄ WYZNACZANIA CEN NA ROZKŁAD KOSZTÓW ZMIANY MENU

N/A
N/A
Protected

Academic year: 2021

Share "ANALIZA WRAŻLIWOŚCI KRZYWEJ PHILLIPSA Z ENDOGENICZNĄ CZĘSTOTLIWOŚCIĄ WYZNACZANIA CEN NA ROZKŁAD KOSZTÓW ZMIANY MENU"

Copied!
8
0
0

Pełen tekst

(1)

Maciej Malaczewski

Uniwersytet Łódzki

Wydział Ekonomiczno-Socjologiczny Katedra Ekonometrii

mmalaczewski@uni.lodz.pl

ANALIZA WRAŻLIWOŚCI KRZYWEJ PHILLIPSA Z ENDOGENICZNĄ

CZĘSTOTLIWOŚCIĄ WYZNACZANIA CEN NA ROZKŁAD KOSZTÓW ZMIANY MENU

Streszczenie: W nowokeynesistowskich modelach klasy DSGE, które są wartościowym narzędziem analizy wpływu polityki pieniężnej na sytuację gospodarczą, istotny element stanowi krzywa Phillipsa. Jej estymacja, w sytuacji gdy w modelu założona jest endogenicz- na aktualizacja cen, nie należy do najprostszych. Celem niniejszej pracy jest rozważenie wrażliwości krzywej Phillipsa z endogeniczną na zmianę przyjętego rozkładu kosztów zmia- ny ceny. Okazuje się, że przyjęcie w rozważaniach rodziny rozkładów o tangensokształtnych dystrybuantach, z wyjątkiem skrajnych przypadków, zapewnia relatywną niewrażliwość ilości powstałych warstw na drobne zmiany w parametrach rozkładu.

Słowa kluczowe: sztywność cen, endogeniczna aktualizacja cen, koszty zmiany menu, krzywa Phillipsa, DSGE.

Wprowadzenie

Estymacja parametrów krzywej Phillipsa z endogeniczną częstotliwością zmiany cen, wyprowadzonej w pracy Bakhshiego, Khana i Rudolf [2007], nie jest łatwa – w krzywej tej występują zmienne nieobserwowalne empirycznie

Niniejszy artykuł powstał w ramach projektu sfinansowanego ze środków na naukę w latach 2010-2012 jako projekt badawczy pt. „Modelowanie i prognozowanie inflacji w Polsce przy pomocy modeli z endogeniczną częstotliwością zmiany cen” (nr N N111 209439, umowa nr 2094/B/H03/2010/39). Autor chciałby w tym miejscu podziękować pozostałym uczestnikom projektu, dr Pawłowi Baranowskiemu, dr Mariuszowi Górajskiemu i dr Grzegorzowi Szafrań- skiemu, za wspólną pracę i liczne uwagi merytoryczne.

(2)

oraz nieskończone ilości opóźnień. Aby podjąć decyzję o tym, jaki ostatecznie rząd opóźnień przyjąć, konieczne jest dokonanie wstępnych wyliczeń. By tego dokonać, niezbędne jest przyjęcie pewnej postaci rozkładu kosztów zmiany ceny.

Celem niniejszej pracy jest zbadanie wrażliwości wyników generowanych przez krzywą Phillipsa z endogeniczną częstotliwością wyznaczania cen na zmiany rozkładu kosztów zmiany menu. Eksperyment polega na obserwacji zmian zachodzących w strukturze równowagowego rozkładu firm w poszcze- gólnych warstwach dla mechanizmu cenotwórczego spowodowanych zadaniem odmiennego rozkładu kosztów menu. Rozkład ten będzie pochodzić z pewnej ogólnej rodziny, zaproponowanej w pracy Dotseya, Kinga i Wolmana [1999]

i stosowanej następnie w kolejnych badaniach innych autorów [np. Bakhshi, Khan, Rudolf, 2007; Baranowski i in. 2014].

1. Krzywa Phillipsa z endogeniczną częstotliwością wyznaczania cen

Najpopularniejsze podejście do opisu mechanizmu cenotwórczego, tzw. sche- mat Calvo [1983], jest oparte na losowych sztywnościach nominalnych – każde przedsiębiorstwo ma jednakowe i stałe w czasie prawdopodobieństwo aktualizacji ceny. Nowa cena jest wyznaczana poprzez maksymalizację strumienia łącznego zdyskontowanego zysku danej firmy. Mikropodstawy tego schematu nie są dla ekonomisty satysfakcjonujące. Kolejne prace doprowadziły do konstrukcji mo- delu opartego na występowaniu kosztów zmiany ceny. Model ten, wprowadzony w pracy Dotseya, Kinga i Wolmana [1999], doczekał się wyprowadzenia krzywej Phillipsa, co umożliwia dołączenie endogenicznego wyznaczania częstotliwości zmiany cen do nowokeynesistowskiego modelu oraz badania empiryczne.

Zakładamy występowanie nieskończonej ilości jednakowych przedsię- biorstw, pomiędzy którymi toczy się konkurencja monopolistyczna. Firmy te w każdym okresie podejmują (lub nie) decyzję o aktualizacji ceny. Każde przed- siębiorstwo staje w obliczu pewnych kosztów zmiany ceny, których rozkład w całej populacji przedsiębiorstw jest dany. Rozkład ten jest stały w czasie i ograniczony przez stałą B, która oznacza maksymalne możliwe do wystąpienia koszty zmiany ceny. Firmy porównują korzyści płynące ze zmiany ceny z wy- stępującymi dla nich kosztami zmiany menu. Jeżeli korzyści przewyższają kosz- ty, firma decyduje się zaktualizować cenę, jeżeli nie − pozostawia cenę aktualną.

Wszystkie przedsiębiorstwa są zatem podzielone na grupy, nazywanymi dalej

„warstwami”, zawierające te firmy, które zmieniały ceny ostatnio w tym samym

(3)

momencie. Przedsiębiorstwa te są identyczne, dlatego ustalana przez nich w danym okresie cena optymalna jest jednakowa. Występowanie stale pewnego, niezero- wego poziomu inflacji oraz ograniczoność kosztów zmiany ceny implikuje wy- stępowanie skończonej ilość warstw J.

Powyżej opisana idea pozwala dokonać wyprowadzenia krzywej Phillipsa, określanej w literaturze skrótem SDPC (state-dependent Phillips curve). Krzywa ta jest następującej postaci:

=

+ +

+ , , +

+

+

(1)

gdzie oznacza wielkość inflacji w chwili t, stanowi procentowe odchyle- nie realnego kosztu krańcowego od poziomu równowagi w chwili t, , jest wartością bezwzględną odchyleń udziału przedsiębiorstw znajdujących się w war- stwie + 1 w momencie + od ich udziału w stanie równowagi, a Ω jest średnią ważoną tych odchyleń, zaś stanowi tutaj egzogeniczny szok.

, , , oraz są parametrami, które powstają poprzez serię przekształceń1 z podstawowych parametrów modelu Dotseya, Kinga i Wolmana [1999].

W estymacji parametrów krzywej (1) istotne jest określenie liczby wyprze- dzeń i opóźnień zmiennych. Ilość wyprzedzeń jest równa liczbie występujących warstw, co zależy od podstawowych parametrów gospodarki oraz rozkładu kosztów zmiany ceny. O ile jednak kalibracja parametrów wydaje się być wyko- nywalna2, o tyle kształt rozkładu kosztów zmiany ceny do tej pory nie stanowił przedmiotu badań znanych w literaturze. Ilość opóźnień zmiennych po prawej stronie w równaniu (1) wymaga wyznaczenia wartości kolejnych parametrów wpływu na poziom bieżącej inflacji, po czym dokonania wyboru odpowiedniego rzędu opóźnień. Tu także jednak niezbędna jest znajomość nie tylko podstawo- wych parametrów badanej gospodarki, ale także rozkładu kosztów zmiany ceny.

Zmiany jego kształtu owocują wyznaczeniem odmiennego równowagowego rozkładu firm po poszczególnych warstwach, co prowadzi do zmiany wszystkich w zasadzie wartości parametrów. Wydaje się zatem istotne określenie stopnia wpływu, jaki może mieć niedokładne oszacowanie rozkładu kosztów zmiany ceny na równowagowy rozkład firm po poszczególnych warstwach.

1 Szczegółowe wzory wyrażające związki wymienionych parametrów z „głębokimi” parametrami modelu Dotseya, Kinga i Wolmana [1999] są dość złożone, dla przejrzystości wywodu pomi- jamy je w niniejszej pracy. Zainteresowany czytelnik może znaleźć je w pracy [Bakhshi, Khan, Rudolf, 2007, s. 2326-2328].

2 Do parametrów tych należą stopa dyskonta konsumpcji, równowagowy poziom inflacji, udział sił pracy w łącznym produkcie oraz marża przedsiębiorstw.

(4)

2. Dystrybuanta rozkładu kosztów zmiany ceny

Za [Dotsey, King, Wolman, 1999], będziemy rozpatrywać rodzinę dystry- buant tangensokształtnych, o następującej ogólnej postaci3:

( )= + tan( ⋅ − ) (2)

Poprzez ograniczenia nakładane na dystrybuanty rozkładu prawdopodo- bieństwa4 powyższą rodzinę funkcji możemy ograniczyć do pewnej dwuparame- trycznej rodziny. W skrajnych przypadkach wartości parametrów, rozkład kosz- tów zmiany ceny przyjmuje znane, specyficzne postaci, np. dla → +∞, dystrybuanta określa rozkład jednostajny na przedziale [0;B], dla np. = 0 i → 0 – rozkład jednopunktowy, dla ∈ (0; 1) i → 0 − rozkład dwu- punktowy (z zerowym oraz maksymalnym kosztem jako jedynymi możliwymi realizacjami wartości zmiennej losowej). Pomiędzy tymi skrajnościami kryje się nieskończenie wiele możliwości, z których każda może być wskazana jedno- znacznie poprzez podanie wartości dwóch parametrów. Parametr odpowiada głównie za położenie dystrybuanty względem osi rzędnych, a parametr za stopień jej wygięcia.

3. Wyniki

Dla , zmieniających się w przedziale od 10 do 2 dokonano wylicze- nia ilości warstw, na jakie dzielą się przedsiębiorstwa. Wartości parametrów gospodarki ustalono następująco: udział pracy w produkcie na poziomie α = 0,5, maksymalną wartość w rozkładzie kosztów menu – na poziomie 3,75% wyna- grodzeń (tj. 0,75% jednostkowych kosztów pracy), 25% narzut marży brutto, a także przyjęto 4% inflacji (w ujęciu rocznym) w stanie równowagi. Dla każdej pary dobranych parametrów wyznaczono równowagowy rozkład przedsię- biorstw po poszczególnych warstwach. Łącznie dokonano wyliczeń dla 298456

3 Postać ta jest zaproponowana w pracy [Dotsey, King, Wolman, 1999, s. 670].

4 To znaczy chcemy, by funkcja F przyjmowała odpowiednie wartości na krańcach przedziału określoności (F(0) = 0, F(B) = 1) oraz by była niemalejąca (co w przypadku rozpatrywanej klasy funkcji oznacza nieujemność pierwszej pochodnej). Powyższe ograniczenia oznaczają, że z ogólnej postaci dystrybuanty dwa parametry, np. i , możemy wyrazić jako funkcje pozostałych:

= arctan 1 − − arctan −

oraz = − arctan − . Dotsey, King i Wolman [1999] używają następujących wartości parametrów: c1 = 0,1964, c2 = 0,0625, c3 = 367,44, c4 = 1,2626.

(5)

różnych par parametrów i . Kalkulacji dokonano za pomocą algorytmu zaproponowanego przez Bakshiego, Khana i Rudolf [2007], udostępnionego przez B. Rudolf i przetłumaczonego przez M. Górajskiego i G Szafrańskiego na język programowania programu Matlab.

Tab. 1 prezentuje ilość warstw występujących w równowagowym rozkła- dzie przedsiębiorstw. Dla wartości parametru przekraczającego 0,1, ilość warstw jest dość stabilna, przy przyjętej kalibracji ilość ta wynosi 8 lub 9. Ozna- cza to, że dla przeciętnych wartości parametrów i model Dotseya, Kinga i Wolmana [1999] wykazuje się względną odpornością i niewrażliwością na dobór dystrybuanty. Przy przyjmowaniu zatem subiektywnej dystrybuanty rozkładu kosztów zmiany ceny wybór jednej z rodziny tangensokształtnych o wskazanych parametrach nie spowoduje dużej zmiany w ilości przyspieszeń w krzywej SDPC.

Oznacza to także, że dobór odpowiedniej ilości przyspieszeń po prawej stronie jest stosunkowo niewrażliwy na wybór rozkładu kosztów zmiany ceny.

Rys. 1-4 zawierają wykresy równowagowych rozkładów przedsiębiorstw po poszczególnych warstwach dla różnych wartości parametrów i . Ze względu na mnogość wszystkich możliwych kombinacji tych parametrów (a zatem na mno- gość różnych dystrybuant rozkładów kosztu menu) pogrupowano je celem zacho- wania przejrzystości i zaprezentowano kilka wybranych5. Ze wszystkich rozkładów przynależnych do danej grupy przedstawiono rozkład medianowy oraz rozkłady stanowiące 5 i 95 percentyl dla każdej warstwy.

Tabela 1. Ilość warstw w zależności od przyjętej parametryzacji dystrybuanty

0,001

−0,1 0,101

−0,2 0,201

−0,3 0,301

−0,4 0,401

−0,5 0,501

−0,6 0,601

−0,7 0,701

−0,8 0,801

−0,9 0,901

−1 1,001

−2 0,001-0,1 3-14 8 8 8 8 8 8 8 8 8 8 0,101-0,2 3-14 8 8 8 8 8 8 8 8 8 8 0,201-0,3 3-15 8 8 8 8 8 8 8 8 8 8 0,301-0,4 7-15 8 8 8 8 8 8 8 8 8 8 0,401-0,5 7-15 8 8 8 8 8 8 8 8 8 8 0,501-0,6 7-15 8 8 8 8 8 8 8 8 8 8 0,601-0,7 7-16 8 8 8 8 8 8 8 8 8 8 0,701-0,8 7-17 8 8 8 8 8 8 8-9 8–9 8-9 8-9 0,801-0,9 8-16 8 8-9 8-9 8-9 8-9 8–9 9 9 9 9 0,901-1 8-17 8-9 8-9 9 9 9 9 9 9 9 9 1,001-2 8-9 9 9 9 9 9 9 9 9 9 9

Źródło: Opracowanie własne na podstawie obliczeń.

5 Szczegółowe wyniki symulacji oraz niezbędne kody programu Matlab są dostępne na życzenie u autora niniejszej pracy. Punktem startowym tego algorytmu jest przyjęcie równomiernego rozkładu firm po warstwach, a następnie iterowanie ruchu przedsiębiorstw pomiędzy warstwa- mi zgodnie z ideą mechanizmu cenotwórczego [Dotsey, King, Wolman, 1999].

(6)

R

Ź

R

Ź

R

Ź Rys.

Źród

Rys.

Źród

Rys.

Źród . 1. M

w dło: O

. 2. M w dło: O

. 3. M w dło: O

Med w st Opra

Med w st Opra

Med w st Opra

dian tanie

cowa

dian tanie

cowa

dian tanie

cowa na or

e rów anie

na or e rów

anie

na or e rów

anie raz wno wła

raz wno wła

raz wno wła

5 i owag sne n

5 i owag sne n

5 i owag sne n

95 gi dl na po

95 gi dl na po

95 gi dl na po

perc la pa odsta

perc la pa odsta

perc la pa odsta

cent aram awie

cent aram awie

cent aram awie

tyl r metró e obli

tyl r metró e obli

tyl r metró e obli

rozk ów iczeń

rozk ów iczeń

rozk ów iczeń

kładu , ń.

kładu , ń.

kładu

ń.

u pr

u pr

u pr

∈ (0, rzed

(0,2

rzed (0,0

rzed ,1; 1

dsięb 2; 0,

dsięb 001

dsięb 1),

biors ,8)

biors

; 0,1

biors

stw

stw 1)

stw (0,0

po

po

po 001;

posz

posz

posz

; 0,1 zcze

zcze

zcze )

ególn

ególn

ególn nych

nych

nych h w

h w

h w warst

warst

warst twac

twac

twac ch

ch

ch

(7)

R

Ź

w b p i ż p t n z r n k r d r t n o s n ż

6

Rys.

Źród

weg bior pun i um żon prze trow nied zmi rozk nie ku t rozk duż rozk tryb nej obra stąd na d żon

6 T ch ko or

. 4. M w dło: O

G go, rstw nkto miar ny d eds wan daw iany kład to j taki kład że w kład buan pro azo d wn

drob ny d

wier hani oszt raz n

Med w st Opra

Gdy sko w je owe rko do ięb ny n wno

y c d ce jest ich dem wah

du nta opo owa nio bne do ro

rdzą izmu tów

niesk dian tanie

cowa

roz onc est z go, owa

roz iors na o –

eny en p t zg

roz m 5 hani kos roz ozyc ane osek e m ozk

ą on u ce zmi koń

na or e rów

anie

zkła cent zbli roz anyc zkła stw

wa nie y ro przy godn zkła

-pe ia i sztó zkła cji D

prz k o w mody

kład

ni bo enotw

iany czon

raz wno wła

ad k trow iżon zkła ch w adu w po

arstw e da ozkł ybi ne z adó erce wz ów

adó Dot zez wzg yfik du d

owie wórc y cen nośc

5 i owag sne n

kos wan ny d ad t war u je o po

wac alej

ładu era z w w, enty

zglę zm ów k

tsey mi ględ kacj dwu

em, czeg ny z ci.

95 gi dl na po

ztów nego

do r ten b rtoś edno

oszc ch z niż u z a po wyw mo ylow

ędn mian kos ya, iary dne je r upun

że m go D zosta

perc la pa odsta

w z o n rów był ści p

osta cze zaw ż 3 zbli ostać woda ożem

wym ną w ny c sztu Kin y po ej o rozk

nkt

mec Dots anie

cent aram awie

zmi na w wno

łby para ajne egól wier ok żon ć p ami my m a wraż

ceny u zm nga ozy dpo kład tow

hani seya e prz

tyl r metró e obli

iany war omie rów ame ego lnyc rają kres nego

odo i Ba zao a ro

żliw y w mian a i W ycyj orno du k wego

izm , Ki zyjęt

rozk ów iczeń

y ce rtoś ern wno etru . W ch ącyc

sy t o d obn aks obse ozkł woś w ta ny Wo jne, ośc kosz o.

cen inga ty ro

kładu

ń.

eny ści

ego omie

u Wid

war ch tem do d ną d

shie erw łade ć ró akim cen olma

, w i m ztów

notw a i W

ozkł u pr

∈ (1;

y jes B, o (ry

erny ro dać rstw firm mu.

dwu do ro ego, wow em

ówn m p ny p

ana wyda mode w z

wórc Wolm

ład rzed

; 2),

st z to ys.

y. D ozkł

(ry wac my, W upu ozk , Kh wać

95 now przy

przy a [1

ają elu zmi

zy C man

dwu dsięb

zbliż rów 2).

Dla ład ys.

ch j , kt prz unkt kład han rel -pe wag ypad

yjm 999

się Do iany

Calv na [1 upun

biors

∈ (0

żon wn

W wy kos 4), est tóre zyp tow du g na i atyw rce gow dku mie p

9] ( ę by otsey y, o

vo je 1999 nkto

stw 0,1; 1

ny d ow

gra ysok sztó , że wó e zm

adk wego

geom Ru wn nty wego

u. D pos (rys yć

ya, o ile

est s 9] w wy,

po 1)

do r ago anic kich ów e ró ówc mie ku p

o (r met udol ie d ylow o ro Dla

stać s. 1 o w

Kin e ro

szcz w pr sko

posz

rozk owy

cy, h wa

zm ówn czas enia przy rys tryc lf [2 duż wym

ozk por ć zb

), w wiel

nga ozkł

zegó rzyp once

zcze

kła y ro

dla arto mian now s ba ały

yję . 3) czn 200 e ró m. O kład

rów bliżo wah

le m a i W ład

ólnym adku entro

ególn

du ozk a ro ości ny c wag ard

cen cia ) ró ego 07]6

óżn Ozn du f wnan

oną hani mni Wo ten

m pr u, g owan

nych

jed kład zkł i pa ceny gow dziej nę r jak ówn o. S

6. W nice nacz firm

nia ą do

ia r iejs olma n nie

rzyp gdy j ny n

h w

dnop d pr

adu aram y je wy r

j sk rela ko r now Spos W pr po za t m na , gd o or rozk sze.

ana e je

padk jako naok

warst

pun rzed u je metr est z

roz kon atyw

roz wag strz rzyp omię to d a do dy rygi kład Pł a [19

est z

kiem o ro koło

twac

nkto dsię dno ru

zbli zkła ncen wni zkła gow zeże

pad ędz doś obó dys inal dów łyni

999 zbli

m me zkła o zer

ch

o- ę- o-

i- ad

n- ie ad wy e- d- zy

ść ór s-

l- w,

ie 9]

i-

e- ad

ra

(8)

Podsumowanie i wnioski

Badanie pokazuje, że przy przyjętej parametryzacji, przybliżającej parame- tryzację gospodarki Polski, równowagowy rozkład przedsiębiorstw po poszcze- gólnych warstwach wykazuje względną niewrażliwość na zmiany w rozkładzie kosztów zmiany ceny. Jak wiadomo, zgodnie ze stanem wiedzy autora niniejszej pracy, w literaturze nie są znane próby szacowania postaci rozkładu kosztów zmiany ceny. Pomijając przypadki przyjęcia skrajnych postaci tegoż rozkładu, wrażliwość równowagowego rozkładu firm po poszczególnych warstwach jest dość niewielka, co oznacza, że ilość opóźnień i przyspieszeń dla zmiennych znajdujących się po prawej stronie równania (1) jest mało zmienna z powodu doboru konkretnej dystrybuanty. Badana była jednak jedynie jedna rodzina roz- kładów, przyjęta w literaturze, a wybrany skok wartości parametrów może oka- zać się zbyt duży dla potrzeb wyciągania ogólnych wniosków.

Literatura

Bakhshi H., Khan H., Rudolf B. (2007), The Phillips Curve Under State-Dependent Pricing, „Journal of Monetary Economics”, Vol. 54, No. 8.

Baranowski P., Górajski M., Malaczewski M., Szafranski G. (2014), Inflacja w modelu z endogeniczną częstotliwością aktualizacji cen, „Ekonomista”, Vol. 1.

Calvo G. (1983), Staggered Prices in a Utility-maximizing Framework, „Journal of Mo- netary Economics”, Vol. 12, No. 3.

Dotsey M., King R., Wolman A. (1999), State-dependent Pricing and the General Equ- ilibrium Dynamics of Money and Output, „The Quarterly Journal of Economics”, Vol. 114, No. 2.

SENSIBILITY OF STATE-DEPENDENT PHILLIPS CURVE TO MENU COST DISTRIBUTION

Summary: Phillips curve is an important part of New-Keynesian DSGE models, which are prominent tool in analysis of effects of monetary policy. Estimation of state- dependent Phillips curve is complex, therefore in this paper we consider sensibility of it to changes in probability distribution of menu cost. We find that when cumulative distri- bution function is assumed to be tangent-shaped number of vintages is relatively invari- ant to small changes in parameters of the distribution.

Keywords: price rigidity, state-dependent, menu cost, Phillips curve, DSGE.

Cytaty

Powiązane dokumenty

Zwiêkszenie ceny pozwolenia na emisjê CO 2 (scenariusz WYS_CO2) skutkuje wzro- stem kosztów produkcji energii, a w konsekwencji wzrostem cen zarówno na rynku dnia nastêpnego (rys..

Rzucamy

From general consideration of the geometry of scattering problem it can be assumed that at low frequencies (VHF band), when the influence of surface scattering effect becomes low

ANALIZA WRAŻLIWOŚCI KRZYWEJ PHILLIPSA Z ENDOGENICZNĄ CZĘSTOTLIWOŚCIĄ WYZNACZANIA CEN NA ROZKŁAD KOSZTÓW. ZMIANY MENU

Drawing on my experience as the surtitle translator of Noda’s Japanese adaptation “back” into English, I discuss the linguistic and cultural metamorphosis of Noda’s reworking

ny jest czasami w celu zróżnicowania dóbr kompensacyjnych w zależności od pre­ ferencji kraju wprowadzającego handel wiązany. Określa on, jaki procent warto­ ści towaru

Jest ponadto kilka listów Jeana Yanoskiego i wzmianki o nim w listach Micheleta do innych osób.. Yanoski był młodym wówczas historykiem pochodzenia polskiego,

This paper is devoted to the cultural and educational activities that promote the culture of artistic associations for old people in Piotrkow – The Association of Friends of