Hi Murali,
So, the solution to your problem will be to explicitly convert your
matrix to a numeric matrix. Maybe matrix (?matrix) will do, or you'll
also have to use as.numeric (?as.numeric).
The strings on the left seem to me to be the row labels, right..?, not
elements in the matrix/table.
As a
On 17/07/2008 9:47 PM, Murali K wrote:
Hello,
I am a newcomer to R and therefore apologize for posting such a
basic question. I am
trying to multiply 2 matrices t(X1)%*%X1, where t(X1) is:
It's hard to say for sure, but it looks as though X1 really isn't a
numeric matrix. The "ones" a
Hello,
I am a newcomer to R and therefore apologize for posting such a
basic question. I am
trying to multiply 2 matrices t(X1)%*%X1, where t(X1) is:
1 2 3 4 5 8 12 13 20 24 26 27 31 33 34 36 37 40 41 42 45 46
47 48 49
ones 1 1 1 1 1 1 1 1 1 11 1 1 1
on 07/17/2008 08:44 PM [EMAIL PROTECTED] wrote:
Hello R users:
I was trying to draw a boxplot?using 2 variables (rain and wind).
Both of them has length 25056. So dummy is a matrix with dimension 2
x 25056. First I tried to draw a full boxplot using the following
code?and gave an error
boxplot
On Fri, Jul 18, 2008 at 9:44 AM, <[EMAIL PROTECTED]> wrote:
> Hello R users:
>
> I was trying to draw a boxplot?using 2 variables (rain and wind). Both of
> them has length 25056. So dummy is a matrix with dimension 2 x 25056. First I
> tried to draw a full boxplot using the following code?and g
What you should be doing is to return a value from the function and
you can then assign this to an object of your choice:
sizeattributes <- c("length", "width", "height")
blue <- myfunction(bluefrenchwidgets, sizeattributes, "norm_")
red <- myfunction(redcanadianwidgets, sizeattributes, "norm_")
Sorry, I forgot to mention that I also tried??? memory.size(T)? and then tried
to use those commands. Nothing got changed.. Can you please tell me where is
the problem in?my code and how to fix it ?
-Original Message-
From: [EMAIL PROTECTED]
To: r-help@r-project.org
Sent: Thu, 17 Jul 20
Hello R users:
I was trying to draw a boxplot?using 2 variables (rain and wind). Both of them
has length 25056. So dummy is a matrix with dimension 2 x 25056. First I tried
to draw a full boxplot using the following code?and gave an error
> boxplot(dummy$Rain~dummy$Wind)
Error: cannot allocat
On Fri, Jul 18, 2008 at 5:09 AM, Rolf Turner <[EMAIL PROTECTED]> wrote:
>
> On 18/07/2008, at 8:42 AM, Steve Murray wrote:
>
>
>
>> So what I'm looking to do is, 'pad out' the shorter file, by adding in the
>> rows with those that are 'missing' from the longer file (ie. if a particular
>> c
I believe that a short answer to your question is that the
"smooth" is a linear combination of B-spline basis functions, and the
coefficients are the weights assigned to the different B-splines in that
basis.
Before offering a much longer answer, I would want to know what
problem yo
Hi,
I have issue a jacktest command and it produces me a list of numbers. I
would like to map them to a dataframe so I could save it to a file.
I tried to convert the mvr.obj to vector but when I check the class, it is
still jacktest.
converted<-as.vector(mvr.obj)
Or any otehr way to map this t
Thanks to everyone who replied. I really appreciate your help, and have
been trying variations on the solutions you've offered. Unfortunately, none
do quite what I'm hoping to do. Uwe's suggestion seems to be the closest (I
think), but the new variables created within the function didn't exist
o
>Davide Massidda wrote:
> Dear all,
>I have installed R on Linux/Ubuntu 8.04.
> Erik Iverson wrote:
> But you don't say how. Are you compiling R yourself or installing the
Ubuntu package from CRAN?
I comple R using R-2.7.1.tar.gz, I don't use the package from CRAN.
[[alternative HTML ve
Davide Massidda wrote:
Dear all,
I have installed R on Linux/Ubuntu 8.04.
But you don't say how. Are you compiling R yourself or installing the
Ubuntu package from CRAN?
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/
On Fri, Jul 18, 2008 at 01:09:25AM +0200, Davide Massidda wrote:
> Dear all,
> I have installed R on Linux/Ubuntu 8.04. When I try to load the tcltk
> package, I get the response:
>
> > library(tcltk)
> Error in firstlib(which.lib.loc, package) :
> Tcl/Tk support is not available on this system
On Thu, 17 Jul 2008, [EMAIL PROTECTED] wrote:
Hi everybody,
This has been causing me some trouble:
I want to combine several lists of pairs into a single?vector but keep the
pairs together. The lists are of the form:
(1)?x1 x2 x3... N?
(2) y1 y2 y3..? N
I would like to keep it this way, simp
Dear all,
I have installed R on Linux/Ubuntu 8.04. When I try to load the tcltk
package, I get the response:
> library(tcltk)
Error in firstlib(which.lib.loc, package) :
Tcl/Tk support is not available on this system
Error in library(tcltk) : . First.lib failed for 'tcltk'
In order to solve thi
Mi Mariana:
on 07/17/2008 12:40 PM Mariana Carla Gambarotta wrote:
Hi, everybody. Sorry to bring up the subject again but I have visited
the revolution computing web page, but its not clear when or what
will be the new release be. Does anybody have information about
that?. Does anybody know if
Hi everybody,
This has been causing me some trouble:
I want to combine several lists of pairs into a single?vector but keep the
pairs together. The lists are of the form:
(1)?x1 x2 x3... N?
(2) y1 y2 y3..? N
I would like to keep it this way, simply appending one list to another, e.g.
(1) x1i x
Hi Ptit,
>> I would like to fit data with the following formula :
>> y=V*(1+alpha*(x-25))
>> where y and x are my data, V is a constant and alpha is the slope I'm
>> looking for.
Priorities first: lm() or ordinary least-square regression is a basically a
method for finding the best-fitting strai
On 18/07/2008, at 8:42 AM, Steve Murray wrote:
So what I'm looking to do is, 'pad out' the shorter file, by adding
in the rows with those that are 'missing' from the longer file (ie.
if a particular coordinate isn't present in the shorter file but is
in the 'longer/master' file),
Thanks for all the excellent replies.
Iain
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and p
Below is an example of a problem I encounter repeatedly when I write
functions. A call works at the command line, but it does not work inside a
function, even when I have made sure that all required variables are
available within the function. The only way I know to solve it is to make
the require
I like what smooth.spline does but I am unclear on the output. I can see from
the documentation that there are fit.coef but I am unclear what those
coeficients are applied to.With spline I understand the "noraml" coefficients
applied to a cubic polynomial. But these coefficients I am not sure ho
Steve
?merge
and the all collection of arguments
Peter Alspach
> -Original Message-
> From: [EMAIL PROTECTED]
> [mailto:[EMAIL PROTECTED] On Behalf Of Steve Murray
> Sent: Friday, 18 July 2008 8:43 a.m.
> To: r-help@r-project.org
> Subject: [R] Matching Up Values
>
>
> Dear all,
>
Below is an example of a problem I encounter repeatedly when I write functions.
A call works at the command line, but it does not work inside a function, even
when I have made sure that all required variables are available within the
function. The only way I know to solve it is to make the re
Dear all,
I have two files, both of similar formats. In column 1 are Latitude values
(real numbers, e.g. -179.25), column 2 has Longitude values (also real numbers)
and in one of the files, column 3 has Population Density values (integers);
there is no column 3 in the other file.
However, the
On Thu, 17 Jul 2008, Mike Lawrence wrote:
Hi all,
I'm looking for an analytic method to obtain the PDF & CDF of the sampling
distribution of a given correlation (rho) at a given sample size (N).
See Fisher (1915) Biometrika.
It is non-central t up to a monotone transformation.
Finding th
This is really odd. Because in Suse I could have compiz on and still
have rotating rgl plots. I did have kde in suse and now gnome in ubuntu
but I wouldn't have thought that would make the difference. It must be a
ubuntu thing. Very odd.
Paul
Duncan Murdoch wrote:
On 7/17/2008 2:54 PM, H. Pa
On 7/17/2008 2:54 PM, H. Paul Benton wrote:
Hello all,
I've switched over my OS to ubuntu from Suse 10.3. In suse I was able to
load up the 'rgl' library in R and then move my plots around after I had
plotted them. Which I assumed was part of the openGL. However, since
switching to ubuntu I a
Hello all,
I've switched over my OS to ubuntu from Suse 10.3. In suse I was able to
load up the 'rgl' library in R and then move my plots around after I had
plotted them. Which I assumed was part of the openGL. However, since
switching to ubuntu I am unable to 'move' my plot around, they are j
?print
?cat
HTH
Stephan
[EMAIL PROTECTED] schrieb:
Hi,
I know this must be a stupid question, and sorry in advance for being such a
noob. But, is there way to get R to display only certain variables when running
a script. I know if you want to see the value of a variable when using the
interfa
on 07/17/2008 12:40 PM Mariana Carla Gambarotta wrote:
Hi, everybody. Sorry to bring up the subject again but I have visited
the revolution computing web page, but its not clear when or what
will be the new release be. Does anybody have information about
that?. Does anybody know if the version t
Depending on your motivation here, you may want to search for 'debugging
in R' or something to that effect to look at the various options
available. I like the debug package available on CRAN. Also see
?browser and ?trace if you want to debug a function.
[EMAIL PROTECTED] wrote:
Hi,
I kn
?print
?cat
On Thu, Jul 17, 2008 at 1:47 PM, <[EMAIL PROTECTED]> wrote:
> Hi,
>
> I know this must be a stupid question, and sorry in advance for being such a
> noob. But, is there way to get R to display only certain variables when
> running
> a script. I know if you want to see the value of a
Monica Pisica wrote:
Hi,
I would like to compare differences in AUC from 2 different models, glm and gam
for predicting presence / absence. I know that in theory the model with a
higher AUC is better, but what I am interested in is if statistically the
increase in AUC from the glm model to th
Hi,
I know this must be a stupid question, and sorry in advance for being such a
noob. But, is there way to get R to display only certain variables when running
a script. I know if you want to see the value of a variable when using the
interface, you just type it in and hit enter, but is there a c
Hi, everybody. Sorry to bring up the subject again but I have visited the
revolution computing web page, but its not clear when or what will be the new
release be. Does anybody have information about that?. Does anybody know if the
version that will be released will be validated?. I have been
On 17-Jul-08 02:32:41, Frank E Harrell Jr wrote:
> Barry Rowlingson wrote:
>> 2008/7/16 Ted Harding <[EMAIL PROTECTED]>:
>>
>>> Hi Folks,
>>> I'd be grateful for good suggestions about the following.
>>>
>>> I'm plotting a family of (X,Y) curves (for different levels
>>> of another variable, Z): s
You seems to be confused here. A function is "invoked" -- i.e. "called" --
via the expression
fun(...) ## *Note the parentheses.* The function body is *not* printed.
The expression
fun ## note: no parentheses
is not a call to the function. Rather, it is by default a call to the
**print** m
Tudor Bodea wrote:
Dear UseRs:
I intend to write a package to handle the basic operations a revenue management
analyst has to deal with on a regular basis (e.g., demand untruncation,
capacity allocation, pricing decisions, etc.). Following the directions posted
by Peter Rossi (Making R Packages
Does this do at least the means for you:
> x <- read.csv(textConnection("Col1 , Col2
+ A, 3
+ , 2
+ , 3
+ B, 4
+ , 5
+ , 4
+ C , 1
+ , 4
+ , 3"), strip.white=TRUE)
> x
Col1 Col2
1A3
2 2
3 3
4
on 07/17/2008 10:57 AM Tudor Bodea wrote:
Dear UseRs:
I intend to write a package to handle the basic operations a revenue management
analyst has to deal with on a regular basis (e.g., demand untruncation,
capacity allocation, pricing decisions, etc.). Following the directions posted
by Peter R
> "MM" == Martin Maechler <[EMAIL PROTECTED]>
> on Wed, 16 Jul 2008 18:02:47 +0200 writes:
> "MC" == Mark Cowley <[EMAIL PROTECTED]>
> on Wed, 16 Jul 2008 15:32:30 +1000 writes:
MC> Dear list,
MC> I am analysing a set of quantitative proteomics data
MC> from 1
Hi, hope this will help:
txt <- "Col1,Col2
A,3
, 2
,3
B,4
, 5
, 4
C,1
, 4
, 3"
## read data
dat <- read.csv(textConnection(txt),na.string="")
## fill in empty cells with correct category
dat$Col1[] <-
Reduce(function
Nice tip on filling missing col1 values. I've always just done the
filling in excel by hand. I then use that same hand to smack the
person that gave me data in excel format :Op
On 17-Jul-08, at 1:19 PM, Erik Iverson wrote:
Iain Gallagher wrote:
Hello list
I have been given some Excel sheet
Try this:
#x <- read.csv('your_file.csv')
x$Col1 <- rep(as.character(x$Col1[x$Col1!="" ]),
each = unique(diff(which(x$Col1 != ""
with(x, sapply(c(sd, mean), function(x)tapply(Col2, Col1, x)))
On Thu, Jul 17, 2008 at 12:50 PM, Iain Gallagher
<[EMAIL PROTECTED]> wrote:
> Hello li
Iain Gallagher wrote:
Hello list
I have been given some Excel sheets with data laid like this:
Col1Col2 A 3 2 3 B 4 5 4 C 1 4 3
I was hoping to import this into R as a csv and then get the mean and
SD for each letter in column 1.
Could someone give me some guidance
Hi all,
I know pnorm() is able to calculate a sequence of probabilities by providing
a vector of means, and a single number for location and standard deviation.
For example, pnorm(0.5,c(0.5,2),1) will produce [1] 0.500 0.0668072.
However, I wonder if its multivariate counterpart mvtnorm can do
See Chapter 22 in N.L. Johnson, S. Kotz, and N. Balakrishnan,
"Continuous Univariate Distributions", Volume 2, Second Edition, 1995.
Maria
On Thu, Jul 17, 2008 at 11:29 AM, Mike Lawrence <[EMAIL PROTECTED]> wrote:
> Hi all,
>
> I'm looking for an analytic method to obtain the PDF & CDF of the s
Dear UseRs:
I intend to write a package to handle the basic operations a revenue management
analyst has to deal with on a regular basis (e.g., demand untruncation,
capacity allocation, pricing decisions, etc.). Following the directions posted
by Peter Rossi (Making R Packages under Windows: A Tut
Hello,
I would like to fit data with the following formula :
y=V*(1+alpha*(x-25))
where y and x are my data, V is a constant and alpha is the slope I'm
looking for.
How to translate this into R-language ?
At the moment, I only know : lm(y ~ x)
Sorry for such a basic question. I thought I could
I have microarray data with gene names in the first column, gene id in the
second and the expression data in the remaining columns. When trying use
read.table I get the error, . . . "more columns than column names". Is
there any way to keep both columns of "names" without having to discard one.
Hello list
I have been given some Excel sheets with data laid like this:
Col1Col2
A 3
2
3
B 4
5
4
C 1
4
3
I was hoping to import this into R as a csv and then get the mean and SD for
each letter in
"Leandro Marino" <[EMAIL PROTECTED]> wrote in message
news:[EMAIL PROTECTED]
But this pearson want to specify the colours. This person want me to create
an pallete of colours like that:
Name of color - Code - Colour
White - 0xFF - color white (like an box with this color)
This page shows
Hi all,
I'm looking for an analytic method to obtain the PDF & CDF of the
sampling distribution of a given correlation (rho) at a given sample
size (N).
I've attached code describing a monte carlo method of achieving this,
and while it is relatively fast, an analytic solution would obvio
Dirk Eddelbuettel wrote:
On 15 July 2008 at 17:23, hadley wickham wrote:
| An alternative to enscript is highlight,
| http://www.andre-simon.de/doku/highlight/en/highlight.html, which does
| come with R highlighting built in.
Another alternative is GNU a2ps which has definitions for R source,
do
Here are a couple of ways:
> dat<-rnorm(1000,1,2)
> hist(dat, col='red')
> tmp <- par('usr')
> clip(0,tmp[2],tmp[3],tmp[4])
> hist(dat, col='green', add=TRUE)
>
> # or
>
> library(TeachingDemos)
> hist(dat, col='red')
> clipplot( hist(dat, col='green', add=TRUE), c(0, par('usr')[2]) )
>
Hope this
Hi
[EMAIL PROTECTED] napsal dne 17.07.2008 16:14:33:
> Hi,
> I need to draw a barplot with the value numbers on the bars. In
> particular I would draw a barplot(beside=T) with on the bars the value
> numbers and in the middle of the bars a parcentage.
>
> --
>
Hi
Maybe there is some neat function for it but
apply(data.frame(t(col2rgb(colours()[1:10]))),1,function(x)
paste(as.character.hexmode(x),collapse=""))
gives you hex values for colours specified.
Then you can use colours()[1:10] to get names (or you can invent your own)
and you can print real
On 7/17/2008 10:25 AM, Leandro Marino wrote:
Hi list,
I will help an person that will use some graphics of R in internet.
But this pearson want to specify the colours. This person want me to create
an pallete of colours like that:
Name of color - Code - Colour
White - 0xFF - col
I have found a solution that appears to work. Instead of calling
mpi.close.Rslaves() to shut down the slaves, I use
mpi.bcast.cmd(q("no")) and rely on the scheduler/mpirun to shut down
anything else. However, because I really don't know what I am doing, I
would appreciate it if anyone who sees some
Hi,
I need to draw a barplot with the value numbers on the bars. In
particular I would draw a barplot(beside=T) with on the bars the value
numbers and in the middle of the bars a parcentage.
--
|
12.34% | 150
Dear All:
I am trying to plot a series of data using histogram plot. But I want to
change the default of setting for histogram plot. For example, I want to set
frequency plot starting with zero.
Many Thanks
Xin
[[alternative HTML version deleted]]
_
I know that!
But i will select about 10 colurs and show him
But he needs in that format!
thanks@
-Mensagem original-
De: Duncan Murdoch [mailto:[EMAIL PROTECTED]
Enviada em: quinta-feira, 17 de julho de 2008 11:32
Para: Leandro Marino
Cc: [EMAIL PROTECTED] Org
Assunto: Re: [R] Colo
Dear R users,
I have had a problem using the "step" function for the selection of the
variables from a gamma model with log link function.
CODE:
step.glm.gam1.infroma <- step(glm.gam1.infroma, direction = "backward", k=2)
ERROR MESSAGE:
Error in glm.fit(x[, jj, drop = FALSE], y, wt, offset = obj
Hi list,
I will help an person that will use some graphics of R in internet.
But this pearson want to specify the colours. This person want me to create
an pallete of colours like that:
Name of color - Code - Colour
White - 0xFF - color white (like an box with this color)
I know
Hi,
I would like to compare differences in AUC from 2 different models, glm and gam
for predicting presence / absence. I know that in theory the model with a
higher AUC is better, but what I am interested in is if statistically the
increase in AUC from the glm model to the gam model is signifi
ok, sorry, my mistake was the C printf.
Thank you for your good answer
Regards
- Message d'origine
De : JS Ubei <[EMAIL PROTECTED]>
À : jim holtman <[EMAIL PROTECTED]>
Cc : r-help@r-project.org
Envoyé le : Jeudi, 17 Juillet 2008, 15h25mn 07s
Objet : [R] Re : float and double precision
Dear Steven,
Your xlim() function is what changes the behavior of the x axis values; If you remove the xlim()
statement from your function, you get a "correct" x axis:
##
YC=82:91
Age=rev(seq(2,11,1))
Num=c(2,0,8,21,49,18,79,28,273,175)
box44=data.frame(YC,Age,Num)
mod1=lm(log(Num+1)~YC, dat
thank you for your quick answer,
I'm far of the digits capacity and my values are not the result of a
computation. I'm developping a R package to acces a specific data source. And I
need precision a few better. How can I do ?
When I try this In R console, this is correct and what I need :
>
Dear Hsin-Ya,
The problem here is that subject is a factor with 14 levels, and sequence again is a factor with 2
levels; so the subject:sequence interaction already uses up another (13*1)=13 d.f.; given that there
are only 28 replicates, it is not surprising that there are no residual degrees o
haven't even peeked into the docs, have you?
try:
strsplit(String, ".", fixed=TRUE)
strsplit(String, "[.]")
vQ
Kurapati, Ravichandra (Ravichandra) wrote:
> Hi
>
>
>
> String<-"130.5"
>
>
>
> Df<-Strsplit(".",":",String)
>
>
>
> Then Df get
>
>"" "" "" ""
>
>
FAQ 7.31
Also read
"What Every Computer Scientist Should Know About Floating-Point
Arithmetic", ACM Computing Surveys, 23/1, 5–48, also available via
http://docs.sun.com/source/806-3568/ncg_goldberg.html.
On Thu, Jul 17, 2008 at 8:47 AM, JS Ubei <[EMAIL PROTECTED]> wrote:
> Hi all,
>
> There is
Hi
Str<-"130.1,2.2,3.1,..."
I want the o/p like
Dataframe1
row
130
2
3
.
.
Dataframe2
Col
1
2
1
.
Hi all,
There is a mistake for wich I need your ligths :
///
I have a small C code :
SEXP testData()
{
SEXP result;
void * rans;
float * my_data;
int my_data_length;
my_data_length = 2;
my_data = new float[my_data_length];
m
On Thu, 2008-07-17 at 17:56 +0530, Kurapati, Ravichandra (Ravichandra)
wrote:
> Hi
>
>
>
> String<-"130.5"
>
>
>
> Df<-Strsplit(".",":",String)
>
>
>
> Then Df get
>
>"" "" "" ""
>
> But I want Df should contains
>
>
>
>Df
>
> "130"
>
> "5"
## add horizontal grid:
barchart(data=df,Ratio~reorder(Compound, as.numeric(Class))|Class,
origin = 0, scales = list(x = "free"),
panel=function(...) {
panel.grid(h=-1, v=0); panel.barchart(...)})
## make bars all the same width:
library(latticeExtra)
update(trel
Hi
String<-"130.5"
Df<-Strsplit(".",":",String)
Then Df get
"" "" "" ""
But I want Df should contains
Df
"130"
"5"
If any body knows how to do it.tel me
Thanks
K.Ravichandra
[[alternative HTML version deleted]]
___
Dear Deepayan and other users,
thanks a lot, that was exactly what i intended to get. However, there are
some aspects of the plot that can be improved.
It would be nice to have the same wisth of the bars in all panels. I saw
the bow.width-argument, but as I understood it is no help to get a fix
Hi everyone
It looks like repeated runs of fastICA produce quite significantly different
mixing matrices (not only in terms of sign and row order). I'm not a
specialist, so would appreciate any advice on whether this should really be
the case:
> res3 =
> fastICA(af[,2:20],4,alg.typ="parallel",fu
Barry Rowlingson wrote:
Why use color when you can use black and label the curves where they
are most separated? This solves problems with color blindness,
xeroxing, and faxing.
Where should I put the labels in this example:
> set.seed(123)
> z=matrix(runif(6*50),50,6)
> matplo
On Thu, Jul 17, 2008 at 2:50 AM, Rune Haubo <[EMAIL PROTECTED]> wrote:
> 2008/7/16 Dimitris Rizopoulos <[EMAIL PROTECTED]>:
>> well, for computing the p-value you need to use pchisq() and dchisq() (check
>> ?dchisq for more info). For model fits with a logLik method you can directly
>> use the foll
Kunzler, Andreas wrote:
Dear List,
I ran into some trouble by calculating differences. For me it is
important that differences are either 0 or not.
So I don't understand the outcome of this calculation
865.56-(782.86+0+63.85+18.85+0)
[1] -1.136868e-13
I run R version 2.71 on WinXP
I could so
G'day Christophe,
On Wed, 16 Jul 2008 18:22:49 +0200
[EMAIL PROTECTED] wrote:
> Thanks for your answer.
My pleasure.
> > I guess writing
> > a regular expression that says "export everything that does not
> > start with a dot but do not export foo and bar" would be not
> > trivial to write (at
Hi Andreas,
It's because you are dealing with binary or floating point calculations, not
just a few apples and oranges, or an abacus (which, by the way, is an
excellent calculating device, and still widely used in some [sophisticated]
parts of the world).
http://en.wikipedia.org/wiki/Floating_po
Hi
[EMAIL PROTECTED] napsal dne 17.07.2008 11:47:42:
> Dear List,
>
> I ran into some trouble by calculating differences. For me it is
> important that differences are either 0 or not.
So you shall not use computers. They do not work with infinite precision.
See Faq 7.31
Regards
Petr
>
> So
On Thu, 17 Jul 2008 11:47:42 +0200
"Kunzler, Andreas" <[EMAIL PROTECTED]> wrote:
> I ran into some trouble by calculating differences. For me it is
> important that differences are either 0 or not.
If that is the case, restrict yourself to calculating with integers. :)
> So I don't understand th
Not sure how to interpret standard errors produced by
model.tables(object, type="effects", se=T)
For example, simplest possible case is single stratum anova with two
equally replicated treatments (n observations each). The effects displayed
are plus and minus half the difference between the tw
Dear List,
I ran into some trouble by calculating differences. For me it is
important that differences are either 0 or not.
So I don't understand the outcome of this calculation
865.56-(782.86+0+63.85+18.85+0)
[1] -1.136868e-13
I run R version 2.71 on WinXP
I could solve my problem by using
ro
On Thu, Jul 17, 2008 at 5:13 PM, Mohammad Ehsanul Karim
<[EMAIL PROTECTED]> wrote:
> Dear List,
>
> Say, we generate data like this-
>
> dat<-rnorm(1000,1,2)
> hist(dat)
library(ggplot)
qplot(dat, geom="histogram", colour = factor(dat < 0))
Hadley
--
http://had.co.nz/
Hi,
I am trying to learn lognormal mixture models with EM.
I was wondering how does one compute the log likelihood.
The current implementation I have is as follows,
which perform really bad in learning the mixture models.
__BEGIN__
# compute probably density of lognormal.
dens <- function(lamb
Hi Daniela,
Spencer (? Graves) is not at home. Seriously, this is a list that many
people read and use. If you wish to elicit a response, then you would be
wise to give a better statement of what your difficulty is.
The function you enquire about is well documented with an example, see
##
libra
Here is something that is close:
> x <- rnorm(1)
> y <- hist(x, plot=FALSE)
> plot(y, col=ifelse(y$mid<0,'red','green'))
>
On Thu, Jul 17, 2008 at 5:13 AM, Mohammad Ehsanul Karim
<[EMAIL PROTECTED]> wrote:
> Dear List,
>
> Say, we generate data like this-
>
> dat<-rnorm(1000,1,2)
> hist(dat)
Dear List,
Say, we generate data like this-
dat<-rnorm(1000,1,2)
hist(dat)
How do i make the histogram, say, red (col = 2) before X = dat = 0, and rest
say, green (col = 3) beyond X = dat = 0 in R?
The resulting histogram could be like this
http://ehsan.karim.googlepages.com/histogram.JPG (e
2008/7/16 Dimitris Rizopoulos <[EMAIL PROTECTED]>:
> well, for computing the p-value you need to use pchisq() and dchisq() (check
> ?dchisq for more info). For model fits with a logLik method you can directly
> use the following simple function:
>
> lrt <- function (obj1, obj2) {
>L0 <- logLik(
> Why use color when you can use black and label the curves where they are
> most separated? This solves problems with color blindness, xeroxing, and
> faxing.
>
Where should I put the labels in this example:
> set.seed(123)
> z=matrix(runif(6*50),50,6)
> matplot(z,type='l',col=1,lty=1)
A lab
See ?plot.lm
[EMAIL PROTECTED] wrote:
running R 2.5.1 on WinXP.
plot(model formula here) produces 4 useful residual diagnostic plots. Two questions about them:
1. How can I isolate just one of the four, and have only that one be shown?
str(plot(model)) didn't give me any insight.
2.
Please ask BioConductor related questions on the corresponding mailing list.
In this case: Rgraphviz is an interface package to graphviz. With the
current Windows binary of Rgraphviz you will need to install graphviz
itself as well - in the older version 2.16 (newer versions of graphviz
do not
This one should (I am tempted to write "must") make its way to fortune
()...
Thankyouthankyouthankyou ...
Emmanuel Charpentier
On Mon, 14 Jul 2008 14:58:13 -0600, Greg Snow wrote :
> For those people who feel the need for a p-value to test normality on
>
100 matches
Mail list logo
