Dear R-helpers,
I am using a vrtest on time series data. My commands are as follows;
read.table("B.txt",sep="\t",fill=TRUE, na.strings = "NA")
require(vrtest)
rm(list=ls(all=TRUE))
datamat <- read.table("B.txt",sep="\t",fill=TRUE, na.strings = "NA")
column <- 1
nob <- nrow(datamat)
y <- log(datamat[2:nob,column])-log(datamat[1:(nob-1),column])
After, the use of last command, I get the following message;
in Math.factor(c(37L, 36L, 42L, 41L, 44L, 38L, 31L, 61L, 66L, 91L, :
log not meaningful for factors
My data is composed of one column.
I would appreciate, if any one could provide me with some hints to get
around the problem.
Best regards,
Ben
BNS CER
"11,14"
"11,11"
"11,25"
"11,22"
"11,31"
"11,15"
"11,00"
"11,88"
"11,94"
"12,44"
"12,90"
"12,94"
"12,97"
"12,80"
"12,78"
"13,14"
"13,72"
"13,12"
"12,93"
"12,62"
"12,62"
"12,22"
"12,30"
"12,12"
"12,16"
"11,89"
"12,30"
"12,37"
"12,50"
"12,50"
"12,41"
"12,72"
"12,54"
"12,44"
"12,73"
"12,94"
"13,04"
"13,16"
"13,60"
"13,58"
"13,65"
"13,67"
"13,80"
"13,87"
"13,33"
"13,32"
"13,42"
"13,51"
"13,63"
"13,90"
"13,68"
"13,57"
"13,43"
"13,61"
"13,84"
"13,48"
"13,39"
"12,83"
"12,54"
"12,59"
"12,57"
"12,15"
"11,93"
"11,98"
"12,09"
"11,91"
"12,27"
"12,12"
"12,00"
"12,35"
"12,69"
"12,83"
"12,59"
"12,64"
"12,88"
"12,90"
"12,37"
"12,97"
"13,37"
"13,41"
"13,56"
"13,59"
"13,52"
"13,53"
"13,33"
"13,22"
"13,36"
"13,43"
"13,40"
"13,45"
"13,43"
"13,36"
"13,38"
"13,27"
"13,10"
"13,02"
"12,81"
"12,70"
"12,69"
"12,66"
"12,60"
"12,64"
"12,67"
"12,67"
"12,69"
"12,81"
"12,67"
"12,36"
"12,40"
"12,39"
"12,64"
"12,87"
"12,85"
"12,85"
"12,87"
"12,92"
"12,90"
"12,88"
"13,01"
"12,70"
"12,32"
"12,53"
"12,19"
"11,98"
"11,94"
"12,00"
"11,80"
"11,92"
"11,88"
"12,00"
"11,80"
"12,09"
"11,80"
"11,82"
"11,65"
"11,42"
"11,16"
"11,40"
"11,03"
"11,04"
"11,10"
"10,61"
"11,10"
"11,50"
"11,45"
"11,35"
"11,50"
"12,03"
"12,20"
"12,13"
"12,65"
"12,79"
"12,52"
"12,79"
"12,79"
"12,56"
"12,58"
"12,33"
"12,10"
"12,55"
"12,00"
"11,74"
"11,76"
"11,85"
"12,03"
"12,27"
"12,81"
"12,47"
"12,06"
"11,92"
"11,59"
"11,80"
"11,47"
"11,59"
"11,50"
"11,61"
"11,63"
"11,35"
"11,14"
"11,00"
"10,88"
"11,00"
"11,01"
"10,90"
"10,89"
"10,74"
"10,45"
"10,30"
"10,56"
"10,85"
"10,60"
"10,83"
"10,45"
"10,29"
"10,30"
"9,75"
"9,55"
"10,23"
"10,72"
"10,22"
"10,77"
"11,25"
"11,74"
"11,43"
"11,20"
"10,83"
"10,78"
"11,00"
"10,63"
"10,15"
"10,80"
"10,78"
"10,16"
"9,45"
"9,50"
"9,91"
"9,25"
"8,68"
"8,94"
"9,30"
"9,75"
"9,10"
"8,28"
"7,95"
"8,26"
"7,60"
"7,96"
"8,78"
"9,64"
"9,60"
"9,61"
"9,76"
"10,00"
"10,34"
"10,60"
"10,54"
"10,61"
"10,53"
"10,73"
"10,43"
"10,63"
"10,51"
"10,39"
"10,46"
"11,30"
"12,21"
"12,06"
"12,33"
"12,58"
"12,68"
"12,80"
"13,53"
"13,55"
"13,30"
"13,25"
"13,53"
"13,40"
"13,06"
"13,01"
"12,83"
"12,91"
"13,07"
"13,09"
"13,90"
"13,43"
"13,61"
"13,22"
"13,36"
"13,23"
"13,38"
"13,28"
"13,18"
"13,90"
"14,31"
"14,34"
"14,58"
"14,08"
"14,45"
"14,18"
"14,15"
"14,35"
"13,90"
"13,66"
"15,13"
"15,15"
"15,55"
"15,89"
"16,08"
"15,74"
"15,63"
"15,69"
"15,66"
"15,28"
"15,22"
"15,44"
"15,08"
"15,02"
"15,23"
"15,50"
"15,25"
"15,17"
"16,20"
"16,95"
"16,65"
"18,05"
"18,65"
"18,93"
"18,98"
"19,85"
"20,26"
"19,65"
"19,08"
"19,78"
"18,83"
"18,55"
"18,38"
"18,80"
"19,57"
"18,88"
"18,65"
"18,88"
"19,83"
"19,43"
"19,56"
"20,15"
"20,18"
"20,18"
"19,25"
"18,98"
"18,90"
"19,10"
"19,59"
"18,83"
"18,85"
"19,13"
"20,05"
"20,31"
"20,65"
"20,25"
"20,45"
"20,80"
"20,90"
"20,63"
"20,65"
"20,30"
"20,10"
"20,05"
"20,25"
"20,00"
"19,98"
"19,90"
"19,60"
"19,74"
"19,83"
______________________________________________
R-help@r-project.org 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.
