[R] (no subject)

Peter Ehlers ehlers at ucalgary.ca
Tue May 18 18:24:55 CEST 2010


Arantzazu,

Your problem is that the data were probably imported
from Excel where you had at least one cell containing "#N/A".
You need to replace those cases in your dataframe with NA.
Then you should be able to do as.numeric(as.character(arcilla)).

  -Peter Ehlers

On 2010-05-18 10:07, Arantzazu Blanco Bernardeau wrote:
>
> Hello
> Well, the problem is, that arcilla is the percentage of clay in the soil sample. So, for linear model, I need to work with that number or value. Now, R thinks that arcilla (arcilla means clay in spanish), is a factor, and gives me the value as a factor, so the output of the linear model is
> Call:
> lm(formula = formula, data = caperf)
>
> Residuals:
>         Min         1Q     Median         3Q        Max
> -1.466e+01 -1.376e-15  1.780e-16  2.038e-15  1.279e+01
>
> Coefficients:
>                Estimate Std. Error t value Pr(>|t|)
> (Intercept)    1.68964    6.33889   0.267 0.790221
> arcilla0.9     1.90228    8.90888   0.214 0.831239
> arcilla10      1.26371    7.96734   0.159 0.874212
> arcilla10.3   15.70081    9.05141   1.735 0.085090 .
> arcilla10.4    7.27517    7.72806   0.941 0.348183
> arcilla10.45   7.03879    9.02600   0.780 0.436853
> arcilla10.5    2.41241    8.90827   0.271 0.786954
> arcilla10.65  15.44298    9.03879   1.709 0.089838 .
> arcilla10.7   19.35651    9.04675   2.140 0.034185 *
> arcilla10.9    3.55947    9.18501   0.388 0.698974
>
> [...]
>
> arcilla9.9     6.31949    7.35724   0.859 0.391892

===> arcilla#N/A   24.17959    8.87201   2.725 0.007274 **

> limo           0.24920    0.04605   5.412 2.76e-07 ***
> CO_gkg1        0.21015    0.03931   5.346 3.73e-07 ***
> C03Ca          0.01711    0.02727   0.628 0.531337
> ---
> Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
>
> Residual standard error: 6.249 on 135 degrees of freedom
>    (50 observations deleted due to missingness)
> Multiple R-squared: 0.9736,    Adjusted R-squared: 0.9014
> F-statistic: 13.47 on 370 and 135 DF,  p-value:<  2.2e-16
>
> So, in the desired linear model, arcilla should be just a line, with the valors of the linear model.
> I hope you understand better more. If not, I could make an english version of the file to send, so you can try the commands.
> Thanks a lot for your help!
>
>
>
> Arantzazu Blanco Bernardeau
> Dpto de Química Agrícola, Geología y Edafología
> Universidad de Murcia-Campus de Espinardo
>
>
[snip]



More information about the R-help mailing list