Nothing
## ------- calculation of Hessian -------------
snqProfitHessian <- function( beta, prices, weights,
scalingFactors = rep( 1, length( weights ) ) ) {
prices <- unlist( prices ) * scalingFactors
normPrice <- sum( t( prices ) %*% weights )
hessian <- beta / normPrice -
beta %*% prices %*% t( weights ) / normPrice^2 -
weights %*% t( prices ) %*% beta / normPrice^2 +
weights %*% t( weights ) *
mean( ( t( prices ) %*% beta %*% prices ) / normPrice^3 )
return( hessian )
}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.