similar to: Syntax for random effect in glmmPQL

Dear R-listers Is there anyone who knows why I get different eigenvectors when I run MatLab and R? I run both programs in Windows Me. Can I make R to produce the same vectors as MatLab? #R Matrix PA9900<-c(11/24 ,10/53 ,0/1 ,0/1 ,29/43 ,1/24 ,27/53 ,0/1 ,0/1 ,13/43 ,14/24 ,178/53 ,146/244 ,17/23 ,15/43 ,2/24 ,4/53 ,0/1 ,2/23 ,2/43 ,4/24 ,58/53 ,26/244 ,0/1 ,5/43) #R-syntax

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 R-listers I'm a new R-user who needs some help with a test that I want to do. I have done a field experiment: four treatments (cont, x, y and xy) at three sites (A, B and C), the response is count data (0 - 15). I've done a Poisson regression: >glm(response~as.factor(treatment)*as.factor(site), family=quasipoisson, offset(max.response), data=dat) The "offset" is the

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

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 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, I'm working with a population projections matrix model from demographic data. The dominant eigenvalue of the matrix is the growth rate of the population (lambda). I would like to estimate confidence intervals for the lambda values with bootstrap. The function for calculating the eigenvalue I wrote like this: mat.fun <- function(d,i){ mat<-matrix(c(0,

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 a problem with specifying the no. of knots in our function which include gam(). I last worked with this in mid September but since then I have reinstalled R and Simon Wood's library(mgcv), which he has changed since then. The statistician (and good R-coder) with whom I co-operate is now unfortunately overloaded with teaching, and I'm in the sprut of my thesis.... I

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 try to plot a logstic regression model as follows: Valkror<- read.table("Valkror.txt", header=T) np.bark<- glm(Npinc~bark, family=binomial, data=Valkror) plot.formula(np.bark$fitted.values~ Valkror$bark) The above looks nice but I want to connect the points with a line so I try: plot.formula(np.bark$fitted.values~ Valkror$bark, type="b") but it gives a

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