Displaying 20 results from an estimated 700 matches similar to: "Zero inflated count models"
2003 Jun 03
3
Update VR_7.1-6
The update of VR by downloading VR_7.1-6.zip and using install.packages
(from local zip files) fails with the following error message:
Error in file(file, "r") : unable to open connection
In addition: Warning message:
cannot open file `VR/DESCRIPTION'
Other packages can be installed without problems, except of dse_2003.4-1
with a similar error message.
Why?
Operating System:
2003 May 08
3
Avoiding loops to spare time and memory
Is it possible to avoid the loop in the following function (or make the
function otherwise more efficient) and can someone point me to a
possible solution? (It would be great if hours could be reduced to
seconds :-).
# ---------------------------------------------
RanEigen=function(items=x,cases=y,sample=z)
{
X=matrix(rnorm(cases*items),nrow=cases,byrow=F)
S=crossprod(X-rep(1,cases) %*%
2005 Mar 14
2
par(new=TRUE) vs par(new=FALSE)
Just out of curiosity (and I hope that my question does not confuse
those who intuitively understood the par(new=TRUE) statement right in
the beginning):
I would like to know whether beginners (learning R) had the same
difficulty as I had: Intuitively I thought that par(new=TRUE) would draw
a new plot (and not into an already existing plot) and that
par(new=FALSE) would not draw a new plot
2005 May 23
3
skewness and kurtosis in e1071 correct?
I wonder whether the functions for skewness and kurtosis in the e1071
package are based on correct formulas.
The functions in the package e1071 are:
# --------------------------------------------
skewness <- function (x, na.rm = FALSE)
{
if (na.rm)
x <- x[!is.na(x)]
sum((x - mean(x))^3)/(length(x) * sd(x)^3)
}
# --------------------------------------------
and
#
2005 Sep 20
1
Problem with read.spss() and as.data.frame(), or: alternative to subset()?
Trying to select a subset of cases (rows of data) I encountered several
problems:
Firstly, because I did not read the help to read.spss() thoroughly
enough, I treated the data read as a data frame. For example,
dr2000 <- read.spss('myfile.sav')
d <- subset(dr2000,RBINZ99 > 0)
and thus received an error message (Object "RBINZ99" not found), because
dr2000 is not a
2006 Aug 13
3
How to reply to a thread if receiving R-help mails in digest form
I receive R-help messages in digest form that makes it difficult to
answer a post. I noticed that my answer is not added to the thread
(instead, a new thread is started) although I use the same subject line
(starting with "Re: ") as the original post. Is there a solution (I
prefer the digest to separate messages for several reasons and don't
want to change my email reader)?
The
2006 Jul 24
1
conflict of package "Matrix" and summary of lme object
After loading the package "Matrix" (version 0.995-12), using the summary
function with an lme (package nlme version 3.1-75) object results in an
error message saying
Fehler in dim(x) : kein Slot des Namens "Dim" f?r dieses Objekt der
Klasse "correlation"
(translated: 'Error in dim(x) : no slot of the name "Dim" for this
object of class
2006 Jul 24
3
standardized random effects with ranef.lme()
Using ranef() (package nlme, version 3.1-75) with an 'lme' object I can
obtain random effects for intercept and slope of a certain level (say:
1) - this corresponds to (say level 1) "residuals" in MLWin. Maybe I'm
mistaken here, but the results are identical.
However, if I try to get the standardized random effects adding the
paramter "standard=T" to the
2005 Jul 19
2
using argument names (of indeterminate number) within a function
Although I tried to find an answer in the manuals and archives, I cannot
solve this (please excuse that my English and/or R programming skills
are not good enough to state my problem more clearly):
I want to write a function with an indeterminate (not pre-defined)
number of arguments and think that I should use the "..." construct and
the match.call() function. The goal is to write
2005 May 24
0
skewness and kurtosis in e1071 correct? (correction)
I'm sorry, but my previous message, as often happens, some brackets were
wrong:
Here are the correct formulas:
sd = 1/n * sum((x-mean(x))^2) # (1)
sd = 1/(n-1) * sum((x-mean(x))^2) # (2)
This also occured in the last paragraph.
Dirk
*************************************************
Dr. Dirk Enzmann
Institute of Criminal Sciences
Dept. of Criminology
Edmund-Siemers-Allee 1
D-20146
2012 Apr 05
4
Appropriate method for sharing data across functions
In trying to streamline various optimization functions, I would like to have a scratch pad
of working data that is shared across a number of functions. These can be called from
different levels within some wrapper functions for maximum likelihood and other such
computations. I'm sure there are other applications that could benefit from this.
Below are two approaches. One uses the <<-
2008 Mar 03
0
Fwd: R: Studdy Missing Data, differentiate between a percent with in the valid answers and with in the different missing answers
Begin forwarded message:
> Subject: R: Studdy Missing Data, differentiate between a percent
> with in the valid answers and with in the different missing answers
> take a look at the Website from Dirk Enzmann:
> http://www2.jura.uni-hamburg.de/instkrim/kriminologie/Mitarbeiter/
> Enzmann/Software/Enzmann_Software.html
> He provides a function called "freq.r" - maybe
2009 Apr 20
1
doing zero inflated glmm for count data with fmr
Hello R users,
Doing My PhD I collected count data which I believe is zero inflated. I have
run a statistical model with lmer and family=poisson and got
summary(model)@sigma=1 so I believe there is no overdispertion. I would
like to use the fmr function from the 'gnlm' library but I just cannot
figure out from the examples in the help page and some forums out there how
to convert the lmer
2001 May 09
3
odesolve check fails
Hi,
I just installed the odesolve package and ran the check command on it. It
failed trying to execute
library(gnlm).
Sure enough, there is no gnlm library on my system and I could not find it
on the CRAN archive either. Am I missing anything obvious or is gnlm some
private library that somehow found its way into the example section of
odesolve?
Thanks in advance,
Andy
2012 May 07
1
Repeating
Dear All,
I have a codes which calculates the result of Ripley's K function of my
data. I want to repeat this process 999 times. However, i am getting an
error when i use the "for i in" function. Is there any way to repeat this
analysis 999 times. Here are the codes i used ;
data4 <- matrix(c(sample(id),data1),203,3)
a <- data4[,1]
random.case=data4[a==0,]
2004 Dec 13
3
Percentages in contingency tables *warning trivial question*
I hesitate to post this question in the light of recent threads, indeed
I have hesitated for several weeks, however I have come to a full stop
and really need some help if I am going to progress. I am a new user of
R for medical statistics. I have attempted to read all the relevant
documents, but would welcome any suggestions as to what I have missed.
I am trying to contruct "table 1"
2003 Oct 29
2
Where is rmutil package?
Pursing my earlier question, when I tried loading Lindsey's gnlm, I got
a
message
Loading required package: rmutil
Warning message:
There is no package called 'rmutil' in: library(package, character.only
= TRUE, logical = TRUE, warn.conflicts = warn.conflicts,
According to the R documentation
http://finzi.psych.upenn.edu/R/doc/html/packages.html
rmutil is in the standard
2003 Oct 29
1
One inflated Poisson or Negative Binomal regression
Hello
I am interested in Poisson or (ideally) Negative Binomial regression
with an inflated number of 1 responses
I have seen JK Lindsey's fmr function in the gnlm library, which fits
zero inflated Poisson (ZIP) or zero inflated negative binomial
regression, but the help file states that for ' Poisson or related
distributions the mixture involves the zero category'.
I had thought
2000 Oct 27
1
- Estimate LD50 with bnlr{gnlm}
Hi,
I'm not yet familar with GLM and still learning.
How can I perform a BNL (to estimate LDp values) with matrixes like this
(N indicates the observed objects):
data from:
Kerr D, Meador J. 1996. "Modeling Dose Response Using Generalized Linear
Models". Environ Toxicol Chem 15(3)
conc respond n
0 0 20
0 1 19
0 1 20
0 0 20
13.5 5 23
13.5 2 20
29 9 20
29 6 20
53 12 20
53 15 20
85
2002 Sep 27
3
Retaining regularly used add-ons
Hi everyone,
How might I go about configuring R to keep
add-on packages loaded from session to session?
Is this undesireable for some reason?
At present, I keep a file called "pckgs.txt" in my
working directory with, e.g.,
library(gnlm)
library(Hmisc)
...
and then type source("add-ons.txt") every time I start
a new session. I suspect there's a more elegant way
to hold