Displaying 20 results from an estimated 2000 matches similar to: "optim() error message"
2003 Dec 07
2
par(las = 1) not possible in polymap(), library(splancs)?
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
2003 Jul 21
2
bold AND italic as font in text()
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
2003 Jan 10
0
Thanks: Re: count levels per factor level
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)])) )
2003 Feb 05
1
simplify a data frame
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
2003 Jan 02
1
replace NA with factor class
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
>
2000 Jun 22
1
NAs introduced when changing data type
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
2003 Jan 02
1
aggregate: "sum" not meaningful for factors
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);
2003 Apr 27
0
Thanks for: no response of the tab key
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:
2003 Sep 09
1
tree mortality risk model using glm()
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
2003 May 05
1
Error in library(MASS) : package/namespace load failed
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
2003 Jan 16
1
problem with as.data.frame.table
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
2003 May 26
0
knots fixed in gam(), library(mgcv)
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
2003 May 26
0
as.matrix() problem in mantel.test()?
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
2003 Jan 09
0
aggregate() nlevels()
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 =
2001 Oct 04
1
sort data.frame and contour()
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)
2001 Nov 21
2
distances from points to line
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
2003 May 22
7
extract half a matrix
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
2003 Feb 25
1
summary(polr.object)
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
2003 Jul 31
2
how as.numeric() !-> factor
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
>
2001 May 27
2
library for mixed GLM?
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