similar to: knots fixed in gam(), library(mgcv)

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

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 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 >

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've used glm(family=binomial(link="logit")) several times, but now I think that a log link is more appropriate. I want to fit a model for probability of tree fall (TF)), with tree diameter (dbh) and soil moisure (soil) as predictors. A large number of trees have been checked every second year whether they stand up (0) or have fallen (1). I assume that the tree fall

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 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 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 get an unexpected error when using the mantel.test function in library(ape). Is there something wrong with my matrices? > dmatr<- read.table("obt.Loiselle.kinship.intens02.DIST.matrix.txt",head=F) > kmatr<- read.table("obt.Loiselle.kinship.intens02.KINSHIP.matrix.txt",head=F) > > mantel.test(as.matrix(kmatr), as.matrix(dmatr), nperm = 500, graph

Dear all, I've worked with optim before but never encountered this error message: Nelder-Mead direct search function minimizer 0.23 0Error: subscript out of bounds The error seems to depend on the initial parameter values. However, strangely (I think), I recieve this error message when I have given very good initial values - the values which were returned when I gave worse values. An

Dear all, I try to calculate the number of occurring levels, per ObjektID, of cpy.busk$TradArt, but this line returns the TOTAL number of levels occurring in cpy.busk$TradArt for each ObjektID. This was difficult to say in English! Did you get it? Is there any function for doing what I want. I have tried searching the archives without success. buskartant<- aggregate(list(buskartant =

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

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 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 have used polr in MASS but I am uncertain about the summary(polr.object) interpretation and would be happy for help on that. This is my summary: > summary(shade.polr) Re-fitting to get Hessian Call: polr(formula = as.ordered(shade) ~ as.factor(objekt), data = sof, weights = as.numeric(frek)) Coefficients: Value Std. Error t value 2.1699520 0.3681840 5.8936612

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