Hi, I get the above error when starting my program with rake (ruby make)
instead of manual startup (which works as expected).
It happens on R initialization with "REmbeddedPostgres --gui=none
--silent" as arguments (this is part of my way of embedding R in a C
application using libR.so). Any hel
On 3/3/08, Richard Rowe <[EMAIL PROTECTED]> wrote:
> This is a homework problem. I know how to do a PCA, you need to learn.
> I suggest you visit your textbook, then check the documentation for R's
> various PCA implementations to work out how to effect the analysis.
Check Rcmdr. There you can
Dear Shu,
Why not store your dataset in a database? Then you can start each loop
by reading the submatrix you need for the analysis. This will require
much less memory. loops from the apply-family with work better than the
for loop.
HTH,
Thierry
-
On 3/3/08, shu zhang <[EMAIL PROTECTED]> wrote:
> Hello R users,
>
> I'm wondering whether it is possible to manage big data set in R? I
This [1] recent thread might be of interest.
Liviu
[1] http://www.nabble.com/How-to-read-HUGE-data-sets--tt15729830.html
_
On Sunday 02 March 2008 (18:50:53), Luca Braglia wrote:
> hello everybody
>
> i was seeking a converter beetween R regression output (eg with
> summary) and the conventional way to present regression output in
> paper: every model as a vertical vector with \beta, t beetween
> parenthesis below the
Dear All,
Is there any existing code for plotting silhouette for kmeans clustering
results?
Many thanks!
Linda
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read
On 3/3/08 8:21 PM, Ericka Lundström wrote:
> I'm trying to emigrate from SPSS to R, thou I have some problems whit
> getting R to distinguish between the different kind of missing.
...
> Is there a smart way in R to differentiate between missing and valid
> and at the same time treat both the
On Mon, 3 Mar 2008, [EMAIL PROTECTED] wrote:
> On 03-Mar-08 03:19:01, Wensui Liu wrote:
>> HI, John,
>> my understanding is that you should use log(...) instead of its
>> original scale. Below is the logic in the case of poisson reg.
>> log(y / offset) = x'b
>> => log(y) - log(offset) = x'b
>> =>
> Calum wrote:
>
>> All of that is very nice so far. The I followed bits and pieces of
>> other peoples posts in the past to plot on a weibull regression...
>>
>> > my_curve.Plac <- survreg( Surv(Survival, Censored==0)~
>> TreatmentGroup, subset=TreatmentGroup=="Placebo", data=TestData,
>>
Hi,
I consulted ?png, and it uses X11. is there any way to save plots into png,
without using X11 ?
On 3/3/08, Prof Brian Ripley <[EMAIL PROTECTED]> wrote:
>
> Please consult the help page for png.
>
> On Mon, 3 Mar 2008, Ng Stanley wrote:
>
> > Hi,
> >
> > I have installed R on a computational
Does
names(obj)
do what you want?
Patrick Burns
[EMAIL PROTECTED]
+44 (0)20 8525 0696
http://www.burns-stat.com
(home of S Poetry and "A Guide for the Unwilling S User")
Nair, Murlidharan T wrote:
>Is there a method to list the components of an object, instead of looking at
>the help for tha
On 3/3/08 6:42 AM, ArunPrasad wrote:
> Hi everyone,
> I am looking for a package in R which can help me in using the
> imputation technique to find the missing values for my regression analysis.
> Any help would be appreciated.
>
> Cheers
> Arun
Have a look at the "Missing data" se
Dear Sir/Madam,
Am using mt.teststat function from 'multtest' package. My variable to
'classlabel' argument having 'NA'.I tried to get rid of this by using
na=.mt.naNUM. But this parameter is not working. I understood from help file
that this haven't fully implemented yet. Could yo suggest
I have a number of non-periodic signals , that have been detrended and denoised.
Each signal is made up of a variable number of cycles non necessarily
of the same length.
We can now fit a regression model to each cycle and then use the
returned coefficients,
together with some other data, to form c
Andreas Maunz wrote:
> Hi, I get the above error when starting my program with rake (ruby make)
> instead of manual startup (which works as expected).
> It happens on R initialization with "REmbeddedPostgres --gui=none
> --silent" as arguments (this is part of my way of embedding R in a C
> ap
Dear R users,
I have a problem since I try to plot my datas with different colors.
plot(tvar, var, xlab="zeit [s]",ylab="Variation [%]", col = ifelse(var <=
varstability, 'green','red'))
this works well!
But since I add a type="l" to my plot, it will color all the plot with
green!!!
Is there an
lin tang wrote:
> hi, Dear R users -
>
> I wonder is there a written R function to convert a number to a text, say
> convert 1 to "one" , 100 to "one hundred". I know in xls. has such a
> function BAHTTEXT, does anybody know is there a similar function in R ?
> Thanks.
>
Hi Lin,
I got bored a
Uwe Ligges wrote:
> Andreas Maunz wrote:
>> Hi, I get the above error when starting my program with rake (ruby
>> make) instead of manual startup (which works as expected).
>> It happens on R initialization with "REmbeddedPostgres --gui=none
>> --silent" as arguments (this is part of my way of em
Hi there!
I'd like to make many tests and plots, for each variables.
i.E.
i <- 1
repeat {
print (t.test(e[i], f[i]))
i <- i + 1
if (i == 21) break
}
but the output is:
Welch Two Sample t-test
data: e[i] and f[i]
...
so I don't know which data was calculated.
Is there a way
Hi
you can also
write.table(probe_gene, file = "c:/foo.xls", sep = "\t", col.names = NA)
or
write.table(tab, "clipboard", sep = "\t" , col.names = NA)
and after opening Excel just press Ctrl-V. Maybe you also find usefull to
add row.names=FALSE into export construction.
Petr
[EMAIL PROTECTED
The other respondants here have missed the point - write.table() as
invoked by the OP does not produce a csv file. The default separator is
" " (a space) in write.table(), so no wonder Excel cocked it up as it
assumes that a csv file is one that has elements separated by a comma
",". We don't need
Does anybody know how one can implement the Somer's and Goodman -Kruskal
ordinal measures of association in R?
With regards
-
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
One option is:
> head(f)
A B C D E
1 0.75135094 -0.06779700 -0.2617316 -0.3701998 0.7332249
2 -0.85657366 0.27442139 -0.4291410 0.1107827 0.6942490
3 1.16125212 -0.18631979 -0.3918744 0.4804644 0.9374566
4 0.32100928 0.03773391 -0.8712896
Andreas Maunz wrote:
> Hi, I get the above error when starting my program with rake (ruby make)
> instead of manual startup (which works as expected).
> It happens on R initialization with "REmbeddedPostgres --gui=none
> --silent" as arguments (this is part of my way of embedding R in a C
> appl
On Mon, Mar 03, 2008 at 02:03:07AM -0800, mysimbaa wrote:
>
> Dear R users,
> I have a problem since I try to plot my datas with different colors.
>
> plot(tvar, var, xlab="zeit [s]",ylab="Variation [%]", col = ifelse(var <=
> varstability, 'green','red'))
> this works well!
>
> But since I add
You want just the points with colors?
var <- rnorm(10)
vars <- rnorm(10)
plot(var, type='l')
points(var, col=ifelse(var < vars, 'green', 'red'), pch= 16)
or
the lines also?
for(i in 1:9)segments(i, var[i], i + 1, var[i+1], ifelse(var[i] <
vars[i], 'green', 'red'))
On 03/03/2008, mysimbaa <[E
>Does anybody know how one can implement the Somer's and Goodman -Kruskal
>ordinal measures of association in R?
install.packages("Hmisc")
library(Hmisc)
?rcorr.cens
HTH,
Tobias
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinf
> "LS" == Linda Smith <[EMAIL PROTECTED]>
> on Mon, 3 Mar 2008 00:57:49 -0800 writes:
LS> Dear All,
LS> Is there any existing code for plotting silhouette for kmeans clustering
LS> results?
Yes. Have you tried at all?
Anyway, here's some commented example code
(with a la
Jim Lemon wrote:
>
> I got bored after the millions, but this should handle the smaller
> numbers, and you can always extend it.
>
> digits2text<-function(x,mult="") {
> units<-c("one","two","three","four","five",
>"six","seven","eight","nine")
> teens<-c("ten","eleven","twelve","thirtee
check also the article by John Fox in Rnews volume 5/1, May 2005,
Programmer's Niche 51-55:
http://cran.r-project.org/doc/Rnews/Rnews_2005-1.pdf
Best,
Dimitris
Dimitris Rizopoulos
Biostatistical Centre
School of Public Health
Catholic University of Leuven
Address: Kapucijnenvoer 35, Leuve
Try
RSiteSearch("numbers2words")
On Sun, Mar 2, 2008 at 9:05 PM, lin tang <[EMAIL PROTECTED]> wrote:
> hi, Dear R users -
>
> I wonder is there a written R function to convert a number to a text, say
> convert 1 to "one" , 100 to "one hundred". I know in xls. has such a
> function BAHTTEXT, d
On Fri, Feb 29, 2008 at 11:12 PM, Louise Hoffman
<[EMAIL PROTECTED]> wrote:
> [snip]
>
> > Seriously. Be specific if you have a problem. (read the posting guide). R
> > can
> > also plot. If you don't like R's plots (which I could not understand) you
> can
> > export data and import them to
On Mon, Mar 3, 2008 at 1:25 AM, Ng Stanley <[EMAIL PROTECTED]> wrote:
> Hi,
>
> I consulted ?png, and it uses X11. is there any way to save plots into png,
> without using X11 ?
See "See also" under help(png) for alternatives. Rule of thumb: If
you get a reply from BR that you don't get the fir
zack holden:
> I need to sort through a vector (x) and identify the point at which 2
> successive values become smaller than the previous value.
x <- c(5,5,7,6,5,4,3)
a=c(diff(x, 1) < 0, FALSE) & c(diff(x, 2) < 0, FALSE, FALSE)
a # FALSE FALSE TRUE TRUE TRUE FALSE FALSE
which(a) # 3 4 5
Hi R-users,
Maybe my last question was not specified enough to get a lot of answers,
so I'm trying to do a new attempt:
I use R 2.6.2 and Windows XP. I've got the following question: I wrote a
lot of text into the margin of a plot.
I have to take the borders of the graphics device and enlarge th
Dear James and Arun,
The Social Sciences task view is out of date (updating it is on my to-do
list), and I'm aware of at least one other package, Amelia, for multiple
imputation of missing data (which isn't mentioned in the more up-to-date
Multivariate task view). There are probably others as well
> "LH" == Louise Hoffman <[EMAIL PROTECTED]>
> on Sat, 1 Mar 2008 00:54:56 +0100 writes:
>> If you still want to then read ?write.table, that can export your data
>> into a spreadsheet-like ascii format which can be used from GNUplot
>> easily.
LH> Very interesting.
James Reilly wrote:
> On 3/3/08 8:21 PM, Ericka Lundström wrote:
> > I'm trying to emigrate from SPSS to R, thou I have some problems whit
> > getting R to distinguish between the different kind of missing.
> ...
> > Is there a smart way in R to differentiate between missing and valid
> > and a
Hi,
anybody know any way to make a bash script to start Rcmdr directly widthout
need to open R and execute library(Rcmdr)?
I try to make this using R CMD BATCH somethink, but dont find the way.
Thanks
Ronaldo
--
> Prof. Ronaldo Reis Júnior
| .''`. UNIMONTES/Depto. Biologia Geral/Lab. de Biolog
Dear list,
I am trying to build a cox model and then perform ROC analysis in order to
retrieve some genes that are correlated with breast cancer. When I calculate
the hazard score taking into account different numbers of genes and their
coefficients ( I am trying to find the pest predictor number
On Friday 29 February 2008 (14:50:53), Silvia Lipski wrote:
> Dear R-users,
>
> I am sorry if I ask for something that has been asked
> before, however, I still could not solve my little
> problem by consulting the previous thread on this
> topic:
>
> I would like to replace several values in a dat
Your code below changes the class of the time index from "Date" to "character".
Perhaps you want "yearmon" class. The following averages all values in
each month producing a series with time class "yearmon". See ?yearmon
library(zoo)
z <- zoo(c(0.1423065, -0.03276228, -0.0652789, -0.04999873,
-0
Utkarsh Singhal wrote:
> Hi R,
>
>
>
> I am getting this error while trying to use 'lrm' function with nine
> independent variables:
>
>
>
>> res =
> lrm(y1994~WC08301+WC08376+WC08316+WC08311+WC01001+WC08221+WC08106+WC0810
> 1+WC08231,data=y)
>
>
>
> singular information matrix in lrm.f
Dear Ronaldo,
You should be able to put
local({
old <- getOption("defaultPackages")
options(defaultPackages = c(old, "Rcmdr"))
})
in an .Rprofile file in a directory reserved for this purpose, and have your
bash script start R from that directory (untested). See ?Startup for details
on how R sta
On Mon, 3 Mar 2008, Ronaldo Reis Junior wrote:
> Hi,
>
> anybody know any way to make a bash script to start Rcmdr directly widthout
> need to open R and execute library(Rcmdr)?
>
> I try to make this using R CMD BATCH somethink, but dont find the way.
For some reason, Rcmdr requires an interacti
Have a look at the Hmisc package
Kind regards,
Marco
-Oorspronkelijk bericht-
Van: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED]
Namens John Fox
Verzonden: maandag 3 maart 2008 14:05
Aan: 'James Reilly'; 'ArunPrasad'
CC: r-help@r-project.org
Onderwerp: Re: [R] Imputation Packages
Dear J
Attached is a histogram chart with a break on y-axis which my friend made
for me using phot shop. Do anyone know how to make such a break using R?
Thanks in advance.
(See attached file: 2007 age distribution.doc)
Hongsheng (Hank) Liao, Ph.D.
Lab Manager
Center for Quantitative Fisheries Ecology
As someone of partly French heritage, I would also ask how this
distribution came to be called "Gaussian". It seems very unfair to de
Moivre, who discovered the distribution at least half a century earlier.
:-)
--Jim Rogers
On Mar 2, 2008, at 7:33 AM, (Ted Harding) wrote:
> Hi Folks,
> Apo
Thanks Martin
Well it works except that "as" seems to not like the "initialize" method
: the following code (that is the same than yours with some initialize
for A B and C) does not compile. It seems that as(c,"A") does not work
if we definie a initialize for A...
--- 8< --
setClas
On 3/3/2008 9:10 AM, Rogers, James A [PGRD Groton] wrote:
> As someone of partly French heritage, I would also ask how this
> distribution came to be called "Gaussian". It seems very unfair to de
> Moivre, who discovered the distribution at least half a century earlier.
> :-)
Just an example of St
Well, here is the discussion I had with Duncan in the meantime...
Anybody else a suggestion what might be the problem?
Greetings,
Andreas
--
http://www.maunz.de
Yoda of Borg are we: Futile is resistance. Assimilate you, we will.
--- Begin Message ---
FFrom [EMAIL PROTECTED] Mon Mar 3 15
On Mon, 3 Mar 2008, Martin Maechler wrote:
>> "LH" == Louise Hoffman <[EMAIL PROTECTED]>
>> on Sat, 1 Mar 2008 00:54:56 +0100 writes:
>
>>> If you still want to then read ?write.table, that can export your data
>>> into a spreadsheet-like ascii format which can be used from GNU
Calum wrote:
>>> Also is it possible to get an R-squared type value for the fit of
>>> this curve from someplace?
>>>
>>> Finally (three questions in one!) the first two censored data points
>>> (1 in each group) are actually lost to follow-ups. Should they be
>>> marked differently from censored?
It is possible to do this with R, but I really don't think it is
desirable because it can produce a very misleading graphic. Why do
you want to introduce a break in the axis?
Hadley
On Mon, Mar 3, 2008 at 8:03 AM, HongSheng Liao <[EMAIL PROTECTED]> wrote:
>
> Attached is a histogram chart with
Calum had a long question about drawing survival curves after fitting a Weibull
model, using pweibull, which I have not reproduced.
It is easier to get survival curves using the predict function. Here is a
simple example:
> library(survival)
> tfit <- survreg(Surv(time, status) ~ factor(ph.ecog
On Mon, Mar 03, 2008 at 10:13:53AM -0300, Ronaldo Reis Junior wrote:
> Hi,
>
> anybody know any way to make a bash script to start Rcmdr directly widthout
> need to open R and execute library(Rcmdr)?
>
> I try to make this using R CMD BATCH somethink, but dont find the way.
$ r -lRcmdr -e'while
On Mon, 03 Mar 2008 22:02:17 +1300, James Reilly wrote
> On 3/3/08 8:21 PM, Ericka Lundström wrote:
> > I'm trying to emigrate from SPSS to R, thou I have some
> problems whit > getting R to distinguish between the different
> kind of missing. ... > Is there a smart way in R to
> differentiate
On Mon, Mar 3, 2008 at 4:03 AM, mysimbaa <[EMAIL PROTECTED]> wrote:
>
> Dear R users,
> I have a problem since I try to plot my datas with different colors.
>
> plot(tvar, var, xlab="zeit [s]",ylab="Variation [%]", col = ifelse(var <=
> varstability, 'green','red'))
> this works well!
>
> Bu
Dear ALL:
How I can find the number of observations less than each value in column one
with indicator 1 in column two. Please see the data below.
For example: number of observations less than 1 with indicator 1 (including
those 1 with indicator 1) =2
number of observations
Dear Carlos,
One approach is to use structural equation modeling (SEM). Some SEM
packages, such as LISREL, Mplus and Mx, allow inequality and nonlinear
constraints. Phantom variables (Rindskopf, 1984) may be used to impose
inequality constraints. Your model is basically:
y = b0 + b1*b1*x1 + b2*b2*
> "BDR" == Prof Brian Ripley <[EMAIL PROTECTED]>
> on Mon, 3 Mar 2008 14:16:15 + (GMT) writes:
[]
BDR> In a later message Louise mentioned the desire to use TeX fonts for
BDR> annotation, to match a LaTeX document. Paul Murrell has pointed out
his
BDR>
Dear all,
I've been trying to install the lme4 package from
http://r-forge.r-project.org/projects/lme4/. However, when wanting to load
the package, I get an message saying that "x-tabs" are masked (see pasted
code). Can anyone point to what has gone wrong?
Kind regards,
Andreas Nord
Sweden
>
Nothing's wrong. It just means that the package or one of its
dependencies, has its own xtabs function, which "hides" the default
xtabs function, since it comes first in the search path. So when you
next write xtabs(...), it is this new xtabs that is being loaded. If
you want to call the or
Dear Ellison:
it did not do it.
I edited my previous email to make my question more clear.
The out put should be: (2,11,33,43,46)
For example:
number of all observations less than 1 with indicator 1 (including those 1 with
indicator 1 but not 1 with indicator 0) =2
num
On Mar 3, 2008, at 11:12 AM, stephen sefick wrote:
> x<-read.zoo("SC2.csv", sep="," , format="%m%m/%d%d/%y%y%y%y %h%h:%m%
> m")
>
> #Error in read.zoo("SC2.csv", sep = ",", format =
> "%m%m/%d%d/%y%y%y%y %h%h:%m%m") :
> index contains NAs Error message
You need header=TRUE in th
> Workshop on user interfaces and interactive graphics in R
> ==
> Tuesday, 1st April 10.00-16.30 Postgraduate Statistics Centre,
> Department of Maths and Statistics, Lancaster University, UK
> =
> This one-day work
Ciao a tutti!
ho un problema con un esercizio in R.
In una lista devo trovare tutte le parole che si differenziano per una
lettera da un'altra parola.
Esempio: per "casa": cosa, rosa...
C'è una formula generale da applicare a tutte le parole?
Grazie dell'aiuto, Flo
-
One must specify that there is a header if there is one, each % item
should only
be listed once, and format is currently only used for "Date" class so
you want to
use the FUN= argument -- although in the development version of zoo format is
availailable more generally. Here are several ways to do t
For me is not very clear, but if I understand:
sapply(sort(unique(data[data[,2]==1,1])),
function(x)sum(data[data[,2]==1 & data[,1] <= x, 1]))
But the output is:
2 6 31 71 86
On 03/03/2008, AbouEl-Makarim Aboueissa <[EMAIL PROTECTED]> wrote:
> Dear Ellison:
>
> it did not do i
__
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 provide commented, minimal, self-contained, reproducible code.
Hello,
I am trying to make a barplot with nested count data which is build like
this: first there are several birds (n)laying 3 clutches composed of 2 eggs
half of the second and third clutch received treatment and this treatment
was tested to influence sex of offspring. I want a barplot showing c
Dear ALL:
Please see below. I hope this will make it more clear.
[1,]11
[2,]11 number of all observations less than 1 with
indicator 1 (including those 1 with indicator 1 but not 1 with
indicator 0)=2
[3,]10
[4,]10
[5,]10
[6,]10
[7,
Hi Everyone,
I need some simple help.
Here are my codes
##will give me 1 probesets
data.sub = data.matrix[order(variableprobe,decreasing=TRUE),][1:1,]
dim(data.sub)
data_output<-write.table(data.sub, file = "c://data_output.csv", sep = ",",
col.names = NA)
perhaps you need something like sum(X[X < 1])
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
-Original Message-
From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On
Behalf Of AbouEl-Makarim Aboueissa
Sent: Monday, M
Hi Christophe --
This is a variant of the problem that Jim Regetz is having in a thread
in R-devel. Here's where the trouble is
> as(c, "A")
Error in .local(.Object, ...) :
argument "value" is missing, with no default
By default, 'as(c, "A")' will create a new instance of it's second argument
I believe I wrote too hastily and that what you want is sum(X < 1) which
will sum the indicator (T/F) function.
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
-Original Message-
From: [EMAIL PROTECTED] [mailto:[EMAIL
G'day Carlos,
On Mon, Mar 3, 2008 at 11:52 AM
Carlos Alzola <[EMAIL PROTECTED]> wrote:
> I am trying to get information on how to fit a linear regression
> with constrained parameters. Specifically, I have 8 predictors ,
> their coeffiecients should all be non-negative and add up to 1. I
> unde
You can try this:
cbind(data.sub, p.value=apply(data.sub, 1, function(x)t.test(x)$p.value))
On 03/03/2008, Keizer_71 <[EMAIL PROTECTED]> wrote:
>
> Hi Everyone,
>
> I need some simple help.
>
> Here are my codes
>
> ##will give me 1 probesets
> data.sub = data
Hi, there,
I cannot get accurate value for calculation.
for example:
ld<-sqrt(1*0.05*0.95*0.05*0.95)
0.05*0.95-ld=-6.938894e-18
0.05*0.95-ld==0 is False.
I met this problem in my program, how can I handle it. Thanks.
xj.
__
R-help@r-project.org mai
HI, Thanks for your help previously. I can go around R and scilab now, except
import my research data.
This is an example dataset (below) of radon radiation levels. How can I
import this txt file using Rcmdr?
How do I import .xls files? My alpha guard radon monitors output .dvd files
that can b
ICANN 2008EXTENDED SUBMISSION DEADLINE MARCH 10 2008
Dear Colleague
Due to numerous requests we decided to extend the submission deadline
for
FAQ 7.31 (You need to understand what floating point numbers are)
On 3/3/08, Xuejun Qin <[EMAIL PROTECTED]> wrote:
> Hi, there,
> I cannot get accurate value for calculation.
> for example:
> ld<-sqrt(1*0.05*0.95*0.05*0.95)
> 0.05*0.95-ld=-6.938894e-18
> 0.05*0.95-ld==0 is False.
>
> I met this
On Mon, 2008-03-03 at 12:52 -0500, Xuejun Qin wrote:
> Hi, there,
> I cannot get accurate value for calculation.
> for example:
> ld<-sqrt(1*0.05*0.95*0.05*0.95)
> 0.05*0.95-ld=-6.938894e-18
> 0.05*0.95-ld==0 is False.
>
> I met this problem in my program, how can I handle it. Thanks.
Answer 1
Xuejun Qin wrote:
> Hi, there,
> I cannot get accurate value for calculation.
> for example:
> ld<-sqrt(1*0.05*0.95*0.05*0.95)
> 0.05*0.95-ld=-6.938894e-18
> 0.05*0.95-ld==0 is False.
>
> I met this problem in my program, how can I handle it. Thanks.
>
>
> xj.
>
> ___
On Mon, Mar 3, 2008 at 8:25 AM, Duncan Murdoch <[EMAIL PROTECTED]> wrote:
> On 3/3/2008 9:10 AM, Rogers, James A [PGRD Groton] wrote:
> > As someone of partly French heritage, I would also ask how this
> > distribution came to be called "Gaussian". It seems very unfair to de
> > Moivre, who disc
If I understand you correctly what you want to do is
do t-test (mu=0) for each column of the data.
Treating the data as a data.frame rather than a matrix
you can do something like this and then pick out the
p-values but with 140 t-tests I don't know what you'll
get in terms of anything meaninful.
R FAQ 7.31. G.
On Mon, Mar 03, 2008 at 12:52:43PM -0500, Xuejun Qin wrote:
> Hi, there,
> I cannot get accurate value for calculation.
> for example:
> ld<-sqrt(1*0.05*0.95*0.05*0.95)
> 0.05*0.95-ld=-6.938894e-18
> 0.05*0.95-ld==0 is False.
>
> I met this problem in my program, how can I handl
Hey Flo,
I'm not entirely sure what you're problem is -- due to the language thing,
probably -- but I think you want a function that will give you (or computes
in an intermediate step) the "edit distance" between strings. There is the
"sdists" function in the cba package that might help (it can g
apparently you want to check the genefilter package...
it defines functions like:
rowttests
colttests
rowFtests
colFtests
rowVars
rowSds
moreover, a quick look at Biobase is recommended...
that would save you lots of time as you wouldn't have to reinvent the
wheel.
b
On Mar 3, 2008, at
Yes, the sociologist Robert Merton.
url:www.econ.uiuc.edu/~rogerRoger Koenker
email[EMAIL PROTECTED]Department of Economics
vox: 217-333-4558University of Illinois
fax: 217-244-6678Champaign, IL 61820
On Mar 3, 2008, at 12
On Mon, Mar 3, 2008 at 9:57 AM, Charilaos Skiadas <[EMAIL PROTECTED]> wrote:
> Nothing's wrong. It just means that the package or one of its
> dependencies, has its own xtabs function, which "hides" the default
> xtabs function, since it comes first in the search path. So when you
> next write x
Dear Doug,
As I recall, according to Stigler, yes -- he wasn't the first to
formulate Stigler's law of eponymy (but I don't recall to whom he
attributed it).
Regards,
John
On Mon, 3 Mar 2008 12:17:59 -0600
"Douglas Bates" <[EMAIL PROTECTED]> wrote:
> On Mon, Mar 3, 2008 at 8:25 AM, Duncan Mur
Hi R,
I am getting this error while trying to use 'lrm' function with nine
independent variables:
> res =
lrm(y1994~WC08301+WC08376+WC08316+WC08311+WC01001+WC08221+WC08106+WC0810
1+WC08231,data=y)
singular information matrix in lrm.fit (rank= 8 ). Offending
variable(s):
WC08101 WC0822
Hi,
Is there any package which provides the functions of create one
dimensional and/or Two dimensional classifiers?
Thanks much.
--
Waverley @ Palo Alto
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read
On 03/03/08 12:52, Xuejun Qin wrote:
> Hi, there,
> I cannot get accurate value for calculation.
> for example:
> ld<-sqrt(1*0.05*0.95*0.05*0.95)
> 0.05*0.95-ld=-6.938894e-18
> 0.05*0.95-ld==0 is False.
>
> I met this problem in my program, how can I handle it. Thanks.
I think what you are ex
--- John Fox <[EMAIL PROTECTED]> wrote:
> Dear Doug,
>
> As I recall, according to Stigler, yes -- he wasn't
> the first to
> formulate Stigler's law of eponymy (but I don't
> recall to whom he
> attributed it).
Possibly a disgruntles M. de Moivre?
>
> Regards,
> John
>
> On Mon, 3 Mar 2
The gp.plot function has a type option that can either be 'p' for points
or 'l' for lines, the function is not that complicated, you could easily
add additional options.
I don't know why your code is not creating the pdf file, you could try
plotting to the screen first (don't change the output and
Thanks for all response. Here is more closer to my question,
p<-0.05
f<-0.05
ld<-sqrt(f* (1 - f) * p * (1 - p))
D <- (p * (1 - f) - ld)/p
d <- ((1 - p) * f - ld)/(1 - p)
haplo<-data.frame(D,d)
> haplo
D d
1 -1.387779e-16 -7.304099e-18
cond1<-pmin(haplo[,1],haplo[,
There was some discussion a while back of an R pgf driver
for latex. Is anyone working on that?
On Mon, Mar 3, 2008 at 2:09 PM, Greg Snow <[EMAIL PROTECTED]> wrote:
> The gp.plot function has a type option that can either be 'p' for points
> or 'l' for lines, the function is not that complicated,
Douglas Bates wrote:
>On Mon, Mar 3, 2008 at 8:25 AM, Duncan Murdoch <[EMAIL PROTECTED]> wrote:
>
>
>>On 3/3/2008 9:10 AM, Rogers, James A [PGRD Groton] wrote:
>> > As someone of partly French heritage, I would also ask how this
>> > distribution came to be called "Gaussian". It seems very unfai
1 - 100 of 159 matches
Mail list logo
