similar to: which model (GLMs)is the best?

Dear friends, summary() doesn't give a good ouput style,e.g.: grass soiltem airtem gheight humidity altitude diluo :38 Min. :15.90 Min. :17.70 Min. : 8.00 Min. : 0.2360 high: 43 huanghuacai:32 1st Qu.:19.32 1st Qu.:22.60 1st Qu.:40.00 1st Qu.: 0.3190 low :119 hucao :46 Median :20.20 Median :25.30 Median :60.00

Dear friends, I used R and SAS to analyze my data through generalized linear model, and there is some difference between them. Results from R: glm(formula = snail ~ grass + gheight + humidity + altitude + soiltemr + airtemr, family = Gamma) Deviance Residuals: Min 1Q Median 3Q Max -1.23873 -0.41123 -0.08703 0.24339 1.21435 Coefficients:

Dear Rusers, My question is about "recode variables". First, i'd like to say something about the idea of recoding: My dataset have three variables:type,soiltem and airtem,which means grass type, soil temperature and air temperature. As we all known, the change of air temperature is greater than soil temperature,so the values in those two different temperaturemay represent different

Dear friends, when i use glm() to fit my data, i use glm(formula = snail ~ vegtype + mhveg + humidity + elevation + soiltem, *family = Gamma(link = inverse),* data =a,)) It shows: error in eval(expr, envir, enclos) : *gamma distribution don't allow negative value*. But i use result<-glm(formula = snail ~ vegtype + mhveg + humidity + elevation + soiltem, family = poisson, data =a) #this

Hi all. I have the problem, that the labels from a factor variable in a tree object, (more exact, the node - labels) is labeled in a plot with letters a,b,c... instead of the numerical levels. when I run >labels(tree.object) It gives me the correct labels (the problem concerns 'landuse': [1] "root" "landuse:1,2,3,4,5,6,7,8,9,10" [3]

Try something like with(df, predict(smooth.spline(x = altitude, y = atm_values), deriv = 1)) Cheers, Andrew -- Andrew Robinson Chief Executive Officer, CEBRA and Professor of Biosecurity, School/s of BioSciences and Mathematics & Statistics University of Melbourne, VIC 3010 Australia Tel: (+61) 0403 138 955 Email: apro at unimelb.edu.au Website: https://researchers.ms.unimelb.edu.au/~apro

Hi everyone, I have a vector with atmospheric measurements (x-axis) that is obtained/calculated at different altitudes (y-axis). The altitude is uniformly distributed every 7 meters. For example my dataframe is: df <- dataframe( *altitude* = c(1005, 1012, 1019, 1026, 1033, 1040, 1047, 1054, 1061, 1068), *atm_values* = c(1.41, 1.40, 1.39, 1.38, 1.37, 1.37, 1.38, 1.36, 1.33, 1.31)

Hi Konstantinos Not exactly derivative but > diff(df[,2]) [1] -0.01 -0.01 -0.01 -0.01 0.00 0.01 -0.02 -0.03 -0.02 May be enaough for you. Cheers Petr > > -----Original Message----- > From: R-help <r-help-bounces at r-project.org> On Behalf Of konstantinos > christodoulou > Sent: Tuesday, January 31, 2023 10:16 AM > To: r-help mailing list <r-help at

Hello, I'm working on 30 natural populations of Aster amellus L. a threatened plant species. My aim is to see if the population size influence on several fitness trait. I build a multi-factorial ANOVA. The independents variables are altitude, humidity, soil component, maintenance .... and population density. Unfortunately my data have some heteroscedasticity even after classical

Hello R-Users, I have recently run into several problems using vgam() in the VGAM package. I am hoping someone might have some solutions... Briefly, I have been trying to fit GAM models for zero-altered negative binomial models. 1. When fitting smoothed parameters (e.g. s(X, df=2)) changing the degrees-of-freedom has no effect on the level of smoothing (e.g. number of knots for the

Dear listers, I have written a function to facilitate the drawing of altitude profiles with x (distance), y (altitude) and a z parameter (altitude magnification). profplot<-function(x,y,z=10,...){ op <- par()$mai par(mai=c(0.95625,0.76875,0.76875,0.95625)) plot(x,y*z, type="l",asp=1,las=1,xlab="",ylab="",yaxt="n",...)

Hi, The script below is my current coding in order to produce a contour plot of temperature across altitude and time. In my case,?time,altitude and temperature?are represented by x, y and z variables. ############################################## Brazilan.Pallete <- colorRampPalette(c("blue","green","yellow","red")) image.plot(x, y, z, col =

I tried the following ordered logistic regression in R: mod1 <- polr(altitude~sp + wind_dir + wind_speed + hr, data=altioot) But when I asked The summary of my regression I got the folloing error message: > summary (mod1) Re-fitting to get Hessian Error in optim(start, fmin, gmin, method = "BFGS", hessian = Hess, ...) : the initial value of 'vmin' is not

I have drawed a picture with persp, it's 3d map with different color, indicate different altitude. In gnuplot, the corresponding command 'splot' will generate a picture beside to indicate the relationship between color and altitude. But in R, how to draw it? I have read the manual of legend, but they are all about how to draw a legend with colored text, not a continuous varing color

Hi i need some help with this exercise: FIles: https://mega.nz/#!JxMFGIwC!qA85SBIBRVagCzYfmLwSvGuNK_qXqCXrakPxXryCpGg #PARZIAL 3: GEO #Data: # Shapefile "INCOME" contains dummy information about revenue # Common Abbreviations in the "INCOME" variable and the centroid altitude #dell common in the variable "ALT" #Richieste # 1 #map of the variable

Hi R community I need to interpolate precipitation data for a natural park. I have precipitation data from some climate stationts. (I know the table is not complete but I only need to show you X,Y, Altitude and PrepJul) X Y Altitude PrepJan PrepFeb PrepMar PrepAp PrepMay PrepJun PrepJul 597706 4093438 41 0 9 77,8 63,1 17,5 0 2,6 597535

Hi all, I'm looking for to interpolate hourly temperature date collected from more than 140 automatic weather station (irregularly spaced) using 4 independent variable: 1-2) geografic coordinates (x,y) (from DEM - 40m) 3) altitude (z) (from DEM - 40m) 4) solar radiation (from a model calculated with grass: r.sun) In addition, I would like to use also "time" variable (e.g.: hours).

Hi All, I have an asterisk box on my DMZ, and I'm using a PF for my firewall, I can make a call but some reasons I have a dead air. Any Ideas? below are my rules... ext_if = "bce0" int_if = "bce1" altitude = "172.16.1.0/24" #### machines #### vbox = "172.16.1.1" uci = "172.16.1.4" voices = "203.172.x.1" ipc =

Hi all, this is my first post in this mailing group. I hope that anyboby could help me parsing a xml file. I found this website http://www.omegahat.org/RSXML/gettingStarted.html but unfortunately my XML file is not as easy as the one in the example. Example: <?xml version="1.0" encoding="UTF-8"?> <?xml-stylesheet

Dear r-helpers, I use function step.gam (package gam, T. Hastie) with several explanatory variables to build a model. Unfortunately, I obviously have too many variables. This message occurs on my 4 core 64bit machine with 8GB RAM in R2.11.1 for Windows (64bit build): Error in array(FALSE, term.lengths) : 'dim' specifies too large an array I read that this message occurs when running out