Interface do R com códigos compilados

Nas instruções a seguir os comandos precedidos de $ devem ser digidados na linha de comando do Linux (ou análogo em outro sistema operacional), e os precedidos por > devem ser digitados no R.

Códigos escritos e compilados em linguagens C, C + + ou Fortran podem ser chamados de dentro do R conforme ilustramos nos exemplos a seguir. Os passos básicos para tal procedimento são:

1. Escrever o código na linguagem desejada (C, C + + ou Fortran)
2. Compilar o código como o comando do linux (no Windows deve-se usar um comando análogo)

   $ R CMD SHLIB

3. Carregar o código compilado no R com o comando

   > dyn.load()

4. Usar a(s) função(ões) do código compilado com os comandos .C(), .Call() ou .Fortran().

Considere escrever uma função em C para calcular valores da função de correlação de Matèrn. Esta função tem como argumento u e parâmetros <latex>$(\phi, \kappa)$</latex> e é dada pela seguinte expressão:
<latex>$\rho(u) = \{2^{\kappa-1} \Gamma(\kappa)\}^{-1} (u/\phi)^\kappa K_\kappa(u/\phi).$</latex> No arquivo clavras01.c encontra-se um código C para calcular valores de uma versão padronizada desta função onde <latex>$\phi=1$</latex>.

• clavras02.c
• clavras.r