# Re: [MORPHMET] help us in multifactor Manova

```My sample size is 147 not either 9 or 7 in any case. I agree that the
computer displays these numbers! don't know why?```
```
Nevertheless, in this analysis the number of specimens (dependent shape
variable) is 147, obviously, with the equal number of observations for each
independent variables (here, we used two factors ind\$Stages and ind\$Location
).

It is interesting to note that I am getting good results using single
factor Manova; dependent shape variable over, either  ind\$stages or
ind\$location, at a given time,  though with/without interactions as well. Even,
it works fine with more factors if I model without any effect of
interaction, and the results are given below

> Paint<-procD.lm(coords ~ Csize + ind\$Stages + ind\$Location, data =
> Paintgeo147, iter = 9999, RRPP = FALSE, print.progress = FALSE, int.first =
> FALSE)

> Paint3factor

Call:

procD.lm(f1 = coords ~ Csize + ind\$Stages + ind\$Location, iter = 9999,

RRPP = FALSE, int.first = FALSE, data = Paintgeo147,
print.progress = FALSE)

Type I (Sequential) Sums of Squares and Cross-products

Randomization of Raw Values used

10000 Permutations

Df      SS       MS      Rsq       F      Z Pr(>F)

Csize          1 0.04688 0.046880 0.064257 13.4431 7.9867  1e-04 ***

ind\$Stages     4 0.17575 0.043938 0.240899 12.5994 8.3726  1e-04 ***

ind\$Location   1 0.01871 0.018713 0.025650  5.3662 3.2145 0.0042 **

Residuals    140 0.48822 0.003487

Total        146 0.72956

---

Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

The main issue, here, is that somehow I don’t know how to use the proper
command to quantify effects of interaction, say for example ind\$Stages:
ind\$Location, Csize: ind\$Stages, Csize: ind\$Locations... etc. ????

Thanks for spending your time and effort!

Looking forward

Truly

Mahen

On Fri, May 18, 2018 at 12:05 AM, Mike Collyer <mlcoll...@gmail.com> wrote:

> OK, this error is indicating you have variables of different size - one
> with 9 and one with 7.  This does not appear to be because of the int.first
> issue, but rather that you introduced a new variable with a different
> number of observations.  Also, I’m not sure what kind of result you expect
> with such a complex model and so few observations.  It might not work,
> anyway, because the number of independent variables exceeds the number of
> observations.  Just a warning.
>
> Mike
>
>
> On May 17, 2018, at 2:19 PM, mahendiran mylswamy <mahenr...@gmail.com>
> wrote:
>
> Hi Mike
>
>
> Thanks indeed for the quick response.
> error messages are given below
>
>
> > PSfact<-procD.lm(coords ~ Csize * ind\$Stages * ind\$Location, data =
> rough, iter = 99, RRPP = FALSE, int.first = TRUE, print.progress = FALSE)
> Error in data.frame(df, SS, MS, Rsq = R2, F = Fs) :
>   arguments imply differing number of rows: 9, 7
> In SS/df : longer object length is not a multiple of shorter object length
>
>
>
>
>
>
> On Thu, May 17, 2018 at 6:22 PM, Mike Collyer <mlcoll...@gmail.com> wrote:
>
>> Mahen,
>>
>> Can you be specific with the error statement?  It might not me related to
>> the int.first use but more so how you set up your formula.  If ind\$Stages
>> and ind\$Location are in the data frame, you are asking R to create a model
>> matrix from dat in an environment within another environment.  This might
>> cause some problems.
>>
>> Please let us know what the error says, and we can see if it is obvious.
>>
>> Cheers!
>> Mike
>>
>> On May 17, 2018, at 8:22 AM, mahendiran mylswamy <mahenr...@gmail.com>
>> wrote:
>>
>> Dear All.
>>
>> I did a single and multi-factor Manova, through Geomorph package in R,
>> using the following commands given below, it is working fine!
>>
>> *Paint1factor <-procD.lm (coords ~ Csize*indep1\$Stages, data =
>> Paintgeo147, iter = 9999, RRPP = FALSE, print.progress = FALSE)*
>>
>> *Paint3factor<-procD.lm(coords ~ Csize + ind\$Stages + ind\$Location, data
>> = Paintgeo3fac, iter = 99, RRPP = FALSE, print.progress = FALSE, int.first
>> = FALSE)*
>>
>> HOWEVER,  when I did the multifactor Manova  to check the interaction
>> effects, using the command, though i inserted int.first = false,  somehow
>> it’s not working often shows error
>>
>> *Paint3factor<-procD.lm(coords ~ Csize * ind\$Stages * ind\$Location, data
>> = Paintgeo3fac, iter = 99, RRPP = FALSE, print.progress = FALSE, int.first
>> = FALSE)*
>>
>>
>> apart from changing, + to *, inserting int.first, any other stuff to be
>>
>>
>>
>> Truly
>>
>> Mahen
>>
>>
>> --
>> ***************************************
>> M Mahendiran, Ph D
>> Scientist - Division of Wetland Ecology
>> Salim Ali Centre for Ornithology and Natural History (SACON)
>> Anaikatti (PO), Coimbatore - 641108, TamilNadu, India
>> Tel: 0422-2203100 (Ext. 122), 2203122 (Direct), Mob: 09787320901
>> Fax: 0422-2657088
>> http://www.sacon.in/staff/dr-m-mahendiran/
>>
>> P Please consider the environment before printing this email
>>
>> --
>> MORPHMET may be accessed via its webpage at http://www.morphometrics.org
>> ---
>> You received this message because you are subscribed to the Google Groups
>> "MORPHMET" group.
>> To unsubscribe from this group and stop receiving emails from it, send an
>> email to morphmet+unsubscr...@morphometrics.org.
>>
>>
>>
>
>
> --
> ***************************************
> M Mahendiran, Ph D
> Scientist - Division of Wetland Ecology
> Salim Ali Centre for Ornithology and Natural History (SACON)
> Anaikatti (PO), Coimbatore - 641108, TamilNadu, India
> Tel: 0422-2203100 (Ext. 122), 2203122 (Direct), Mob: 09787320901
> Fax: 0422-2657088
> http://www.sacon.in/staff/dr-m-mahendiran/
>
> P Please consider the environment before printing this email
>
>
>

--
***************************************
M Mahendiran, Ph D
Scientist - Division of Wetland Ecology
Salim Ali Centre for Ornithology and Natural History (SACON)
Anaikatti (PO), Coimbatore - 641108, TamilNadu, India
Tel: 0422-2203100 (Ext. 122), 2203122 (Direct), Mob: 09787320901
Fax: 0422-2657088
http://www.sacon.in/staff/dr-m-mahendiran/

P Please consider the environment before printing this email

--
MORPHMET may be accessed via its webpage at http://www.morphometrics.org
---
You received this message because you are subscribed to the Google Groups
"MORPHMET" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to morphmet+unsubscr...@morphometrics.org.
```