model<-lm(weight ~ group,data=dataset)

by(dataset$weight,dataset$group,shapiro.test)
dataset$group: 1

    Shapiro-Wilk normality test

data:  dd[x, ]
W = 0.99387, p-value = 0.2839

-----------------------------------------------------------------
dataset$group: 2

    Shapiro-Wilk normality test

data:  dd[x, ]
W = 0.98957, p-value = 0.00698

-----------------------------------------------------------------
dataset$group: 3

    Shapiro-Wilk normality test

data:  dd[x, ]
W = 0.99635, p-value = 0.5562

-----------------------------------------------------------------
dataset$group: 4

    Shapiro-Wilk normality test

data:  dd[x, ]
W = 0.97931, p-value = 4.405e-05

-----------------------------------------------------------------
dataset$group: 5

    Shapiro-Wilk normality test

data:  dd[x, ]
W = 0.99338, p-value = 0.1075

-----------------------------------------------------------------
dataset$group: 6

    Shapiro-Wilk normality test

data:  dd[x, ]
W = 0.98036, p-value = 0.0287

-----------------------------------------------------------------
dataset$group: 7

    Shapiro-Wilk normality test

data:  dd[x, ]
W = 0.99308, p-value = 0.06546

-----------------------------------------------------------------
dataset$group: 8

    Shapiro-Wilk normality test

data:  dd[x, ]
W = 0.97932, p-value = 2.634e-05

-----------------------------------------------------------------
dataset$group: 9

    Shapiro-Wilk normality test

data:  dd[x, ]
W = 0.99666, p-value = 0.5945

Df   Sum Sq Mean Sq F value Pr(>F)
Group        8 13069322 1633665   97.77 <2e-16 ***
Residuals   3101 51816501   16710

Kruskal-Wallis chi-squared = 656.16, df = 8, p-value < 2.2e-16