Re: [R] Trouble embedding functions (e.g., deltaMethod) in other functions

Ah Yes, good point. However, it still does not correct the situation. I still get NaN's. -- View this message in context: http://r.789695.n4.nabble.com/Trouble-embedding-functions-e-g-deltaMethod-in-other-functions-tp4662178p4662187.html Sent from the R help mailing list archive at Nabble.com.

[R] Trouble embedding functions (e.g., deltaMethod) in other functions

Dear R community, I've been writing simple functions for the past year and half and have come across a similar problem several times. The execution of a function within my own function produces NaN's or fails to execute as intended. My conundrum is that I can execute the function outside of my fu

Re: [R] faulty F-tests

Thanks so much. I see my foolish ways now. -- View this message in context: http://r.789695.n4.nabble.com/aov-and-anova-making-faulty-F-tests-tp4660407p4660439.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org mail

[R] aov() and anova() making faulty F-tests

Dear useRs, I've just encountered a serious problem involving the F-test being carried out in aov() and anova(). In the provided example, aov() is not making the correct F-test for an hypothesis involving the expected mean square (EMS) of a factor divided by the EMS of another factor (i.e., inste

[R] RGL axis labels parallel to axes

Hello useRs, I'm having trouble trying to produce axis labels running parallel to the axes in 3d space. Is their a way to do this when producing 3d graphs using the rgl package? Thanks for your help, Patrick -- View this message in context: http://r.789695.n4.nabble.com/RGL-axis-labels-paral

Re: [R] RGL 3D curvilinear shapes

;), lit = F, back = "line") lines3d(x1, y1, z1) lines3d(x2, y2, z2) lines3d(x3, y3, z3) triangles3d(x= t1[,1], y= yfit1, z= t1[,2], lit = F, col = "white", front = "fill", lwd = 10) triangles3d(x= t2[,1], y= yfit2, z= t2[,2], lit = F, col = "white", front =

[R] RGL 3D curvilinear shapes

Dear useRs, I'm trying to simply fill in the area under a curve using RGL. Here' the set up: x <- c(0.75,75.75,150.75,225.75,300.75,375.75,450.75,525.75,600.75,675.75, 0.5,50.5,100.5,150.5,200.5,250.5,300.5,350.5,400.5,450.5, 0.25,25.25,50.25,75.25,100.25,125.25,150.25,175.25,200.2

[R] pmodels in DRC

Dear R users, I'm a little lost on how to define pmodels for the DRC package. My goals are to produce isoboles of binary toxicity data. any tips? I really just need to know what pmodels refers to. Cheers, Pat -- View this message in context: http://r.789695.n4.nabble.com/pmodels-in-DRC-tp4

Re: [R] drc results for binomial data

P.S. Here are the results I get for; #---data 1--# Model fitted: Generalised log-logistic (log(ED50) as parameter) (3 parms) Parameter estimates: Estimate Std. Errort-value p-value Hillslope:(Intercept) -6.05453 2.1

Re: [R] drc results for binomial data

I've been having trouble making sense of the drc results for my binomial response toxicity data. Firstly, the standard errors are far too large for how well the data fit the log-logistic model, particularly compared to other methods of LC50 estimation (e.g. probit and trimmed spearman-karber). Seco

[R] Trimmed Spearman-Karber

Hello All, I am searching for a package which allows me to run an LC50 (LD50) estimate using the trimmer spearman-karber method. Has any one developed code for this analysis? In my case, probit or logit analyses are less accurate. Please help! I will post code if I end up finding/writing it.