similar to: spatial library

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 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, 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, ...) :

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

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

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

Hello all, I am currently working with rpart to classify vegetation types by spectral characteristics, and am comming up with poor classifications based on the fact that I have some vegetation types that have only 15 observations, while others have over 100. I have attempted to supply prior weights to the dataset, though this does not improve the classification greatly. Could anyone supply some

I have a vague recollection of seeing a R-function that will output a table or matix in a format suitable for the tabular enviroment in LaTeX. But I cannot find it. Any hints? Bendix ---------------------- Bendix Carstensen Senior Statistician Steno Diabetes Centre Niels Steensens Vej 2 DK-2820 Gentofte Denmark tel: +45 44 43 87 38 mob: +45 28 25 87 38 fax: +45 44 43 73 13 bxc at novo.dk

Hi, i have some problems to get means, but the data have got NA's which should not be replaced or delete ! tapply(data,var,mean) # makes me problem,because when one value is NA mean is NA,too. Another attempt is this, because the summary does what i want , but is not subsetable to summary[4] if i use tapply,table and the the function below gives me again all means as NA ? xmean

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

Hi, I would like to plot a map of US counties using different colors. map() seems to be the function to use, e.g. library(maps); map('usa'); map('county', 'colorado', add=T,fill = T, col=c(1:5)) plots Colorado counties using colours 1 to 5. However, I want each color to represent a certain value - a value to be picked from a data frame. This code should show a

When I read a table (the dimnames of the table are numbers), R reads the column names as: X plus the number (X1,X2,X3,...) and the rownames as the numbers (1,2,3,...) How can I suppress this action? The code and the output: ^^^^^^^^^^^^^^^^^^^^^^^^ > r2<-read.table("h:\\matrix",header=T,dec=",",na.strings="-99") > dimnames(r2) [[1]] [1] "500"

Rather than spend a tedious interactive session picking several (in fact most) of the library packages off the CRAN site one by one, it would be handy to be able to grab the lot in one go (they seem to add up to about 24MB, so it's not a huge download). Is there a way to do this? Thanks for the help. Ted. -------------------------------------------------------------------- E-Mail: (Ted

I'm using R to manipulate some data for testing nebulisers. I have a data frame as follows: Volt ABal 0.4208 A1 0.3982 AF 0.4353 B1 0.3954 BF 0.4355 C1 0.3975 CF 0.4345 I also have a calibration curve loaded Voltage NaF LogNaF 1 0.3588 20 1.30103 2 0.3990 100 2.00000 3 0.4391 500 2.69897 If I try to