similar to: tree mortality risk model using glm()

Dear all, As part of a larger function, I am randomly removing rows from a data frame. The number of removed rows is determmined by a Poisson distribution with a low mean. Sometimes, the random number is 0, and that's when the problem starts: My data frame: > temp occ x y dbh age 801 0 2977.196 3090.225 6 36.0 802 0 2951.892 3083.769 8 40.6 803 0 2919.111

Dear all, Last autumn there was some discussion on the list of the warning Warning message: fitted probabilities numerically 0 or 1 occurred in: (if (is.empty.model(mt)) glm.fit.null else glm.fit)(x = X, y = Y, when fitting binomial GLMs with many 0 and few 1. Parts of replies: "You should be able to tell which coefficients are infinite -- the coefficients and their standard errors will

To whom it may concern, I get the following message when I run my function: Warning message: multi-argument returns are deprecated in: return(call.fn, repl, time, from, to, last.year, occup.m, ant.occ.m, > version platform i386-pc-mingw32 arch i386 os mingw32 system i386, mingw32 status major 1 minor 8.0

Dear Lockwood, As you can see, I'm a beginner... But thank you very much! Sincerely, Tord Quoting "J.R. Lockwood" <lockwood at rand.org>: > how about > > buskartant$buskartant <- sapply( group.list, function(x) > sum(!is.na(unique(x))) ) > > OR > > buskartant$buskartant <- sapply( group.list, function(x) > length(unique(x[!is.na(x)])) )

Dear all, Is it possible to somshow plot text as italic AND bold. I tried font=c(2,3) in text(), but it doesn't work. It seems like the latter value is used. Thanks in advance! Sincerely, Tord ----------------------------------------------------------------------- Tord Sn?ll Avd. f v?xtekologi, Evolutionsbiologiskt centrum, Uppsala universitet Dept. of Plant Ecology, Evolutionary Biology

Dear all, I want my PhD thesis which I hand in tomorrow to look even nicer: Does polymap in the splancs library not allow horizontal plotting of y-labels? I have tried polymap(studyarea, xlab = "x (m)", ylab = "y (m)", las = 1) but it doesn't change the labels? Mayby some function in library(spatstat) support las? Thanks! Sincerely, Tord

Dear all, For the past three hours I have tried simplify a data frame. I would be really happy if someone could help solving this, I'm sure simple, problem. I want to "aggregate" the data frame: ObjektID BalteNummer Baltessegment S.13 S.13.1 S.13.1.2 S.13 S.13.1 S.13.1.3 S.13 S.13.2 S.13.2.1 S.13 S.13.2 S.13.2.2 S.13 S.13.2 S.13.2.3 S.13 S.13.3 S.13.3.6 S.13 S.13.3 S.13.3.7

Dear all, I have a tree data matrix. For some trees I lack info about tree species, but I want to set them to be spruce. For some reason the tree species names on the remaining (non-NA) rows are changed into numbers (that I do not recognise). I guess that ifelse is not the correct function to use, but I have not found any better one in my searches. Thanks in advance! Sincerely, Tord >

Hi all, I am trying to change data type in the Data Editor but when I do, all numbers in the column changes to NA and back in the Console I get the warning message "NAs introduced by coercion". What's wrong? The table consists of 95 columns and 470 rows. Tord ----------------------------------------------------------------------- ******New e-mail address: Tord.Snall at

Dear all, I'm new to matrix operations in R. I couln't find a solution to the following problem among earlier help mails or in An introd to R, I guess because the question is really basic. I want to extract all above the diagonal, i.e. from 1 2 3 4 1 0 26 49 49 2 26 0 44 40 3 49 44 0 21 4 49 40 21 0 I want 26 49 44 49 40 21 Thanks in advance! Sincerely, Tord

Dear all, I try to summarise my data per category using aggregate, but for some reason I get the error message "sum" not meaningful for factors even though my vector is numeric. The data set is shown below. Could someone please give a hint. Thanks in advance! Sincerely, Tord > names(test) [1] "ObjektID" "tallstubbyta" > is.factor(test$ObjektID);

Dear all, I think that what I want is an as.data.frame.table-object, but see error message below. I have a data frame with one tree per row, diaclass tells if it is a small, mid or large tree > cpy.tradart[1:5, ] ObjektID diaclass 1 AX.Grb.1 bigdia 2 AX.Grb.1 middia 3 AX.Grb.1 middia 4 AX.Grb.1 smalldia 5 AX.Grb.1 middia > I want a data frame telling no of trees per diameter

Dear all, I have installed 1.7.0 on my Win XP and want to use MASS. But I get this error message: > library(MASS) Error in namespaceExport(ns, exports) : undefined exports: abbey, accdeaths, Aids2, Animals, anorexia, austres, bacteria, beav1, beav2, biopsy, birthwt, ... Error in library(MASS) : package/namespace load failed I therefore again downloaded VR.zip. When I extracted the files I

Dear all, I am taking a course in GLM given by a devoted SAS user. He has given us a homework where a mixed GLM with a logistic link and binomially distributed observations should be fitted. I know of the library nlme for mixed effect models but as I understand it, one cannot choose between different links and distributions in the functions provided there. I have so far managed very well with

Dear G?tz Wiegand and Corey Moffet, Thanks for answering while I went to get a cup of coffee. This list is really excellent! The answer, for you other beginners, is in R a tab is represented "\t" Tord > >Dear all, >I want to produce a tab separated text file but R 1.6.2 does not respond to >the tab key on my keyboard. > >When I paste:

Dear all, I have divided two vectors: Np.occup97.98<- as.data.frame(cbind(site = levels(sums$site), Np.occup97.98 = sums$Ant.Nptrad97.98/Ant.trad$Ant.trad97.98)) > Np.occup97.98 site Np.occup97.98 1 erken97 0.342592592592593 2 erken98 0.333333333333333 3 rormyran 0.48471615720524 4 valkror 0.286026200873362 However, at a later stage of the analysis I want >

Dear all, I try to rearrange my ref. database (now in Excel!! :( ) for importing it into a reference manager program (RIS format). My file basically look like this [3,4]-matrix: rbind(c("a", "b", "c", "d"), c("e", "f", "g", "h"), c("i", "j", "k", "l")) [,1] [,2] [,3]

Dear all, I would like to know if there is a function for sorting a data.frame based on one column. sort() does it for a single vector but not for the whole data. order() does it if there are relations between the columns but in my case it is simply a data.frame. Further, I would like to plot other data, pop.size<- c(800, 800, 1500, 1000, 158, 300, 740, 250, 2000, 1500, 250, 700)

Dear all, I have discovered that there are many things that I used to do in my GIS which are easily done directly in R, for example calculating interpoint distances using geoR and pick out points inside a polygon using splancs. I now wonder, is there a function to create a line object like a watercourse and then calculate the distances between many points in space and this line? I couldn't

Dear all, I get an error message when I use polr() in MASS. These are my data: skugg grupp frekv 4 1 gr3 0 5 2 gr3 3 6 3 gr3 6 10 1 gr5 1 11 2 gr5 12 12 3 gr5 1 > > summary(polr(skugg ~ grupp, weights=frekv, data= skugg.cpy1.dat)) Error in optim(start, fmin, gmin, method = "BFGS", hessian = Hess, ...) :