Question

我正在查看R中的cov的source_code，并遇到一段我不太了解的代码段。

此处的协方差goes的数学定义。

if (method == "pearson") 
    .Call(C_cov, x, y, na.method, method == "kendall")
else if ...

帮助手册说明了.Call功能：

CallExternal {base} R Documentation
Modern Interfaces to C/C++ code
Description
Functions to pass R objects to compiled C/C++ code that has been loaded into R.

我想知道在哪里可以找到如何计算C ++或C等协方差的源代码。

感谢。

Answer 1

.Call用于将变量传递给C例程。 C_cov是一个变量（在我们很快就会看到的stats命名空间中）告诉.Call在哪里找到它应该用来计算协方差的例程。

如果您在命令行键入C_cov，则会获得

Error: object 'C_cov' not found

那是因为它对你隐藏了。你必须做一些侦探工作。

getAnywhere('C_cov')
# 4 differing objects matching ‘C_cov’ were found
# in the following places
#   namespace:stats
# Use [] to view one of them

这告诉我们C_cov名称空间中有一个名为stats的变量（您的输出可能与此略有不同）。让我们试试吧。

stats::C_cov
# Error: 'C_cov' is not an exported object from 'namespace:stats'

显然C_cov不适合公众消费。没关系，我们无论如何都可以得到它：

stats:::C_cov # use three colons to get unexported variables.
# $name
# [1] "cov"
# # blah, blah, blah ...
# $dll
# DLL name: stats
# Filename: C:/Program Files/R/R-3.0.1/library/stats/libs/x64/stats.dll
# Dynamic lookup: FALSE
# # blah, blah, ...

这就是我们想要的信息。它告诉我们它所在的例程和库的名称。现在我们只需要转到C源并跟踪它：.../src/library/stats/src/cov.c

R源代码.Call函数

1 个答案: