R을 통한 ANOVA 분석에 대해 질문드립니다.

R에 대해 정확히 one-way ANOVA까지밖에 , 그것도 정말 얕게 아는데, 갑자기 ANCOVA라는 기법까지 하게 되었습니다.

(원리는 모르지만)
A <-aov(A~group, data=asdf)
anova(A) 를 하면 그룹 간의 분산이 그룹 사이의 차이 때문에 나온거다... 이 정도는 압니다.

ANCOVA라는 방식은 여러 개의 독립 변인과 하나의 종속 변인이 있을 때, 그 중 하나의 독립 변인을 주요 변인으로 삼고, 다른 변인은 통제 변인으로 두는 방식이라고 하더군요. 여기서 two-way ANOVA라는 친구는 두 개의 독립 변인을 모두 주요 변인으로 삼는 경우라고 들었습니다.

그래서 구글링을 많이 해보니
제가 만약 y이라는 변수에 대해서 a, b, c, d 라는 네 가지 변수를 독립 변수로 놓고 a, b, c를 통제 변인으로 둔 뒤 범주병 변수 d에 대해서 d1 d2 d3 세 그룹에 대해 ANCOVA를 돌리고 싶습니다.
이때 x1 = aov(y~d) 를 하고 anova( x1 )를 하면 제가 아는 anova 분석 결과가 나오지요.
x2=lm(y~a+d)를 하면 a를 공변수로 두고 anova(x2)를 해서 d에 대한 분산분석을 수행하는 ancova가 나온다고 구글링을 통해 배웠습니다.

그렇다면 x3=lm(y~a+b+c+d)를 하고 anova(x3)를 한 수행한다면 자동으로 d1 d2 d3에 대한 ancova를 수행했다고 볼 수 있는건가요?
총체적으로 R, anova()함수, 분석에 대해 아는게 없다보니 힘드네요ㅠㅠ 내공 백만점입니다.
R에 대한 책이나 사이트도 함께 추천해주시면 감사하겠습니다!