[R] convhulln {geometry} output from .call

2007-07-13 Thread Daniel E. Bunker
t completed. Verifying that all points are below 6.9e-15 of all facets. Will make 100 distance computations. > -- Daniel E. Bunker TraitNet Associate Director BioMERGE Associate Director Post-Doctoral Research Scientist Columbia University Department of Ecology, Evolution and Environmental Biolo

[R] aovlmer.fnc, lmer and binomial response

2007-04-24 Thread Daniel E. Bunker
take it this means that it will not work with family="binomial"? If so, how can I modify it to assess binomial response data? Any advice or guidance would be greatly appreciated! -Dan -- Daniel E. Bunker BioMERGE Associate Director Post-Doctoral Research Scientist Columbia

Re: [R] line width (all elements) in Trellis

2006-11-06 Thread Daniel E. Bunker
Great! Thanks to both Paul and Deepayan for these efficient approaches! -Dan Deepayan Sarkar wrote: > On 11/6/06, Paul Murrell <[EMAIL PROTECTED]> wrote: >> Hi >> >> >> Deepayan Sarkar wrote: >> > On 11/6/06, Daniel E. Bunker <[EMAIL PROTECTED]>

[R] line width (all elements) in Trellis

2006-11-06 Thread Daniel E. Bunker
d=c(rep(2,7)) trellis.par.set("reference.line", tpar12) barchart(yield ~ variety | site, data = barley, groups = year, layout = c(1, 6), ylab = "Barley Yield (bushels/acre)", scales = list(x = list(abbreviate = TRUE, minlength = 5))) -- Daniel E. B

Re: [R] adding error bars to lattice plots

2006-10-13 Thread Daniel E. Bunker
uperpose? Again, any advice would be greatly appreciated. Thanks, Dan Deepayan Sarkar wrote: > On 10/12/06, Daniel E. Bunker <[EMAIL PROTECTED]> wrote: >> Dear R users, >> >> About a year ago Deepayan offered a suggestion to incorporate error bars >> into a dotplot us

[R] adding error bars to lattice plots

2006-10-12 Thread Daniel E. Bunker
Dear R users, About a year ago Deepayan offered a suggestion to incorporate error bars into a dotplot using the singer data as an example <>. When I try to utilize this code with a grouping variable, I get an error stating that the su

[R] Convergence warnings from zeroinfl (package pscl)

2006-07-20 Thread Daniel E. Bunker
e-04 modacpas:intact -0.009994 0.002735 -3.6546 2.576e-04 log(theta) -1.405175 0.062412 -22.5145 2.993e-112 Theta = 0.2453 -- Daniel E. Bunker BioMERGE Associate Director Post-Doctoral Research Scientist Columbia University Department of Ecology, Evolution and Environmental Bio

[R] xYplot, lmline and abline

2006-06-06 Thread Daniel E. Bunker
tes, pch=c(1,2,3,4), data=data1, panel = function(x, y) { panel.xyplot(x,y) panel.lmline(x,y, span=1) } ) -- Daniel E. Bunker BioMERGE Associate Director Post-Doctoral Research Scientist Columbia University Department of Ecology, Evolution and En

[R] vecortizing uniroot() for numerical solutions

2005-06-30 Thread Daniel E. Bunker
#but R rejects the multiple values for d2. #I could use a numerical solver such as XPPAUT, but would rather be able to do it all in R. #Any thoughts on vectorizing this would be greatly appreciated. #Thanks! -Dan -- Daniel E. Bunker Associate Coordinator - BioMERGE Post-Do

[R] vectorized approach to cumulative sampling

2005-04-07 Thread Daniel E. Bunker
Dan # loop approach old=c(1:10) p=runif(1:10) target=20 newsum=0 new=NULL while (newsumtarget){new[length(new)]=target-sum(new[-length(new)])} new newsum=sum(new); newsum target -- Daniel E. Bunker Associate Coordinator - BioMERGE Post-Doctoral Research Scientist Columbia University Department of Ecol

[R] Efficient selection and alteration of dataframe records

2005-02-03 Thread Daniel E. Bunker
umdflength=length(sumdf$s.n) sumdf[sumdflength+1,1]=n+1; sumdf[sumdflength+1,2]=S; sumdf[sumdflength+1,3]=sumcrop; }#d. REPEAT SPECIES DELETION LOOP #housekeeping rm(delsp, plotdf, plotdf.del, plotdf.old, plotdft, traitdf.del) gc() #8. plot results, fit line