I will swallow my pride and post to this list for the second time in 24
hours--I have a paper due to a reviewer and I am desperate.
Has anybody written code to move the output from summary() called on
the results of a zeroinfl() model (from the pscl package) into a form
suitable for
Hi all,
I'm trying to discover the options available to me for logistic and linear
regression. I'm doing some tests on a dataset and want to see how different
flavours of the algorithms cope.
So far for logistic regression I've tried glm(MASS) and lrm (Design) and
found there is a big
glm is not, and never was. part of the MASS package. It's in the stats package.
Have you sorted out why there is a big difference between the results you get
using glm and lrm?
Are you confident it is due to the algorithms used and not your ability to use
the software?
To be helpful, if
Hi Phil,
So far for logistic regression I've tried glm(MASS) and lrm (Design) and
found there is a big
difference.
Be sure that you mean what you say, that you are saying what you mean, and
that you know what you mean when making such statements, especially on this
list. glm is not in
Hi Chris,
My ideal would be to gather the information onto the clipboard so I
could paste it into Excel and do the formatting there, but any approach
would be better than what I have now.
I would never use Excel for this since there are far superior tools
available. But it is very easy to
Hi,
Thank you for your reply.
I have looked into the epicalc package already and it will not help me
this time.
I realize that I wasn't clear enough. In my example I want to relate
all my data to Country, i.e. the distribution of Hair, Eye color, Sex
(they are dichotomous in my example) and Age
Chris,
If you want to go the clipboard route, use write.table():
d1 - coef(summary(fm_zip2))$count
d2 - coef(summary(fm_zip2))$zero
d - rbind(d1, rep(NA, 4), d2)
write.table(d, file=clipboard, col.names=NA, sep=\t)
Now paste into Excel.
-Peter Ehlers
Mark Difford wrote:
Hi Chris,
My
Hi Bill,
Thanks for you comments. You may be right in that my ability to use the
software may be the problem. I was using lm to fit a model with 'target'
values of 0 or 1. I then discovered there was a lrm model as well, so just
replaced lm with lrm and expected it to be fine. Then I found that
Hello,
Please look here:
http://crantastic.org/search?q=biology
http://crantastic.org/search?q=biologyRegards,
Carlos.
On Fri, Oct 30, 2009 at 6:12 PM, mau...@alice.it wrote:
Is there any R package that implements the same capability of MatLab
toolbox called SimBiology ?
We are expecially
Thank you Ravi
Ravi Varadhan wrote:
Simply do this:
lm(y ~ I(X1 - X2))
Ravi.
Ravi Varadhan, Ph.D.
Assistant Professor,
Division of Geriatric Medicine and Gerontology
School of Medicine
Johns Hopkins
Thank you David
David Winsemius wrote:
On Oct 30, 2009, at 7:47 PM, CE.KA wrote:
Sorry there was a mistake:
I could not see what was different?
Hi R users
I want to do a multiple linear regression with R
With a normal model (Y=Cste+A1*X1+A2*X2) the program would be
Renaud,
Package bigmemory can help you with shared-memory matrices, either in RAM or
filebacked. Mutex support currently exists as part of the package, although
for various reasons will soon be abstracted from the package and provided
via a new package, synchronicity.
bigmemory works
OK, I think I've figured it out, the predict of lrm didn't seem to pass it
through the logistic function. If I do this then the value is similar to
that of lm. Is this by design? Why would it be so?
1 / (1 + Exp(-1 * 3.38)) = 0.967
tdm wrote:
Anyway, do you know why the lrm predict
tdm wrote:
OK, I think I've figured it out, the predict of lrm didn't seem to pass
it through the logistic
function. If I do this then the value is similar to that of lm. Is this
by design? Why would it
be so?
Please take some time to read the help files on these functions so that you
at
On Oct 31, 2009, at 7:29 AM, tdm wrote:
OK, I think I've figured it out, the predict of lrm didn't seem to
pass it
through the logistic function. If I do this then the value is
similar to
that of lm. Is this by design?
Yes, at least for certain meanings of this. When working with
I would have used the search strategy stochastic. Despite the risk
that such a strategy could bring up anything to do with statistics, it
in fact brings me to Stefano Iacus' package sde (and refers to the
accompanying book, Simulation and Inference for Stochastic
Differential Equations
Estaré ausente de la oficina desde el 30/10/2009 y no volveré hasta el
12/11/2009.
Responderé a su mensaje cuando regrese.
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
--- On Fri, 10/30/09, Tony Greig tony.gr...@gmail.com wrote:
From: Tony Greig tony.gr...@gmail.com
Two questions:
1 - Say I have average speed and directions for tide and I
would like to
plot them on a polar plot, but with different colors so I
can indicate the
two directions. I'm
Do you mean to apply the same calculation to each element?
? apply
mydata - data.frame(aa = 1:5, bb= 11:15)
mp5 - function(x) x*5
mp5data - apply(mydata, 2, mp5)
mp5data
This is functionally equivelent to a double if loop.
mydata - data.frame(aa = 1:5, bb= 11:15)
newdata -
Hi,
I am new to R. My problem is with the ordered logistic model. Here is my
question:
Generate an order discrete variable using the variable
wrwage1 = wages in first full calendar quarter after benefit application
in the following way:
*
wage*1*Ordered *=
1 *if*0 *· wrwage*1 * *1000
2
On 10/30/2009 9:20 PM, Stropharia wrote:
Dear R users,
Does anyone know of a package that implements the Johnson-Neyman procedure -
for testing differences among groups when the regression slopes are
heterogeneous (in an ANCOVA model)? I did not get any matches for functions
when searching
This is rather obviously homework, and you have not read the Posting
Guide, and you have not addressed the question of academic integrity
policies that are probably in force at your university.
--
David
On Oct 31, 2009, at 1:30 PM, md. jakir hussain talukdar wrote:
Hi,
I am new to R. My
Hi
Consider
a - 1:5; b - 6:10
x - list(a=a,b=b)
y - list(a,b)
I can get c(a,b) from x using names(x). Is it also possible to get
a and b from y? (The command names(y) gives NULL.)
Thanks in advance.
Chirok
__
R-help@r-project.org mailing list
On Oct 31, 2009, at 12:33 PM, John Kane wrote:
Do you mean to apply the same calculation to each element?
? apply
mydata - data.frame(aa = 1:5, bb= 11:15)
mp5 - function(x) x*5
mp5data - apply(mydata, 2, mp5)
mp5data
It would have been much more compact (and in the spirit of functional
I don't think so. You can only get the name of something that has a
name, and in this case it does not.
-Ista
2009/10/31 Chirok Han chirok@gmail.com:
Hi
Consider
a - 1:5; b - 6:10
x - list(a=a,b=b)
y - list(a,b)
I can get c(a,b) from x using names(x). Is it also possible to get
a
list() does not assign names but data.frame(a, b) does and can be used
at least in the case shown where a and b have the same length. And,
of course, a data frame is a list.
a - 1:5; b - 6:10
names(data.frame(a, b))
[1] a b
2009/10/31 Chirok Han chirok@gmail.com:
Hi
Consider
a - 1:5;
Strange, this normally works, but in a recent run, I have a data set in an
xts format, that has a lot of NA's in a few of the variables in the leading
and trailing positions, due to some lagging calculations. Before running an
analysis, I use
newdata-na.omit(orginaldata)
and normally a
What class is 'orginaldata' (and its columns if a data frame)?
Note that packge xts defines methods for na.omit:
library(xts)
...
methods(na.omit)
[1] na.omit.data.frame* na.omit.default*na.omit.ts*
[4] na.omit.xts*
so this is possibly something peculiar to package xts. See the footer
Hi all,
I am trying to figure out a way to improve my code's efficiency by avoiding the
use of loop.
I want to calculate a conditional mean(?) given time.
For example, from the data below, I want to calculate
sum((value|choice==1)/sum(value)) across time.
Is there a way to do it without using
one approach is the following:
# say 'DF' is your data frame, then
with(DF, {
ind - choice == 1
n - tapply(value[ind], time[ind], sum)
d - tapply(value, time, sum)
n / d
})
I hope it helps.
Best,
Dimitris
parkbomee wrote:
Hi all,
I am trying to figure out a way to improve
what is the dimension of your data?
you might try projecting the points into planes defined by 3 cluster centers.
plot, for each cluster, a density plot or histogram of distances to
the cluster center, and perhaps overlay the density curve for points
not in the cluster.
albyn
Quoting
oops, I see the answer in your question: 15...
Quoting Albyn Jones jo...@reed.edu:
what is the dimension of your data?
you might try projecting the points into planes defined by 3 cluster centers.
plot, for each cluster, a density plot or histogram of distances to
the cluster center, and
On Oct 30, 2009, at 3:29 PM, PDXRugger wrote:
David,
You are correct. I think the frist two assumptions can be thrown
out and
only the latter two (c,d) can be considered. So how would i combine
Acres
for matching Bldgids based on assumptions c,d?
I also think your requested output
On 10/30/2009 10:02 PM, Joe King wrote:
Is there anyway to make the lines in the dispersion command come forward in
a plot and allow the fill in the dispersion parameter be transparent so all
of the lines I am using to note confidence intervals are shown?
Hi Joe,
I have probably missed
On 10/31/2009 04:49 AM, Tony Greig wrote:
Hi,
Two questions:
1 - Say I have average speed and directions for tide and I would like to
plot them on a polar plot, but with different colors so I can indicate the
two directions. I'm using polar.plot from the plotrix library. How can I add
a second
Dear R-Help Team,
as a R novice I have a (maybe for you very simple question), how do I get
the following solved in R:
Let R be a n x n matrix:
\mid R\mid^{-\frac{1}{2}}
solve(A) gives me the inverse of the matrix R, however not the ^(-1/2) of
the matrix...
Thank you very much in advance!
On 10/30/2009 10:35 PM, Huidong TIAN wrote:
Dear friends,
I will be very happy if anyone tell me the way to change work directory
permanently?
I mean not use the function setwd() which can only change temporary, when
you close the console, it will the old directory.
Sys.setenv(R_USER = '') also
On Oct 31, 2009, at 4:39 PM, Kajan Saied wrote:
Dear R-Help Team,
as a R novice I have a (maybe for you very simple question), how do
I get
the following solved in R:
Let R be a n x n matrix:
\mid R\mid^{-\frac{1}{2}}
solve(A) gives me the inverse of the matrix R, however not the
This is pretty much equivalent:
tapply(DF$value[DF$choice==1], DF$time[DF$choice==1], sum) /
tapply(DF$value, DF$time, sum)
And both will probably fail if the number of groups with choice==1 is
different than the number overall.
--
David.
On Oct 31, 2009, at 5:14 PM, Dimitris
Thank you both.
However, using tapply() instead of a loop does not seem to improve my code much.
I am using this inside of an optimization function,
and it still takes more than it needs...
CC: bbom...@hotmail.com; r-help@r-project.org
From: dwinsem...@comcast.net
To:
40 matches
Mail list logo