similar to: Plotting a curve for a Holling Type III Functional Response

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 function without error, so it's difficult for me, as a novice functioneer, to figure out

I'm getting an error message when using the new version of sqldf, > library(sqldf) > str(kdv) 'data.frame': 71 obs. of 3 variables: $ dpss: num 0.117 0.144 0.164 0.166 0.165 ... $ npdp: num 0.1264 0.0325 0.0109 0.0033 0.0055 ... $ logk: num 1.12 1.29 1.41 1.41 1.42 ... > test=sqldf("select * from kdv") Error in get("fun", env = this, inherits =

Hi, I am trying to run a generalized linear regression using a negative binomial error distribution. However, I want to use an overdispersion parameter that varies (dependent on the length of a stretch of road) so glm.nb will not do. >From what I've read I should be able to do this using GLM by specifying my own quasi family and describing the variance function using varfun, validmu,

Dear Ben, you answerd to Nancy Shackelford about Clarks 2Dt function. Since the thread ended just after your reply, I would like to ask, if you have an idea how to use this function in R I defined it the following way: function(x , p, u) { (p/(pi*u))*(1+(x^2/u))^(p+1) } and would like to fit this one to my obeservational data (count) [,1] [,2] [1,] 15 12 [2,] 45 13 [3,]

Hi I'm trying to replicate Smith et al.'s (http://www.sciencemag.org/content/330/6008/1216.abstract) findings by fitting their Gompertz and logistic models to their data (given in their supplement). I'm doing this as I want to then apply the equations to my own data. Try as a might, I can't quite replicate them. Any thoughts why are much appreciated. I've tried contacting the

Hi there, I have been having trouble running negative binomial regression (glm.nb) using library MASS in R v2.15.0 on Mac OSX. I am running multiple models on the variables influencing the group size of damselfish in coral reefs (count data). For total group size and two of my species, glm.nb is working great to deal with overdispersion in my count data. For two of my species, I am getting a

Hi Susan, In x86-64 the REX prefix must be used to access an extended register (r8-r15 and their aliases), but cannot be used when accessing the high byte of the ABCD regs (AH, BH, CH, DH). In your test case you have hardcoded %vreg1 to R8B, and %vreg15 to AH, and the test case contains a copy between these registers. The copy simultaneously must have a REX prefix, and cannot have a REX prefix,

Dear R-users: The following code produces two cones in two panels. What I would like to have is to have them in one, and to meet in the origin. Does anyone have any good ideas how to do this? Thanks for your help Jaakko library(lattice) A<-matrix(ncol=2, nrow=64) for(i in 0:63) { A[i+1,1]<-sin(i/10) A[i+1,2]<-cos(i/10) }

Does anyone have a way of producing solid shapes other than spheres in rgl? I am using rgl to produce a simple visualisation of a forest model results using "lollipops". Its just a bit of fun, but as many of the trees are pines I would like to depict their crowns as cones. If there is a solution I need it to work under windows. Here is the example. library(rgl) library(misc3d)

Dear r users, I have two kinds of data I'm trying to represent in Rgl. The first are measurements of soils: heights, diameters, and centerpoints of each. I would typically represent these in 3d space as cylinders or cones. The second are radar data of the area under the solids , xyz plus amplitude that I would like to represent as a point cloud. Is there a currently implemented method to

Hi, I am looking for a package which allows me to draw arrows (cylinders with cones on them, basically) in 3d images produced with rgl. Is such a thing available? I did notice the 'compositions' package, which contains 'arrow3D'. However, it relies on 'tensorA', which fails to install on my machine with the error message copied below. If anyone can help me fix this

Dear R experts, I have been encountering problems with the "optim" routine using "for" loops. I am determining the optimal parameters of several nested models by minimizing the negative Log-Likelihood (NLL) of a dataset. The aim is to find the model which best describes the data. To this end, I am simulating artificial data sets based on the model with the least number

Hi, We are preparing to make the release of OpenSSH 4.8 soon, so we would greatly appreciate testing of snapshot releases in as many environments and on as many operating systems as possible. The highlights of this release are: * Added chroot(2) support for sshd(8), controlled by a new option "ChrootDirectory". Please refer to sshd_config(5) for details, and please use this

On Fri, 4 Jun 2004, Makc wrote: > Here's what you can find in the net on YUV scheme: Hmm? What prompted this? > The statement I care to make here, is simply that there ain't neither > such thing as "luminance", which details "the eye is more sensitive to", "ain't neither"? Double negative, so I take your statement to mean that there *is* such

Well I suspect this list will pick up quite a bit since the article in Linux Magazine. So I'll try to get my questions in early. I looked at Asterisk about a year ago but the TPTB wanted a prebuilt solution (at that time) and I couldn't find anyone local enough to satisfiy them. We still don't have a new phone system because no one offered a good solution at a reasonable price. I

R 1.4.0 is now sinking into feature freeze, and we'd like some help testing it under Windows. Duncan Murdoch and Uwe Ligges have already volunteered. I've put a copy of SetupR140.exe at http://www.stats.ox.ac.uk/pub/bdr/RWin Also in that directory is a file TESTING suggesting some tests, not all of which have been run on this build. As the graphics internals have been changed a lot

Hello, I have a log likelihood function that I was able to optimize using mle2. I have two years of the data used to fit the function and I would like to fit both years simultaneously to test if the model parameter estimates differ between years, using likelihood ratio tests and AIC. Can anyone give advice on how to do this? My likelihood functions are long so I'll use the tadpole

Hello, I am interested in fitting a generalized nonlinear regression (gnlr) model with negative binomial errors. I have found Jim Lindsay's package that will do gnlr, but I have having trouble with the particular model I am interested in fitting. It is a threshhold model, where below a certain value of one of the parameters being fitted, the model changes. Here is a sample: Cones:

Hello, I want to know if there are some functions or packages to solve differential and integral equation using R. Thanks. Shao chunxuan. [[alternative HTML version deleted]]

Item Subject: cc:Mail Text Hi Mike, I got almost the same problem as you when I was using 2.04b on HP-UX 10.20 machine. I can not map to \\servername\sharename to my hpux machine from Win95 successfully when the sharename is different from the client login name. But it's successful on my WinNT. I can not resolve this problem but have a