change the last line m.pr<-predict(mydt.nb, mydata[,-1], type="class")
into m.pr<-predict(mydt.nb, mydata[,-1], type="raw") see ?naiveBayes and ?predict.naiveBayes Your "y" is a continous dependent variable. HTH, weiwei On 1/19/07, leah martell <[EMAIL PROTECTED]> wrote: > Hello > > I have a rather simple code and for some reason it produces an error > message. If someone can tell me why and how to fix it, I would be very > greatful. Thank you in advance. > > > ##### create data > set.seed(10) > n <- 200 # number of training points > n.test <- 200 # number of test points > p<-2 # dimension of input space > z <- matrix(rnorm((n+n.test)*p),ncol=p) > x <- matrix(0,nrow=n+n.test,ncol=p) > for (i in 1:p) > x[,i] <- z%*%rnorm(p) > truecoef <- c(1,2) > prob1 <- exp(x%*%truecoef)/(1+exp(x%*%truecoef)) > # prob is the true probability of class 1 > y <- rbinom(n+n.test,1,prob1) > # separate the data into train and test sets > mydata <- data.frame(y=y[1:n],x=x[1:n,]) > mydata.test <- data.frame(y=y[n+(1:n.test)],x=x[n+(1: n.test),]) > ########################## > library(e1071) > mydt.nb<-naiveBayes(y~ ., data=mydata) > m.pr<-predict(mydt.nb, mydata[,-1], type="class") > > > regards, > Leah > > [[alternative HTML version deleted]] > > ______________________________________________ > R-help@stat.math.ethz.ch mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. > -- Weiwei Shi, Ph.D Research Scientist GeneGO, Inc. "Did you always know?" "No, I did not. But I believed..." ---Matrix III ______________________________________________ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.