Re: [R] ANCOVA error again

From: Birgit Lemcke <birgit.lemcke_at_systbot.uzh.ch>
Date: Mon, 21 Apr 2008 15:43:03 +0200

Hello R users!

I got again an error message.

I used this code:

	with (FemMal85_Sex, {
							 ModelFemMal85<-				glm 
(Sex~outLatTep_like_other*outLatTep_like_conduplicate*outLatTep_keeled_w 
inged*spathellae_conspicuous*spathellae_inconspicuous_absent             
   	 

*InfSpath_persistence*InfSpath_caducuous*bractsSpacing_lax*bractsSpacing _imbricate*InfType_sparsely_paniculate*InfType_racemose*InfType_panicula

te*InfType_globose*bracApexShape_truncate                       
*bracApexShape_rounded                       
*bracApexShape_obtuse                       
*bracApexShape_acute                       
*bracApexShape_acuminate                   
*bracApexShape_apiculate                 
*bracApexShape_aciculate                   
*BracUpperMarg_like_rest*BracUpperMarg_memebranous*BracUpperMarg_honeyco mbed_cells*InfSpathText_coriaceous*InfSpathText_hyaline*InfSpathText_cha rtaceous*InfSpathText_cartilaginous*InfSpathText_membranous*spikShapeSid e_linear*spikShapeSide_oblong*spikShapeSide_square*spikShapeSide_ellipti cal*spikShapeSide_ovate*spikShapeSide_obovate*spikShapeSide_obtriangular *spikShapeSide_orbicular*spikShapeSide_undifferentiated*SpikApexShape_tr
uncate*SpikApexShape_rounded*SpikApexShape_obtuse*SpikApexShape_acute*Sp 
ikApexShape_undifferentiated*BracShape_linear*BracShape_oblong*BracShape 
_square*BracShape_elliptical*BracShape_ovate*BracShape_obovate*BracShape 
_orbicular*BracText_bony*BracText_coriaceous*BracText_hyline*BracText_ch 
artaceous*BracText_cartilaginous                         
*BracText_membranous                         
*BracText_centrChartaceousMargMembranous               
*TepText_bony*TepText_coriaceous*TepText_chartaceous                     
       *TepText_cartilaginous  

*TepText_membranous*InfLengthMin*InfLengthMax*InfWidthMin*InfWidthMax*Sp athellaeLengthMin*SpathellaeLengthMax*SpikLengthMin*SpikLengthMax*FlowNu
mbSpikMin*FlowNumbSpikMax*BracLengthMin*BracLengthMax*FlowLengthMin*Flow 
LengthMax*InfSpathLengthToSpikMin*InfSpathLengthToSpikMax*TepInOutMin*Te 
pInOutMax*BracLengthtoFlowMin*BracLengthtoFlowMax*BracMargMin*BracMargMa 
x*BracAwnToBodyMin*BracAwnToBodyMax,
na.action=na.exclude,family=binomial)})

and got this error message:

Traceback:

  1. terms.formula(formula, data = data)
  2. terms(formula, data = data)
  3. model.frame.default(formula = Sex ~ outLatTep_like_other * outLatTep_like_conduplicate *........... * BracAwnToBodyMax, drop.unused.levels = TRUE)
  4. model.frame(formula = Sex ~ outLatTep_like_other * outLatTep_like_conduplicate *........... * BracAwnToBodyMax, drop.unused.levels = TRUE)
  5. eval(expr, envir, enclos)
  6. eval(mf, parent.frame())
  7. glm(Sex ~ outLatTep_like_other * outLatTep_like_conduplicate *............* BracAwnToBodyMax, family = binomial)
  8. eval.with.vis(expr, envir, enclos)
  9. eval.with.vis(ei, envir)
  10. source("/Users/birgitlemcke/Job/Doktorarbeit/R/Protokolle_Codes/ Protokoll21.04.08.R")

Possible actions:

1: abort (with core dump, if enabled)
2: normal R exit
3: exit R without saving workspace
4: exit R saving workspace

Selection:

........... I deleted here some of the 85 variables

What does this message mean?

Thanks a lot in advance.

B.

Am 21.04.2008 um 14:50 schrieb John Fox:
> Dear Brigit,

>

> My guess is that you forgot to specify the argument family=binomial in
> the call to glm().
>

> Had you included the commands that you used as well as the error that
> was produced, it wouldn't be necessary to guess.
>

> I hope this helps,
> John
>

> On Mon, 21 Apr 2008 14:23:13 +0200
> Birgit Lemcke <birgit.lemcke_at_systbot.uzh.ch> wrote:
>> R version 2.6.2 PowerBook G4
>>
>> Hello R User,
>>
>> I try to perform an ANCOVA using the glm function.
>> I have a dataset with continuous and categorical data (explanatory
>> variables) and my response variable is also binary categorical.
>>
>> Fehler: NA/NaN/Inf in externem Funktionsaufruf (arg 4)
>> Zusätzlich: Warning messages:
>> 1: In Ops.factor(y, mu) : - nicht sinnvoll für Faktoren (makes no
>> sense for factors)
>> 2: In Ops.factor(eta, offset) : - nicht sinnvoll für Faktoren
>> 3: In Ops.factor(y, mu) : - nicht sinnvoll für Faktoren
>>
>> My dataset contains NA`s but if I try to use na.exclude, I got the
>> same Error message.
>>
>> I thought the function should use with my dataset. What am I doing
>> wrong?
>>
>> Thanks in advance for your help.
>>
>> Birgit
>>
>>
>> Birgit Lemcke
>> Institut für Systematische Botanik
>> Zollikerstrasse 107
>> CH-8008 Zürich
>> Switzerland
>> Ph: +41 (0)44 634 8351
>> birgit.lemcke_at_systbot.uzh.ch
>>
>> 175 Jahre UZH
>> «staunen.erleben.begreifen. Naturwissenschaft zum Anfassen.»
>> MNF-Jubiläumsevent für gross und klein.
>> 19. April 2008, 10.00 Uhr bis 02.00 Uhr
>> Campus Irchel, Winterthurerstrasse 190, 8057 Zürich
>> Weitere Informationen http://www.175jahre.uzh.ch/naturwissenschaft
>>
>> ______________________________________________
>> 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.
>

> --------------------------------
> John Fox, Professor
> Department of Sociology
> McMaster University
> Hamilton, Ontario, Canada
> http://socserv.mcmaster.ca/jfox/

Birgit Lemcke
Institut für Systematische Botanik
Zollikerstrasse 107
CH-8008 Zürich
Switzerland
Ph: +41 (0)44 634 8351
birgit.lemcke_at_systbot.uzh.ch

175 Jahre UZH
«staunen.erleben.begreifen. Naturwissenschaft zum Anfassen.» MNF-Jubiläumsevent für gross und klein.
19. April 2008, 10.00 Uhr bis 02.00 Uhr
Campus Irchel, Winterthurerstrasse 190, 8057 Zürich Weitere Informationen http://www.175jahre.uzh.ch/naturwissenschaft



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. Received on Mon 21 Apr 2008 - 13:49:16 GMT

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Mon 21 Apr 2008 - 15:30:30 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.

list of date sections of archive