# There is a preloaded quarterly dataset called 'dailytrades' with 60
# observations. Each observation corresponds to a day and contains the
# total number of buyer-initiated trades ('B') and seller-initiated
# trades ('S') on that day. To know more, type ?dailytrades
xdata <- dailytrades
#--------------------------------------------------------------
# Using generic function pin()
#--------------------------------------------------------------
# Define initial parameters:
# initialset = (alpha, delta, mu, eps.b, eps.s)
initialset <- c(0.3, 0.1, 800, 300, 200)
# Estimate the PIN model using the factorization of the PIN likelihood
# function by Ersan (2016)
estimate <- pin(xdata, initialsets = initialset, verbose = FALSE)
# Display the estimated PIN value
show(estimate@pin)
# Display the estimated parameters
show(estimate@parameters)
# Store the initial parameter sets used for MLE in a dataframe variable,
# and display its first five rows
initialsets <- estimate@initialsets
show(head(initialsets, 5))
Run the code above in your browser using DataLab