0

这可能是 PICNIC，但是在尝试构建函数时，与在控制台中执行操作时相比，我得到了一组相当奇怪的行为。

我正在尝试使用 cut2 生成一组间隔并将它们分配给一个变量以供以后重用。在函数外部使用会产生适当的向量，但在函数内部，赋值的强制方式不同。我尝试将 as.vector 包裹在它周围，但这产生了一个 char[210]。

谁能告诉我我做错了什么？

谢谢，斯蒂芬 PS 全局分配是这样的间隔可以被修改，并在以后重用

设置

library("Hmisc")
library("caret")
# functions in use ----------------------------------------------------------------
# functions for splitting data according to Max Kuhn's preferences in caret vignette, code condensed to save space
splitDataset<-function(dataset=rawdata, nPrimaryKeyCol=1, nOutcomeCol=1) {
  end<-as.numeric(ncol(dataset))
  stopifnot(is.numeric(nPrimaryKeyCol),nPrimaryKeyCol<=end,is.numeric(nOutcomeCol),nOutcomeCol<=end,(nPrimaryKeyCol+nOutcomeCol)<=end)
  predstart<-nPrimaryKeyCol+1
  predend<-ncol(dataset)-nOutcomeCol
  assign(x="keys",value=dataset[, 0:nPrimaryKeyCol],envir = parent.frame())
  assign(x="outcomes",value=dataset[,(predend+1):end],envir = parent.frame())
  assign(x="predictors",value=dataset[,predstart:predend], envir = parent.frame())
}
    partitionDataset<-function(proportion=0.7){
  require("caret")
  assign(x="inTrain",  value=createDataPartition(outcomes,p=proportion, list=FALSE),  envir = parent.frame())
  assign(x="trainKeys", value=keys[inTrain], envir = parent.frame())
  assign(x="trainPredictors",value=predictors[inTrain,],envir = parent.frame())
  assign(x="trainOutcomes", value=outcomes[inTrain],envir = parent.frame())
  assign(x="testKeys",value=keys[-inTrain],envir = parent.frame())
  assign(x="testPredictors",  value=predictors[-inTrain,],  envir = parent.frame())
  assign(x="testOutcomes",value=outcomes[-inTrain], envir = parent.frame())
  assign(x="trainPredictors.Bad", value=subset(trainPredictors,trainOutcomes=="bad"), envir = parent.frame())
}

# this is the problem function
equalFreqBins.derive<-function(characteristic,deriveDataset,g=20){

  stopifnot(is.numeric(deriveDataset[,characteristic]),is.numeric(g))
  dnam<-paste0("interval.",characteristic)
# this is the assignment
  intervals<-cut2(deriveDataset[,characteristic],g,onlycuts=TRUE)
  assign(
    x=dnam,
    value=intervals,
    envir=as.environment(".GlobalEnv")
    )
# return to make it easier to look at variable
  return(str(intervals))
  }

# data loadup----------------------------------------------------------------
data(GermanCredit)
GermanCredit$outcome<-GermanCredit$Class
GermanCredit$Class<-NULL
levels(GermanCredit$outcome)<-c("bad","good")
basedata<-GermanCredit
splitDataset(basedata,nPrimaryKeyCol=0,nOutcomeCol=1)
partitionDataset(proportion=0.7)

问题再现

# this is outputting a factor of 210
equalFreqBins.derive(characteristic="Age",deriveDataset=trainPredictors.Bad,g=20)
# this is outputting a num [1:20]
intervals<-cut2(trainPredictors.Bad[,"Age"],g=20,onlycuts=TRUE)
str(intervals)
4

1 回答 1

1

的概要cut2是

cut2(x, cuts, m, g, levels.mean, digits, minmax=TRUE, oneval=TRUE, onlycuts=FALSE)

所以当你打电话时

intervals <- cut2(deriveDataset[,characteristic], g, onlycuts = TRUE)

您的g变量，作为第二个未命名的参数，被视为cuts参数。解决方案：使用命名参数调用您的函数：

intervals <- cut2(deriveDataset[,characteristic], g = g, onlycuts = TRUE)
于 2013-04-30T11:35:42.600 回答