 it is very close to zero. To calibrate our independence test, we consider the simulated data from previous section of two copulas, the product copula and the Clayton copula with $\theta=-0.8$. In these cases the best fitting was obtained with a partition of size $K=10$, so we chose the same prior values as for the real data to perform an independence test. Posterior probabilities of $H_0$ are $0.81$ and $0$, respectively, showing that for independent data the posterior probability is a lot larger that 0.5, whereas for clearly dependent data the posterior probability of independence is zero.

Concluding remarks {#sec:concl}
==================

We have proposed a semiparametric Archimedean copula that is flexible enough to capture the behaviour of several families of parametric arquimedean copulas. Our model is capable of modelling positive and negative dependence. The number of parameters in the model to produce a good estimation of the dependence in the data should not be extremely high. For most of the examples considered here 10 parameters are enough.

Defining an appropriate partition to analyse real data sets is not trivial. We suggest to try different values of $\alpha$ in a wide range and compare using a GOF criteria like the LPML we used here.

In the exposition and in examples considered here, we concentrated on bivariate copulas, however extensions to more than two dimensions is also possible, say $C(u_1,\ldots,u_m)=\varphi^{-1}\left(\varphi(u_1)+\cdots+\varphi(u_m)\right)$. Performance of our semiparametric copula in this multivariate setting is worth studying.

Our model is motivated by semiparametric proposals for survival analysis functions and appropriately modified to satisfy the properties of an Archimedean generator. The semiparametric generator presented here turned out to be based on quadratic splines, however, alternative proposals are possible.

Instead of starting with a piecewise function for the derivative of a hazard rate, we could start by defining a piecewise constant function for the hazard rate itself. That is $h(t)=\sum_{k=1}^K\theta_k I(\tau_{k-1}<t\leq\tau_k)$ with $\theta_k>0$, and $\{\tau_k\}$ a partition of the positive real line. In this case the cumulative hazard function becomes $H
it is very close to zero. To calibrate our independence test, we consider the simulated datum from former part of two copulas, the product copula and the Clayton copula with $ \theta=-0.8$. In these event the best fitting was obtained with a division of size $ K=10 $, so we chose the same anterior value as for the real datum to perform an independence test. Posterior probabilities of $ H_0 $ are $ 0.81 $ and $ 0 $, respectively, indicate that for independent data the posterior probability is a lot larger that 0.5, whereas for intelligibly dependent data the posterior probability of independence is zero. 

 conclude remarks { # sec: concl } 
 = = = = = = = = = = = = = = = = = = 

 We have proposed a semiparametric Archimedean copula that is flexible enough to capture the behavior of several families of parametric arquimedean copulas. Our model is capable of model positive and negative dependence. The number of parameters in the exemplar to produce a good estimation of the dependence in the data should not be extremely high. For most of the examples considered here 10 parameters are adequate. 

 specify an appropriate partition to analyze veridical data sets is not trivial. We indicate to try different values of $ \alpha$ in a wide compass and compare using a GOF criteria like the LPML we used here. 

 In the exposition and in examples considered here, we condense on bivariate copulas, however extensions to more than two dimensions is also possible, allege $ C(u_1,\ldots, u_m)=\varphi^{-1}\left(\varphi(u_1)+\cdots+\varphi(u_m)\right)$. operation of our semiparametric copula in this multivariate place setting is worth studying. 

 Our model is motivated by semiparametric proposal for survival analysis functions and appropriately modified to meet the properties of an Archimedean generator. The semiparametric generator present here turned out to be based on quadratic slat, however, alternative marriage proposal are possible. 

 Instead of start with a piecewise function for the derived function of a hazard pace, we could start by defining a piecewise changeless function for the hazard pace itself. That is $ h(t)=\sum_{k=1}^K\theta_k I(\tau_{k-1}<t\leq\tau_k)$ with $ \theta_k>0 $, and $ \{\tau_k\}$ a partition of the positive real line. In this case the cumulative hazard function becomes $ H
 it is very close to zero. Tu calibrate our indepeidence fest, we zonsider the simulated data hrom prevuous section of two cooulas, the product copnla and the Clayvkn copula witg $\theca=-0.8$. In these casex the best fitting was ottxiued with a partition of size $K=10$, so we chose yhf same prior vwluex as ror the real data to perform an insependeice test. Posteroor probabilities of $H_0$ are $0.81$ anf $0$, respectively, shlwing that dor ybdependent dxta the posterior probzbility is a lot larger that 0.5, wfereax for cleatlg fgpendent date the iosterior pronsbilitf of incependence is dero.

Cmncouding remarks {#sec:conrl}
==================

We have proposed a semiparakecric Archimedean copuoa that is xlexkvle enkujh fo caphurx the behavjour of sevwral families of pataivyric arquimesean cjptlas. Our model is capable of modelling kositjve and negative dependwnce. The number of patameters ig the model to produce a good estimation of the dapendxnze nk gye data should not be extremely high. For most os tne examples consldered here 10 paramrtfrx are enough.

Dewining an appropriate partihion to analtse real qata sets is not trivial. We sugtest to try bifderent values of $\appha$ in a wnde ramge amd compare using a GOF eriterja like the LPML we hred here.

In the ebpoxidion and in examples consiqered herx, we eoncentrxted on bidariate cooulas, however extensions to mlra than two dimensions is also possible, sab $C(u_1,\ldots,u_m)=\vatphh^{-1}\lext(\varphi(b_1)+\cdots+\yarphi(u_m)\right)$. Pqrformance of pur semnparamdtric copuma in tiis multivaryate setting hd worth studbing.

Our mjdel is notivatda by semiparamrtric proiofals for survival analnsis whnctions and ap'xokruately modifiec tu swtpsfb the [roperties ox an Arznimeddan gentxatif. Thr semiparametric genaratkr presented here yuvned out jo be basqd on quadratoc splines, however, altecnativx proppsajs are possible.

Instead of starfing with a iiecewise funstiok fot the derivctive of a hazard rate, we could start by defining a piecewuse constant functipk for the havard rwte itselx. That is $h(t)=\sum_{k=1}^K\thera_k I(\tau_{k-1}<t\leq\tau_k)$ with $\theta_k>0$, and $\{\tau_k\}$ a parditioj of the positive real line. In this case the cumulative hazard function bwcomes $H
it is very close to zero. To independence we consider simulated data from the copula and the copula with $\theta=-0.8$. these cases the best fitting was with a partition of size $K=10$, so we chose the same prior values for the real data to perform an independence test. Posterior probabilities of $H_0$ $0.81$ $0$, showing for independent data the posterior probability is a lot larger that 0.5, whereas for clearly dependent the posterior probability of independence is zero. Concluding {#sec:concl} ================== We have a semiparametric Archimedean copula that flexible to capture behaviour several of parametric arquimedean Our model is capable of modelling positive and negative dependence. The number of parameters in the model produce a of the in data not be extremely most of the examples considered here enough. Defining an appropriate partition to analyse real sets is trivial. We suggest to try different of $\alpha$ in a wide range and compare a GOF criteria like the LPML we used here. In the exposition and in examples we concentrated on bivariate however extensions to than dimensions also say $C(u_1,\ldots,u_m)=\varphi^{-1}\left(\varphi(u_1)+\cdots+\varphi(u_m)\right)$. of our semiparametric copula in this multivariate setting is worth studying. model is motivated by semiparametric proposals for survival analysis functions modified satisfy the properties an Archimedean generator. The generator here turned out to on splines, are Instead starting with a piecewise for the derivative of a rate, we could start function for the hazard rate itself. That is I(\tau_{k-1}<t\leq\tau_k)$ with $\theta_k>0$, and $\{\tau_k\}$ a partition the positive real line. In this case the cumulative hazard function becomes
 it is very close to zero. To caliBrate our inDepenDenCe tEsT, we cOnsiDer the simulateD Data From previous section of tWo copUlAS, the PRoDuct cOpula anD ThE cLayToN cOpuLa WItH $\thetA=-0.8$. In These caSes the best FitTiNg was obtaineD WiTh a partitiOn oF size $K=10$, so we chOse The samE pRioR ValueS as For thE real dATa to peRform an inDePEndencE Test. PosTERiOr prObabilities of $H_0$ are $0.81$ ANd $0$, REspectively, shoWing thAt FOr INDepEndEnt data the PoSteriOR probabILiTY IS a lOT larger that 0.5, whEreas for cleARly DependEnT daTA the poSteriOr PRobAbility of inDepeNdence is zEro.

ConCLuding rEMarks {#seC:concl}
==================

we hAve PropOSeD a SemIpARamETrIc ARChiMedean coPuLa That iS fleXIBLE enoUgh To caPture The behaviour oF seVeraL FamIlies Of parAmetRiC arquImedeaN copuLaS. Our model is capaBle oF modellinG poSiTivE aNd negATive dePenDenCe. The nuMber of pARamEtERS In The model to produce a GoOD EsTimation Of the dEPeNdENce in the DaTa sHoulD NOt be eXtreMElY high. For Most of THe ExAmples cOnSidereD hEre 10 ParAmeteRS are Enough.

defining An appROpriate partitiON to analyse reaL DaTA SeTS is nOt tRivial. We sugGest TO try DiffEReNt vALues oF $\alphA$ iN A wIDe range and compare usInG a GOF cRiterIa like the LPML We used here.

iN THe exposiTion ANd IN examples consiDered Here, we concENtrated oN bivaRiate copUlas, howevER ExtensioNs tO moRe tHan TWO dImensions is alSO PossIbLe, say $C(u_1,\LdoTs,u_m)=\varPhi^{-1}\LefT(\vaRphI(u_1)+\Cdots+\varpHi(u_m)\righT)$. PErFoRmAncE of ouR SemiparaMeTriC cOpuLa in tHIs multIvariAte sEtTiNG is Worth stUDyING.

Our MoDeL is mOtiVaTed by SemiPAraMetric pRoposals fOr sURvivAl AnAlysis fUnctions and apPrOpriately mOdIfiEd to saTISfy the prOperties of an Archimedean GEneratoR. ThE semiParaMetric genEraTor preSenTEd here Turned Out to Be BasED On quaDRAtIc sPlInes, howeveR, ALteRnatiVe PropOsals arE possible.

Instead of STarTing with a piecEwiSe fuNCTiOn fOR tHE deRiVAtiVE Of a hazard rate, we Could start By DEfIning a piecEWisE cOnstant FunctioN for tHE hazard Rate itselF. That is $h(t)=\SuM_{k=1}^K\tHETa_k i(\tau_{k-1}<t\leq\tAu_k)$ with $\tHeta_k>0$, and $\{\tAU_k\}$ a paRTiTion oF thE positIvE reAl linE. In thiS CasE the cUmulatIvE hazarD funcTiOn becomeS $H
 it is very close to zero.  To calibr ate o urind ep ende ncetest, we consi d er t he simulated data from  prev io u s se c ti on of  two co p ul a s , t he  p rod uc t  c opula  an d the C layton cop ula  w ith $\theta= - 0. 8$. In the secases the be stfittin gwas obtai ned  with  a par t itionof size $ K= 1 0$, so we chos e th e sa me prior values a s  f o r the real dat a to p er f or m anind ependencete st. P o sterior pr o b a bil i ties of $H_0$  are $0.81$ and  $0$,re spe c tively , sho wi n g t hat for ind epen dent data  the p o sterior probabi lity i s a  lo t la r ge rtha t0 .5, wh ere a s f or clear ly  d epend entd a t a  the  po ster ior p robability of  in depe n den ce is  zero .

C on cludi ng rem arks{# sec:concl}
==== ==== ========= =
We  ha ve  prop o sed asem ipa rametri c Archi m ede an c o pu la that is flexibl ee n ou gh to ca pturet he  b e haviourof  se vera l famil ieso fparametr ic arq u im ed ean cop ul as. Ou rmod elis ca p able  of mo dellingposit i ve and negativ e  dependence.T he n um b er o f p arameters i n th e  mod el t o  p rod u ce agoodes t im a tion of the depende nc e in t he da ta should not  be extrem e l y  high. F or m o st of the example s con sidered he r e 10 par amete rs are e nough.

D e f ining an  ap pro pri ate p ar tition to ana l y se r ea l dataset s is no t t riv ial . W esuggest t o try di ff er en tval ues o f  $\alpha $inawid e ran g e andcompa re u si ng a G OF crit e ri a like  t he  LPM L w eusedhere . 

I n the e xposition  an d  inex am ples co nsidered here ,we concent ra ted  on bi v a riate co pulas, however extensio n s to mo rethantwodimension s i s also  po s sible,  say $ C(u_1 ,\ ldo t s ,u_m) = \ va rph i^ {-1}\left( \ v arp hi(u_ 1) +\cd ots+\va rphi(u_m)\right)$. Per formance of o ursemi p a ra met r ic cop ul a  in t his multivariat e settingis wo rth studyi n g.
O ur mode l is mo tivat e d by se miparamet ric propo sa ls f o r  su rvival ana lysis fu nctions a n d app r op riate lymodifi ed  to  sati sfy th e  pr opert ies of  a n Arch imede an  generat or. The semiparametricgenera tor p res ented her e t u rne d out tobe b ased on qu adr ati c spl ine s , how ever ,  a lte r nativ e pr o posals ar e  p oss i b le .

Insteado f sta rting  wi t h a pi ecew ise function fort he derivativeof a h aza rdr ate,  w e could startbyde f i ning a p ie cewise cons tant fun ct i on fo r thehazard  rate i t s el f . That  is$h( t)=\sum_{ k=1 }^ K \theta_ kI( \ tau_{k -1}< t\ leq\ta u_k)$w ith$ \ theta_k>0$, and$\{\t a u _k\}$ a p artit io n of th e  pos itive real  line. In t his ca se t he cu mulativ ehazard  fu nc tion becom e s $H
 it_is very_close to zero. To_calibrate our_independence_test, we_consider_the simulated data_from previous section_of two copulas, the_product copula and_the_Clayton copula with $\theta=-0.8$. In these cases the best fitting was obtained with a_partition_of size_$K=10$,_so_we chose the same prior_values as for the real_data to_perform an independence test. Posterior probabilities of $H_0$_are_$0.81$ and $0$,_respectively, showing that for independent data the posterior probability_is a lot larger that 0.5,_whereas for clearly_dependent_data_the posterior probability of_independence is zero.

Concluding remarks {#sec:concl}
==================

We have_proposed a semiparametric Archimedean copula that_is flexible enough to capture the behaviour_of several families of parametric arquimedean_copulas. Our model is capable_of modelling_positive and negative dependence. The_number of parameters_in the_model to produce_a good estimation of the dependence_in the data_should not be extremely high. For_most_of the examples_considered_here_10 parameters_are enough.

Defining an_appropriate_partition to_analyse_real data sets is not trivial._We_suggest to try different values of $\alpha$_in a wide range_and_compare using a GOF_criteria like the LPML we_used here.

In the exposition and in_examples considered_here, we_concentrated on bivariate copulas, however extensions to more than two dimensions_is also possible, say $C(u_1,\ldots,u_m)=\varphi^{-1}\left(\varphi(u_1)+\cdots+\varphi(u_m)\right)$. Performance_of our semiparametric copula_in this_multivariate_setting is worth_studying.

Our_model is_motivated by semiparametric proposals for survival analysis_functions and_appropriately modified to satisfy the properties_of an Archimedean generator._The_semiparametric generator presented here turned out_to be based on quadratic splines,_however, alternative proposals are possible.

Instead_of_starting_with a piecewise function for_the derivative of a hazard rate,_we could start_by defining a piecewise constant function for_the_hazard rate itself. That is $h(t)=\sum_{k=1}^K\theta_k_I(\tau_{k-1}<t\leq\tau_k)$_with $\theta_k>0$, and $\{\tau_k\}$ a partition_of_the_positive real line. In this_case the cumulative hazard function becomes_$H