search for: harrup

...are listed there are 4 semio categories, 1,2,3 and 4 but only 2,3,and 4 are listed in the summary (semio2,semio3,semio4). And I cant for the life of me work out why category one (semio1) is being ignored, missing etc. Any help or suggestions would be most appreciated. Thanks in advance Lara lara.harrup at bbsrc.ac.uk Call: glm(formula = total ~ (temperature + humidity + semio)^2, family = poisson) Deviance Residuals: Min 1Q Median 3Q Max -22.212 -5.132 -2.484 3.200 18.793 Coefficients: Estimate Std. Error z value Pr(>|z|) (Inter...

...d and if there are any other accuracy measure that can be used with Random Forest in addition to the pseudo r^2 value? as this work will be presented to an entomology/ecology audience where machine learning is a bit outside their (and my) statistics comfort zone. Many thanks in advance Lara lara.harrup@bbsrc.ac.uk [[alternative HTML version deleted]]

Hi I am trying to explore the use of random forests for regression to identify the important environmental/microclimate variables involved in predicting the abundance of a species in different habitats, there are approx 40 variable and between 200 and 500 data points depending on the dataset. I have successfully used the randomForest package to conduct the analysis and looked at the %IncMSE

...otstrap them to find out more about how they perform. Or am I going about this all the wrong way? is there away to calulate vifs from regsubsets or pass it directly to something that calculates them e.g. VIF in the car package? Any help will be most appreciated, many thanks in advance Lara lara.harrup at bbsrc.ac.uk

...[1] 0.8216521 $auc.app [1] 0.8455208 $auc.samples list() my question is how do I access the $auc.samples list() or write it to a csv file or something similar to get the 2.5% and 96.5% quartile range from? or is there another way to get the 95% CI for the AUC? Many thanks in advance Lara Lara.harrup@bbsrc.ac.uk ############ library(Daim) mylda <- function(formula,train,test){ model <- lda(formula,train) predict(model,test)$posterior[,"pos"] } x<-Daim(OUTBREAK2 ~ COHESION.2 + COHESION.23 + ED.2 + PLAND.12 + PLAND.2 + PLAND.25 + PLAND.26 + ALTITUDE_MEAN + SLOPE_MEAN +...

...S' Not sure if this is a general problem with loading an S dialect package and I need to add some thing else before loading the library if it is in S rather than R? or if it is a problem with the 'glar' package. Any suggestions/tips most appreciated, thanks in advance. Lara lara.harrup at bbsrc.ac.uk

Hi I have some experimental data where I have counts of the number of insects collected to different trap types rotated through 5 different location (variable -location), 4 different chemical attractants [A, B, C, D] were applied to the traps (variable - semio) and all were trialled at two different CO2 release rates [1, 2] (variable CO2) I also have a selection of continuous variables

Hi I have been using a Zero-Inflated negative binomial model fitted using the pscl zeroinfl command but I would like to extract a goodness of fit measure are there any suitable pseudo R^2 measures available for this type of analysis to try and assess the amount of variation in the data explained by the model? I have tried with the pR2 command in pscl (for computing various pseudo R2

Hi all I have developed a zero-inflated negative binomial model using the zeroinfl function from the pscl package, which I have carried out model selection based on AIC and have used likelihood ratio tests (lrtest from the lmtest package) to compare the nested models [My end model contains 2 factors and 4 continuous variables in the count model plus one continuous variable in the zero-inflated