Displaying 20 results from an estimated 400 matches similar to: "factor as seq() in for loop"
2010 Feb 05
2
glm models with more than one response
Hi everyone,
I am trying to construct a glm and am running into a couple of questions.
The data set I am using consists of 6 categories for the response and 6
independent predictors representing nutrient concentrations at sample point
locations. Ultimately I'd like to use the probabilities for each response
category in a simulation model such that these probabilities are used to
define a
2003 Jun 25
3
Help
hi
my name is Tomer,im using R version 1.4
i've encontered a couple of problems with my my R programing.
1. i cant understand what is the whights in the nnet$wts command(which
weight fit to the nodes).
2.my OS is Windows XP , i want to operate R from the command line and i want
it to operate a file automaticly(the commands written in a text file) ,is it
possible?( i know that in unix it
2003 Nov 04
1
glm offset and interaction bugs (PR#4941)
Full_Name: Charles J. Geyer
Version: 1.8.0
OS: i686-pc-linux-gnu (Suse 8.2)
Submission from: (NULL) (134.84.86.22)
Two bugs (perhaps related, perhaps independent) revealed by the same
Poisson regression with offset
mydata <- read.table(url("http://www.stat.umn.edu/geyer/5931/mle/seeds.txt"))
out.fubar <- glm(seedlings ~ burn01 + vegtype * burn02 +
offset(log(totalseeds)),
2009 Oct 27
1
Using a variable in the formula
I am sure this question has come up, but searching hasn't given me any
results.
So I need to enter this line:
mx1 <- randomForest(X1 ~ elevation + slope + vegtype, data = moths.train)
But the problem is that X1 is currently hard coded. I would instead like to
be able to put in the value of X1 through a list. For example:
list <- list("X1", "X2", "X3")
2006 Oct 15
1
gamma distribution don't allow negative value in GLMs?
Dear friends,
when i use glm() to fit my data, i use
glm(formula = snail ~ vegtype + mhveg + humidity + elevation + soiltem, *family
= Gamma(link = inverse),* data =a,))
It shows: error in eval(expr, envir, enclos) : *gamma distribution don't
allow negative value*.
But i use
result<-glm(formula = snail ~ vegtype + mhveg + humidity + elevation +
soiltem, family = poisson, data =a) #this
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
2007 Mar 19
5
order of values in vector
Dear all,
I would like to get the order of the values in a vector. I have tried
rank(), order() and searched the archive, though without success.
Here is an example of a try
x= c(20,30,50,40,60,10)
cbind(sort.list(x),x)
x
[1,] 6 20
[2,] 1 30
[3,] 2 50
[4,] 4 40
[5,] 3 60
[6,] 5 10
but I was hoping to get this:
x
[1,] 2 20
[2,] 3 30
[3,] 5 50
[4,] 4 40
[5,] 6 60
[6,] 1 10
I'm
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 Jan 22
3
Error when using polr() in MASS
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, ...) :
2003 Nov 23
4
remove 0 rows from a data frame
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
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
2002 Jan 04
4
line up a matrix
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]
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
2006 Dec 27
3
counties in different colours using map()
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
2003 Oct 08
2
binomial glm warnings revisited
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
2003 Oct 20
4
warning from return() in 1.8 but not in 1.7.0 (PR#4687)
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
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
>
2009 Oct 23
3
opposite estimates from zeroinfl() and hurdle()
Dear all,
A question related to the following has been asked on R-help before, but
I could not find any answer to it. Input will be much appreciated.
I got an unexpected sign of the "slope" parameter associated with a
covariate (diam) using zeroinfl(). It led me to compare the estimates
given by zeroinfl() and hurdle():
The (significant) negative estimate here is surprising, given
2010 Aug 19
1
GLMM random effects
Hello,
I have a couple questions regarding generalized linear mixed models
specifically around fitting the random effects terms correctly to account
for any pseudo-replication.
I am reading through and trying to follow examples from Zuur et al. Mixed
Effects Models and Extensions in Ecology with R, but am still at bit unsure
if I am specifying the models correctly.
Background information:
Our
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
>