Dear R members
Could someone indicate me how to simulate Likert-type data using the GAMLSS
library and how to fit that data?
Lets say, 50 random numbers in a variable that goes from -4 to +4 in steps of
1 (i.e., -4, -3, -2, -1, 0, 1, 2, 3, 4) with a mean of 2 and SD of 0.15 and
belonging to a no
Takeshi N wrote:
> Hi,there,
>
> I could install "Rgraphviz" in R version 2.6. I could not install this in R
> version 2.7 or 2.8.
> Please try in version2.6.
Rgraphviz requires that additional software be installed, for R 2.9.*,
see the SystemRequirements field on the page at
http://www.bioc
Hey Nathan,
You might like the DBSCAN algorithm.
http://en.wikipedia.org/wiki/DBSCAN
There's an implementation in the 'fpc' package.
http://cran.r-project.org/web/packages/fpc/index.html
-Mose
On Sun, Jun 14, 2009 at 7:36 PM, Dylan
Beaudette wrote:
> On Sun, Jun 14, 2009 at 7:26 PM, Nathan S.
-BEGIN PGP SIGNED MESSAGE-
Hash: SHA1
Dylan Beaudette wrote:
> On Sun, Jun 14, 2009 at 4:39 PM, Nathan S.
> Watson-Haigh wrote:
>> -BEGIN PGP SIGNED MESSAGE-
>> Hash: SHA1
>>
>> Is there a library which is capable of identifying distinct clusters of size
>> n
>> from a series of X
Dear Group,
I have made significant improvements to our Sweave template, have made
the template self-contained (i.e., you can run it yourself and it will
find the datasets it needs), and have included the output pdf file.
This is at http://biostat.mc.vanderbilt.edu/SweaveTemplate . You will
Michael wrote:
I have to use logistic regression...
Good. That's what lrm is dedicated to.
Frank
On Sun, Jun 14, 2009 at 8:04 PM, Frank E Harrell
Jr wrote:
Also it would be useful to compare glm with the lrm function in the Design
package, for speed and memory use.
Frank
David Winsemius
Thanks Gabor, this is quite clever and it's nice to see another way of doing
it (without ifelse).
Mark
On Sun, Jun 14, 2009 at 6:51 PM, Gabor Grothendieck wrote:
> Note that TRUE and FALSE become 1 and 0 when used in arithmetic
> formulae so:
>
> result <- with(DF,
>(type %in% c("A", "B
I have to use logistic regression...
On Sun, Jun 14, 2009 at 8:04 PM, Frank E Harrell
Jr wrote:
> Also it would be useful to compare glm with the lrm function in the Design
> package, for speed and memory use.
>
> Frank
>
>
> David Winsemius wrote:
>>
>> On Jun 14, 2009, at 9:06 PM, Michael wrote:
Thanks very much Marc for the very useful ideas.
I do hope there will be an even more concise way to do this. At present
I'm not using figure environments in Sweave as in stat reports I don't
usually need floats.
I tried re-defining \includegraphics to include \centerline but got into
infin
On Sun, Jun 14, 2009 at 7:26 PM, Nathan S.
Watson-Haigh wrote:
> -BEGIN PGP SIGNED MESSAGE-
> Hash: SHA1
>
> Dylan Beaudette wrote:
>> On Sun, Jun 14, 2009 at 4:39 PM, Nathan S.
>> Watson-Haigh wrote:
>>> -BEGIN PGP SIGNED MESSAGE-
>>> Hash: SHA1
>>>
>>> Is there a library which is
Your error message is because "if" wants a single value and you are giving it a
vector.
Typically you want to use functions "all" or any" to correct this error message
(look them up ?all ?any) and eg if(any(is.na(...))) But in this case to
accomplish the task you're after I don't even t
On Sun, Jun 14, 2009 at 4:39 PM, Nathan S.
Watson-Haigh wrote:
> -BEGIN PGP SIGNED MESSAGE-
> Hash: SHA1
>
> Is there a library which is capable of identifying distinct clusters of size n
> from a series of XY coordinates?
>
> Failing this, I'd like to be able to to something like:
> Using
On Jun 14, 2009, at 9:06 PM, Michael wrote:
Hi all,
I am getting the following error message:
mymodel = glm(response ~ . , family=binomial, data=C);
Error: cannot allocate vector of size 734.2 Mb
In addition: Warning messages:
1: In array(0, c(n, n), list(levs, levs)) :
Reached total alloc
That graphic was created as part of a Ph D thesis in the mid 1990's. I
sent the OP off-list a link to that work and the FBI webpage that
describes access to such data.
--
David
On Jun 14, 2009, at 9:11 PM, spencerg wrote:
Hi, Mario:
The following suggest to me that these data may not
Dear list fellows,
I want to study time series and use R to analyse time series of fishing
data from several species (landings and cpue) investigating the
correlation between them and with environmental factors (water
temperature, wind, etc.).
Searching at Amazon I found three books with exa
Hi, Mario:
The following suggest to me that these data may not be available
in any contributed R package:
library(RSiteSearch)
fbi <- RSiteSearch.function('FBI Homicide')
hits(fbi) # 0
hom <- RSiteSearch.function('Homicide')
hits(hom) # 4
HTML(hom)
The 4 items exposed via the l
Hi all,
I am getting the following error message:
> mymodel = glm(response ~ . , family=binomial, data=C);
Error: cannot allocate vector of size 734.2 Mb
In addition: Warning messages:
1: In array(0, c(n, n), list(levs, levs)) :
Reached total allocation of 1535Mb: see help(memory.size)
2: In ar
Note that TRUE and FALSE become 1 and 0 when used in arithmetic
formulae so:
result <- with(DF,
(type %in% c("A", "B", "C")) *
(1 * (status == "a") + 2 * (status == "b") + 3 * (status ==
"c")) +
(type %in% c("D", "E", "F")) *
(9 * (status == "a") +
On Jun 14, 2009, at 2:28 PM, Frank E Harrell Jr wrote:
Is there an elegant way to get Sweave to center graphics by default?
I'd like to use \centerline{\includegraphics{}} etc. to save some
vertical space that \begin{center} ... \end{center} uses, and I'd
like to avoid centering with each
On Jun 14, 2009, at 6:33 PM, Chuck Cleland wrote:
On 6/14/2009 6:18 PM, Mark Na wrote:
Here's what I need to do (dummy example, my data are more
complicated):
If type = A or B or C
and status = a then count = 1
and status = b then count = 2
and status = c then count = 3
Else
-BEGIN PGP SIGNED MESSAGE-
Hash: SHA1
Is there a library which is capable of identifying distinct clusters of size n
from a series of XY coordinates?
Failing this, I'd like to be able to to something like:
Using a sliding window of size n along the x-axis I'd like to determine the
distanc
Hi Spencer,
You have brought up an important topic. I agree with you that a comparative
evaluation of R packages providing similar capabilities is a worthwhile task.
Your example of "gam" is a very good one. Let me give you a couple of other
examples of comparative evaluation that I am worki
Well, suppose I wanted to suggest changes to some documentation, or write an
alternate help file for some function. Where would I put it?
Let's say I type, in R, ?median. Now suppose I have suggestions. If I look at
http://wiki.r-project.org/rwiki/doku.php?id=&idx=rdoc:base
I don't see the m
On Sun, Jun 14, 2009 at 6:41 PM, Ted
Harding wrote:
> That is an interesting and potentially very fruitful idea. I'd like
> to propose a "quasi-readonly" mode in which wiki users could interpolate
> the comments at chosen points within the original text of a help-page
> (but without any alteration
Hi All,
I have been trying to get this LOO-Cross Validation method to work on R for
the past 3 weeks but have had no luck. I am hoping someone would kindly help
me.
Essentially I want to code the LOO-Cross Validation for the 'Local Constant'
and 'Local Linear' constant estimators. I want to fin
On 14-Jun-09 22:23:24, Gabor Grothendieck wrote:
> In PHP and also in MySQL the manual has a wiki capability
> so that users can add notes at the end of each page, e.g.
>
> http://www.php.net/manual/en/functions.variable-functions.php
>
> http://dev.mysql.com/doc/refman/4.1/en/update.html
>
> Th
On 6/14/2009 6:18 PM, Mark Na wrote:
> Hi,
> I've been struggling most of the morning with an IF ELSE problem, and I
> wonder if someone might be able to sort me out.
>
> Here's what I need to do (dummy example, my data are more complicated):
>
> If type = A or B or C
> and status = a then c
In PHP and also in MySQL the manual has a wiki capability
so that users can add notes at the end of each page, e.g.
http://www.php.net/manual/en/functions.variable-functions.php
http://dev.mysql.com/doc/refman/4.1/en/update.html
That would combine documentation and wiki into one. Here it would
i
Hi,
I've been struggling most of the morning with an IF ELSE problem, and I
wonder if someone might be able to sort me out.
Here's what I need to do (dummy example, my data are more complicated):
If type = A or B or C
and status = a then count = 1
and status = b then count = 2
and
I agree that the documentation is a primary source, made now more
accessible with the availability of the RSiteSearch package that allows
more structures searches of the help pages in contributed packages than
previously available. On the other hand, it can sometimes be difficult
to get t
At 11:02 PM +0200 6/14/09, Martin Batholdy wrote:
hi,
How can I compute a reliability score of a scale consisting only of
dichotomous items?
RSiteSearch("cronbach's alpha") suggests the ltm, psych, psy,
psycho, epicalc, packages (among others).
There are a number of ways of estimating re
I certainly don't have anything against the WIKI, but I think that the
documentation
is where the action is, especially for newbies. It's the natural first step
when you want to learn about a function or when you get an error message you
don't understand.
Peter
Peter L. Flom, PhD
Statistical
John Sorkin wrote:
Perhaps help pages should have links to relevant portions of the WIKI.
John
John David Sorkin M.D., Ph.D.
Chief, Biostatistics and Informatics
University of Maryland School of Medicine Division of Gerontology
Baltimore VA Medical Center
10 North Greene Street
GRECC (BT/
hi,
How can I compute a reliability score of a scale consisting only of
dichotomous items?
thanks for any help!
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-projec
Or you could do what I suggested below. I had assumed you had a
~/bin/ directory and knew how to make a symbolic link. If that
assumption is incorrect, let us know and we can tell you how to do
it.
As Bill suggests, removing the present old version would be
preferable, but I can understand you
Glenn Woodart wrote:
Dear Thomas
Thank you for your excellent help. I agree method 3 would be the best one,
and this is the one I have been trying to get working (like you described).
So far my skills with pointers, what is returned etc in pure C code is
limited, so my attempts have failed so fa
On 15/06/2009, at 7:48 AM, Grześ wrote:
Hello!
I wont to use a function is.na()
I have two vectors:
a=c(1,NA,3,3,3)
b=c(0,0,0,0,0)
and when I use is.na function it's ok:
is.na(a)
[1] FALSE TRUE FALSE FALSE FALSE
but I would create sth like this:
for i in 1:length(a){
if (wsp[i] ==
Grześ wrote:
>
> Hello!
> I wont to use a function is.na()
>
> I have two vectors:
>> a=c(1,NA,3,3,3)
>> b=c(0,0,0,0,0)
> and when I use is.na function it's ok:
>> is.na(a)
> [1] FALSE TRUE FALSE FALSE FALSE
>
> but I would create sth like this:
>
> for i in 1:length(a){
> if (wsp[i]
Hello!
I wont to use a function is.na()
I have two vectors:
> a=c(1,NA,3,3,3)
> b=c(0,0,0,0,0)
and when I use is.na function it's ok:
> is.na(a)
[1] FALSE TRUE FALSE FALSE FALSE
but I would create sth like this:
for i in 1:length(a){
if (wsp[i] == is.na(a)) {b=43}
}
or like this
if(
2009/6/14 Juliet Hannah :
> Hi All,
>
> I am trying to understand panel functions. Let's use this example.
>
>
> library(lattice)
> time<-c(rep(1:10,5))
> y <-time+rnorm(50,5,2)
> group<-c(rep('A',30),rep('B',20))
> subject<-c(rep('a',10),rep('b',10),rep('c',10),rep('d',10),rep('e',10))
> myData <-
Hi All,
I am trying to understand panel functions. Let's use this example.
library(lattice)
time<-c(rep(1:10,5))
y <-time+rnorm(50,5,2)
group<-c(rep('A',30),rep('B',20))
subject<-c(rep('a',10),rep('b',10),rep('c',10),rep('d',10),rep('e',10))
myData <-data.frame(subject,group,time,y)
head(myData)
Maybe once the list is going, people who are willing
to submit their packages to such scrutiny and critique
could join in?
Peter
-Original Message-
>From: Patrick Burns
>Sent: Jun 14, 2009 12:30 PM
>To: stephen sefick
>Cc: r-h...@stat.math.ethz.ch, John Sorkin
>Subject: Re: [R] a prop
Is there an elegant way to get Sweave to center graphics by default?
I'd like to use \centerline{\includegraphics{}} etc. to save some
vertical space that \begin{center} ... \end{center} uses, and I'd like
to avoid centering with each <>=
Thanks
Frank
--
Frank E Harrell Jr Professor and Chai
Perhaps help pages should have links to relevant portions of the WIKI.
John
John David Sorkin M.D., Ph.D.
Chief, Biostatistics and Informatics
University of Maryland School of Medicine Division of Gerontology
Baltimore VA Medical Center
10 North Greene Street
GRECC (BT/18/GR)
Baltimore, MD 2
The documentation for the "nlme" package was improved a few years
ago by an informal process of this nature. When I first started
following r-help, I answered many questions suggesting in part the the
person read some portion of Pinheiro and Bates (2000). At that time,
none of the help p
That make sense to just do the R-core documentation. I was suggesting
using my package because I understand it and could have a dialog about
the design choices for the documentation, which would also improve the
package itself. I am starting my PhD program tomorrow and will
probably not have enou
I had only envisioned the scope of the
list to be the packages that R-core are
responsible for.
Expanding to contributed packages would
expand both potential usefulness and
certain complexity exponentially.
Pat
stephen sefick wrote:
I have a package StreamMetabolism. I believe that documenta
I neglected to say that I don't have the
time to be one of the "drivers" for the
proposed list.
Pat
Patrick Burns wrote:
Proposal
That a new mailing list be established
that pertains exclusively to R documentation.
The purpose of the list would be to discuss
weak sections of the documentation
Have you actually looked at the help page for layout()?
--
David Winsemius
On Jun 14, 2009, at 7:35 AM, luscinia wrote:
Hey everybody,
I am trying to generate multiple boxplots on one page by using the
layout ()
function. This works really well but the plots are always badly
proportion
Hi,there,
I could install "Rgraphviz" in R version 2.6. I could not install this in R
version 2.7 or 2.8.
Please try in version2.6.
ram basnet wrote:
>
> Dear R users,
>
> I am not so used to this R software. I have to use the package "
> Rgraphviz" but found some problem in the installatio
It woks fine.
Thanks very much.
Olivier
--
Olivier G. Nuñez
Email: onu...@iberstat.es
Tel : +34 663 03 69 09
Web: http://www.iberstat.es
El 14/06/2009, a las 4:30, Michael Lawrence escribió:
I would recommend just
Hey everybody,
I am trying to generate multiple boxplots on one page by using the layout ()
function. This works really well but the plots are always badly proportioned
(3 times wider than high) and much too small to be easily readable. I have
tried reducing margins to allow for bigger plots, but
I have a package StreamMetabolism. I believe that documentation is
the toughest part. I find it to be straight foward, but then I wrote
the package. Lets try one? I don't m
ind helping.
Stephen Sefick
On Sun, Jun 14, 2009 at 7:23 AM, John Sorkin wrote:
> Issues
> (1) Will people use the new
Issues
(1) Will people use the new listserver?
(2) Will developers respond to postings
Question
(1) Are there any guidelines for creating documentation? If not should they be
developed? It seems to me that every help page should include (a) examples, and
when approprate (b) a reference to an arti
Patrick Burns wrote
>Proposal
>
>That a new mailing list be established
>that pertains exclusively to R documentation.
>The purpose of the list would be to discuss
>weak sections of the documentation and
>establish fixes for those weak spots.
>
>
>Pro
>
>If it works, there would be better documen
Proposal
That a new mailing list be established
that pertains exclusively to R documentation.
The purpose of the list would be to discuss
weak sections of the documentation and
establish fixes for those weak spots.
Pro
If it works, there would be better documentation.
It would be an excellent
Hi,
As Uwe said, this is not really the place to get support for Ubuntu.
They have excellent forums. But the problem seems to be that GRUB was
not correctly installed. Run an ubuntu live cd and reinstall GRUB from
that. The command might be grub-install, but i'm not sure
cheers,
Paul
Le
57 matches
Mail list logo
