Displaying 3 results from an estimated 3 matches for "prof1".
Did you mean:
prof
2007 Nov 30
2
nls and discontinuous functions
...e and could not find a solution; therefore I am sending you this email seeking for help. I need to solve a discontinuous nonlinear function, but R 2.6 keeps giving me an error. Below is a simple example:
f <- function (t,k,l,u) { r0<-346.1 x<-t-l if (x<0) r0 else (r0-u)*exp(-k*x)+u }prof1<-read.table("c:\\data.txt", header=TRUE)results<-nls(y1~f1(t,k,l,u), data=prof1, start=list(k=0.05, l=4, u=154))
Errors I am getting are:
Error in qr(.swts * attr(rhs, "gradient")) : dims [product 3] do not match the length of object [10]In addition: Warning messa...
2009 Jul 30
0
temporal and spatial pseudorreplication with lmer
...urements (3 seasons during 3 years and 2 season during the last year, factor named percod, with 11 levels) of a bird presence/absence (pres.f) on 14 artificial ponds (charca.f). The ponds are integrated in 4 greater spatial units (zepa). I have 11 explanatory variables, some of them are continuous (prof1,prof2,prof4,capac,perim_veget,div_veg,div_herp) and some are factors (ungul.f,ganado.f,turb.f,peces.f). Not all artificial ponds have been monitored in all seasons, so I have some missing values in the dependent variable.
I have thought that the correct model should be, ignoring the spatial pseu...
2008 Nov 30
1
methods not found inside function?
...ts: object
Use showMethods("summary") for currently available ones.
summary(gm1)
## Generalized linear mixed model fit by the Laplace approximation
## Formula: cbind(incidence, size - incidence) ~ period + (1 | herd)
## Data: cbpp
## etc ...
library(profileModel)
debug(profileModel)
prof1.gm1 <- profileModel(fitted=gm1,
objective=deviance)
## step into the function, and then try:
Browse[1]> summary
function (object, ...)
UseMethod("summary")
<environment: namespace:base>
Browse[1]> summary(gm1)
Length Class Mode
1 mer...