r语言初学者,最近在看判别分析,看书上的例题,用的source(“discriminiant.distance.R”),这个题之前讲了一堆这个的程序,是不是把
discriminiant.distance.R的程序写在一个r脚本文件里面,我放在了c盘,然后在r命令框里面写的
source('C:/discriminiant.distance.R'),出现了错误,Error in file(filename, "r", encoding = encoding) : cannot open the connectionIn addition: Warning message:In file(filename, "r", encoding = encoding) : cannot open file 'C:/discriminiant.distance.R': No such file or directory,请问这个是怎么回事,哪块出问题了,以下是粘贴书上的程序放在了一个r文件里面,然后再放在了c盘
discriminiant.distance <- function
(TrnX1, TrnX2, TstX = NULL, var.equal = FALSE){
if (is.null(TstX) == TRUE) TstX <- rbind(TrnX1,TrnX2)
if (is.vector(TstX) == TRUE) TstX <- t(as.matrix(TstX))
else if (is.matrix(TstX) != TRUE)
TstX <- as.matrix(TstX)
if (is.matrix(TrnX1) != TRUE) TrnX1 <- as.matrix(TrnX1)
if (is.matrix(TrnX2) != TRUE) TrnX2 <- as.matrix(TrnX2)
nx <- nrow(TstX)
blong <- matrix(rep(0, nx), nrow=1, byrow=TRUE,
dimnames=list("blong", 1:nx))
mu1 <- colMeans(TrnX1); mu2 <- colMeans(TrnX2)
if (var.equal == TRUE || var.equal == T){
S <- var(rbind(TrnX1,TrnX2))
w <- mahalanobis(TstX, mu2, S)
- mahalanobis(TstX, mu1, S)
}
else{
S1 < -var(TrnX1); S2 <- var(TrnX2)
w <- mahalanobis(TstX, mu2, S2)
- mahalanobis(TstX, mu1, S1)
}
for (i in 1:nx){
if (w > 0)
blong <- 1
else
blong <- 2
}
blong
} 6tOo__O>