Hi, The below very strange: # a) aov function av <- aov(Sepal.Length ~ Species, data=iris) # Error in parse(text = x) : # unexpected symbol in "Sepal(Sepal.Length+Species)Length" av <- aov(iris[, 1] ~ iris[, 5]) # summary(av) # Df Sum Sq Mean Sq F value Pr(>F) # iris[, 5] 2 63.2 31.6 119 <2e-16 *** # Residuals 147 39.0 0.3 # --- # Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 # b) lm function lm1 <- lm(Sepal.Length ~ Sepal.Width, data=iris) # Error in parse(text = x) : # unexpected symbol in "Sepal(Sepal.Length+Sepal.Width)Length" lm1 <- lm(iris[, 1] ~ iris[, 2]) summary(lm1) # Call: # lm(formula = iris[, 1] ~ iris[, 2]) # # Residuals: # Min 1Q Median 3Q Max # -1.556 -0.633 -0.112 0.558 2.223 # # Coefficients: # Estimate Std. Error t value Pr(>|t|) # (Intercept) 6.526 0.479 13.63 <2e-16 *** # iris[, 2] -0.223 0.155 -1.44 0.15 # --- # Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 # # Residual standard error: 0.82 on 148 degrees of freedom # Multiple R-squared: 0.0138, Adjusted R-squared: 0.00716 # F-statistic: 2.07 on 1 and 148 DF, p-value: 0.152 I'm wrong or is it a bug? platform i386-pc-mingw32 arch i386 os mingw32 system i386, mingw32 status Patched major 2 minor 8.1 year 2009 month 01 day 22 svn rev 47680 language R version.string R version 2.8.1 Patched (2009-01-22 r47680) Thanks, -- ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ Jose Claudio Faria Estatistica - prof. Titular UESC/DCET/Brasil joseclaudio.faria at gmail.com joseclaudio.faria at terra.com.br ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\
There is no bug.
Both of your examples work fine for me. You may have a corrupt
version of ``iris'' somewhere in your search path before
``datasets''.
What does find("iris") tell you? What does names("iris")
tell you?
cheers,
Rolf Turner
On 9/04/2009, at 8:16 AM, Jose Claudio Faria wrote:
> Hi,
>
> The below very strange:
>
> # a) aov function
> av <- aov(Sepal.Length ~ Species, data=iris)
>
> # Error in parse(text = x) :
> # unexpected symbol in "Sepal(Sepal.Length+Species)Length"
>
> av <- aov(iris[, 1] ~ iris[, 5])
> # summary(av)
> # Df Sum Sq Mean Sq F value Pr(>F)
> # iris[, 5] 2 63.2 31.6 119 <2e-16 ***
> # Residuals 147 39.0 0.3
> # ---
> # Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05
'.' 0.1 ' ' 1
>
>
>
> # b) lm function
> lm1 <- lm(Sepal.Length ~ Sepal.Width, data=iris)
>
> # Error in parse(text = x) :
> # unexpected symbol in "Sepal(Sepal.Length+Sepal.Width)Length"
>
> lm1 <- lm(iris[, 1] ~ iris[, 2])
> summary(lm1)
>
> # Call:
> # lm(formula = iris[, 1] ~ iris[, 2])
> #
> # Residuals:
> # Min 1Q Median 3Q Max
> # -1.556 -0.633 -0.112 0.558 2.223
> #
> # Coefficients:
> # Estimate Std. Error t value Pr(>|t|)
> # (Intercept) 6.526 0.479 13.63 <2e-16 ***
> # iris[, 2] -0.223 0.155 -1.44 0.15
> # ---
> # Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05
'.' 0.1 ' ' 1
> #
> # Residual standard error: 0.82 on 148 degrees of freedom
> # Multiple R-squared: 0.0138, Adjusted R-squared: 0.00716
> # F-statistic: 2.07 on 1 and 148 DF, p-value: 0.152
>
> I'm wrong or is it a bug?
>
> platform i386-pc-mingw32
> arch i386
> os mingw32
> system i386, mingw32
> status Patched
> major 2
> minor 8.1
> year 2009
> month 01
> day 22
> svn rev 47680
> language R
> version.string R version 2.8.1 Patched (2009-01-22 r47680)
>
> Thanks,
> --
> ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\
> Jose Claudio Faria
> Estatistica - prof. Titular
> UESC/DCET/Brasil
> joseclaudio.faria at gmail.com
> joseclaudio.faria at terra.com.br
> ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\
>
> ______________________________________________
> R-help at 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.
######################################################################
Attention:\ This e-mail message is privileged and confid...{{dropped:9}}
I'm also suspicious of code using parse to construct expressions. However this problem arose because the terms class is a subclass of the formula class and AlgDesign defines a model.matrix method for the formula class but not for the terms class. Hence model.matrix(terms(...)) invokes model.matrix.formula, which causes the problem because AlgDesign::model.matrix.formula isn't general enough to handle a terms object. If AlgDesign had a model.matrix.terms function that called model.matrix.default or was a copy of model.matrix.default then the problem would go away. Bill Dunlap TIBCO Software Inc - Spotfire Division wdunlap tibco.com>> library(AlgDesign) > >> aov(Sepal.Length ~ Species, data=iris) > > Error in parse(text = x) : > unexpected symbol in "Sepal(Sepal.Length+Species)Length"