Dear Justin
Ties means that you have identical values in
"Year5.lm$residuals". Please remark that you can have a large
R^2, but your residuals are not normally distributed. A large
R^2 shows a strong linear relationship, but that does not say
anything about the error distribution (see example belo
On Mon, 2006-07-10 at 13:27 +0800, Guo Wei-Wei wrote:
> Dear all,
>
> I'm working on a data.frame named en.data, which has n cases and m columns.
> I generate the correlation matrix of en.data by
>
> > cor(en.data)
>
> I find that there is no p-value on each correlation in the correlation
> matr
you can use function rcor.test() from package 'ltm', e.g.,
help(rcor.test, package = "ltm")
###
library(ltm)
dat <- data.frame(matrix(rnorm(1000), 100, 10))
rcor.test(dat)
rcor.test(dat, method = "kendall")
rcor.test(dat, method = "spearman")
I hope it helps.
Best,
Dimitris
Dimitris Riz
Dear Chi Kai,
Three years ago there was a similar thread.
At that time David Firth offered a solution for the quasipoission
problem with negative observations, see:
http://finzi.psych.upenn.edu/R/Rhelp02a/archive/16143.html
I remember that his code gave you slightly different answers than the
exa
Hi, Gavin, your program is excellent. Thank you very much!
And I have two further questions.
1. Since it is very possible that the data contains missing value and
the program will failed against missing values, I have to delete all
the cases contained NA. Can it be done pairwisely?
2. Can the pr
On 10-Jul-06 Spencer Graves wrote:
> [...]
Excellent overview of many string reasons for R, Spencer! Many
thanks for this well-thought-out advocacy piece.
> [...]
> At one time, the market dominance of both SAS and IBM was
> sustained by the perception that 'nobody ever got fired'
> for using the
Hello everyone
I am very much a newcomer to R and I need to call R processing routines
from java and pass back the results to java. It seems sjava is the way
to do this, but there seems to be a lot of complexity surrounding
configuring R for this purpose. I need to do this under windows XP. Has
Thank you, Dimitris.
The function rcor.test() is very nice. It can pass arguments to cor()
and solve my first problem of pairwised case deletion.
Best regards,
Wei-Wei
2006/7/10, Dimitris Rizopoulos <[EMAIL PROTECTED]>:
> you can use function rcor.test() from package 'ltm', e.g.,
>
> help(rcor.
Brian Lunergan wrote:
> Uwe Ligges wrote:
Where did you find them? At least not in their current versions on
the official repository on CRAN master for R-2.3.x, I hope.
>>>
>>> FracSim_0.2.zip
>>> http://cran.stat.ucla.edu/bin/windows/contrib/2.3/
>>>
>>> RDCOMClient_0.91-0.zip
>>> http:
Does anyone have any code for drawing weighted histograms (a la
Manet/Mondrian) in R?
Thanks,
Hadley
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/postin
Hi
I'm having difficulty with a multiple plot. What I want is 12
plots, all of which are the same size and shape, differing
only in colour. Each one is a square, so there is an asp=1
in the plot call. I'm working in an Sweave environment so I
am free to choose the height and width of the plot.
(R)SJava is part of the Omegahat project, not the R project. In theory at
least, Omegahat has its own mailing lists.
In any case, the R posting guide asks you to contact the maintainer of a
contributed package before this list, so please do so now.
There are easier-to-setup alternatives to usin
[apologies for possible multiple posting]
Hi Joris
great suggestion! I never thought to use layout() in this way.
If I have
postscript(file="~/f.ps",width=5,height=8)
nf <- layout(matrix(
c(00,01,00,02,00,03,
04,05,00,06,00,07,
Hello all,
I posted a question to this list last week and received no response. I am
unsure if this means no-one knows the answer or if I posed the question badly.
I'm going to assume I posed the question badly and try again. I am new to R so
it is quite likely it's a very naive question, howev
On Mon, 2006-07-10 at 16:22 +0800, Guo Wei-Wei wrote:
> Hi, Gavin, your program is excellent. Thank you very much!
>
> And I have two further questions.
>
> 1. Since it is very possible that the data contains missing value and
> the program will failed against missing values, I have to delete al
Hi,
I am looking forward if some of you could help me with this. I have been
looking R-functions to conduct tree-ring analysis in R. I know that exists
a COFECHA version in R developed by C. Bigler. But I don't know if is there
another collection of functions to do dendroclimatological analysi
On Mon, 2006-07-10 at 11:56 +0100, Richard Birnie wrote:
> Hello all,
Hi Richard,
Sorry, I know nothing about pvclust and have never used it, but here are
a couple of general suggestions/observations.
You are asked to contact the package maintainer *not* R-Help for
questions relating to contribu
Dear All,
Is there a better way to include NA's of a factor in the output of table()
than using as.character()?
Admittedly, I do not understand the help page for table concerning the
exclude argument applied to factors. I tried in different ways, but could
not get NA to be included in the table, i
fcv <- factor(c('a', NA, 'c'), exclude=NULL)
table(fcv, exclude=NULL)
---
Jacques VESLOT
CNRS UMR 8090
I.B.L (2ème étage)
1 rue du Professeur Calmette
B.P. 245
59019 Lille Cedex
Tel : 33 (0)3.20.87.10.44
Fax : 33 (0)3.20.87.10.31
ht
If we specify exclude=NULL in factor then it need not also be
specified in table:
fcv <- factor(c('a', NA, 'c'), exclude=NULL)
table(fcv)
On 7/10/06, Jacques VESLOT <[EMAIL PROTECTED]> wrote:
> fcv <- factor(c('a', NA, 'c'), exclude=NULL)
> table(fcv, exclude=NULL)
> ---
Thanks for the responses,
Now I know why I got no response I'll know better in the future. Apologies all
round for any inconvenience caused that was never my intention, just
inexperience.
Richard
Dr Richard Birnie
Scientific Officer
Section of Pathology and Tumour Biology
Welcome Brenner Build
>
>Hello.
>
>Is it possible to choose the distance in the kmeans algorithm?
>
>I have m vectors of n components and I want to cluster them using kmeans
>algorithm but I want to use the Mahalanobis distance or another distance.
>
>How can I do it in R?
>If I use kmeans, I have no option to choose t
On Mon, 2006-07-10 at 12:48 +0100, Gavin Simpson wrote:
> On Mon, 2006-07-10 at 16:22 +0800, Guo Wei-Wei wrote:
> > Hi, Gavin, your program is excellent. Thank you very much!
> >
> > And I have two further questions.
> >
> > 1. Since it is very possible that the data contains missing value and
>
R 2.3.1
Windows XP
Question: How can I obtain the values of BIAS and STD. ERROR from a
bootstrap.
Background:
I am running a bootstrap:
result2<-boot(1:400,regSEvssample,R=5000)
and obtain the following results:
Bootstrap Statistics :
originalbiasstd. error
t1* 1.876602 -0.00013
Look at the code of boot:::print.boot , which does the calculations.
On Mon, 10 Jul 2006, John Sorkin wrote:
> R 2.3.1
> Windows XP
>
> Question: How can I obtain the values of BIAS and STD. ERROR from a
> bootstrap.
>
> Background:
> I am running a bootstrap:
> result2<-boot(1:400,regSEvssample,
Thank you Jacques and Gabor!
Your solution does work and it led me to try also:
table(factor(fcv, exclude=NULL)) # shows a, c, NA
Greetings,
Heinz
At 08:34 10.07.2006 -0400, Gabor Grothendieck wrote:
>If we specify exclude=NULL in factor then it need not also be
>specified in table:
>
>fcv
Hi,
I have calculated chi-square goodness of fit test,Sample coming from
Poisson distribution.
please copy this script in R & run the script
The R script is as follows
## start
#
No_of_Frauds<-
c(4,1,6,9,9,10,2,4,8,2,3,0,1,2,3,1,3,
I have a data frame with named columns and I would like to know if it is
possible to retrieve a column name once selected:
print(colnames(df)) # assumes to print "col1" "col2"
print.name(df$col1) # would like to print "col1"
print.name(df$col2) # would like to print "col2"
So what the print.n
Hi,
I have calculated chi-square goodness of fit test,Sample coming from
Poisson distribution.
please copy this script in R & run the script The R script is as follows
## start
#
No_of_Frauds<-
c(4,1,6,9,9,10,2,4,8,2,3,0,1,2,3,1,3,
On Mon, 10 Jul 2006, Laurent Deniau wrote:
> I have a data frame with named columns and I would like to know if it is
> possible to retrieve a column name once selected:
Not really. df$col1 is a new object which does not know where it came
from.
If you wanted to do this before selection, then
On Fri, 7 Jul 2006, Valentin Dimitrov wrote:
>
> I do not need a accelerated failure model, but a
> proportional hazard model with a f0= weibull,
> exponential, loglogistic or lognormal baseline
> distribution. The hazard function is
> lambda(t)=exp(Xi*beta)*lambda0(t),
> where lambda0 is the basel
On Sat, 8 Jul 2006, S?ren Merser wrote:
Hi
It seems that survfit() doesn't accept the argumnet 'error' as below
survfit(fit, error='greenwood')
Error in survfit.coxph(fit, error = "greenwood") :
unused argument(s) (error ...)
Isn't is allowed to do that for a coxph object?
Looking
"priti desai" <[EMAIL PROTECTED]> writes:
> Hi,
> I have calculated chi-square goodness of fit test,Sample coming from
> Poisson distribution.
> please copy this script in R & run the script
> The R script is as follows
>
> ## start
>
Use drop = FALSE. For example using builtin data frame BOD we
can display the Time column with its heading:
BOD[, "Time", drop = FALSE]
On 7/10/06, Laurent Deniau <[EMAIL PROTECTED]> wrote:
> I have a data frame with named columns and I would like to know if it is
> possible to retrieve a colum
On Mon, 10 Jul 2006, Gabor Grothendieck wrote:
> The problem can be reduced to this:
>
> x <- 1
> x[1] <<- 2 # error
>
> The following are ok:
>
> x <- 1
> x[1] <- 3
>
> x <- 1
> x <- 4
>
> x <- 1
> x <<- 5
>
> Does anyone know why? Is this a bug in <<- ?
No, it's a feature. The fact that x<<-5
> "roger" == roger koenker <[EMAIL PROTECTED]>
> on Sun, 9 Jul 2006 12:31:16 -0500 writes:
roger> On Jul 9, 2006, at 11:06 AM, Douglas Bates wrote:
>> Your matrix Mc should be flagged as invalid. Martin and
>> I should discuss whether we want to add such a test to
>>
Prof Brian Ripley wrote:
> On Mon, 10 Jul 2006, Laurent Deniau wrote:
>
>> I have a data frame with named columns and I would like to know if it is
>> possible to retrieve a column name once selected:
>
>
> Not really. df$col1 is a new object which does not know where it came
> from.
>
> If y
Hi,
I'm drawing a very simple plot with both axes logarithmic (default base 10).
Example:
vec=c(1,10,100,1000,1,10,100,1000)
plot(vec,vec,log="xy")
The axes on the plot now show the technical notation like 1E+3 but I would
prefer to have it the notation 10 ^3 i.e. with the expone
On 7/10/2006 12:43 PM, [EMAIL PROTECTED] wrote:
> Hi,
>
> I'm drawing a very simple plot with both axes logarithmic (default base 10).
> Example:
> vec=c(1,10,100,1000,1,10,100,1000)
> plot(vec,vec,log="xy")
>
> The axes on the plot now show the technical notation like 1E+3 but I
You can always draw the axes by hand.
?par with axes =FALSE
?axis
?plotmath for mathematical notation in R (for exponents)
-- Bert Gunter
Genentech Non-Clinical Statistics
South San Francisco, CA
"The business of the statistician is to catalyze the scientific learning
process." - George E. P.
With some help from those with expertise on this list, I managed to
produce a plot using trellis that looked like I wanted it to look.
Now, I need to take the same plot and make the lines on it color, but
I want to specify the color for the lines myself.
I've managed to make the key use the colors
Use col= to specify colors, e.g.
library(lattice)
x <- 1:12
xyplot(x ~ x, group = gl(3,4), col = 1:3, type = "l", auto.key = TRUE)
If this is not sufficiently close to your problem
1. cut your example down to a *minimal* size and
2. provide it as *self contained* and
3. *reproducible*
Carlos and Uwe:
I apologize to all, but my server has crashed and I need to order
a new computer! I'm at home now, but when I can get to my office
I'll try to send each of you a copy.
I will be submitting my packages to CRAN soon. Uwe, perhaps the
quad3d stuff should just be incorporated into
The col command worked.
Thanks!
Jamie
On 7/10/06, Gabor Grothendieck <[EMAIL PROTECTED]> wrote:
> Use col= to specify colors, e.g.
>
> library(lattice)
> x <- 1:12
> xyplot(x ~ x, group = gl(3,4), col = 1:3, type = "l", auto.key = TRUE)
>
> If this is not sufficiently close to your p
Try looking at the tkwait.window function, it may be what you need. Create
your tkwindow with all its components, then call tkwait.window on the toplevel
window and the calling function will wait until that window goes away before
continuing execution.
Hope this helps,
--
Gregory (Greg) L.
On 7/10/06, Jamieson Cobleigh <[EMAIL PROTECTED]> wrote:
> With some help from those with expertise on this list, I managed to
> produce a plot using trellis that looked like I wanted it to look.
> Now, I need to take the same plot and make the lines on it color, but
> I want to specify the color f
There are a couple of different options that you can try (each
generalize a bit differently to other situations).
If you just want lines at the tick marks to go the entire width (and be
generated automatically to match the labels), then look at ?par and look
at 'tck', you can use this in the barp
I am a very novice R user, a social scientist (linguist) who is trying
to learn to use R after being very familiar with SPSS. Please be kind!
My concern:
I cannot figure out a way to get an accurate count of observations of
one column of data split by a factor when there are NAs in the data.
I kn
Wouldn't something like table(status) give you want you want? E.g.:
R> status <- factor(c("A", "B", "A", NA, "A", "B"))
R> table(status)
status
A B
3 2
Andy
From: Jenifer Larson-Hall
>
> I am a very novice R user, a social scientist (linguist) who
> is trying to learn to use R after being
Dear Dr. Graves,
Thanks for your information and kindness! I really appreciate it.
According to my findings so far, POLYCLASS doesnot seem to allow for
variables being forced into the model. Even the commercialized MARS
doesnot have this function either. However, one ad hoc approach that
they sug
"Jenifer Larson-Hall" <[EMAIL PROTECTED]> writes:
> I am a very novice R user, a social scientist (linguist) who is trying
> to learn to use R after being very familiar with SPSS. Please be kind!
>
> My concern:
> I cannot figure out a way to get an accurate count of observations of
> one column
Hi Folks,
I'm trying to locate the source code for a (typical)
R graphics device, in order to study how it's done.
The underlying reason is that I'm thinking of trying
to create a graphics device for 'pic' (the diagram
drawing component of [g]troff).
I thought the xfig device would be a good pla
Dear R Users,
Is it possible to make R print the largest item in each row of a matrix X
with red font? Example:
1247
8431
...
Therefore 7 and 8 should be in red color.
I would appreciate any suggestion
Robert McFadden
That's correct. The option is not there.
You can force it in a starting model - but it can still be deleted,
At some stage Insightful implemented a commercial set of polynomial
spline functions in S-Plus that does have that option. They got it
as far as a beta-library [quoting from my own homepage
Hello,
I am new to R and still feeling my way thru it.
I am trying to plot the values from this file below on the X-axis of a plot. I
have attached the graph to the email...the one i am trying to recreate.
Exonstart end
5'UTR 2254006022540121
1 225401222254
On 7/10/2006 3:52 PM, (Ted Harding) wrote:
> Hi Folks,
>
> I'm trying to locate the source code for a (typical)
> R graphics device, in order to study how it's done.
>
> The underlying reason is that I'm thinking of trying
> to create a graphics device for 'pic' (the diagram
> drawing component o
Dear Gang,
I'm not clear about the form of the model that you're trying to fit,
nor what the form of the likelihood would be for such as model, but the
short answer is that sem() requires that each covariance among observed
variables be for the same observations and hence for the same number
(N) o
On 10-Jul-06 Ted Harding wrote:
> Hi Folks,
>
> I'm trying to locate the source code for a (typical)
> R graphics device, in order to study how it's done.
>
> [...]
>
> When, in R, I type "xfig" to see the R code, I get the
> line
>
> .Internal(XFig(file, old$paper, old$family, old$bg, old$fg
One option is
library(R2HTML)
?HTML.cormat
The thing you're after is traffic highlighting (via CSS or HTML tags).
If HTML.cormat() doesn't do exactly what you want, modify the source
code. (By the way, I haven't used R2HTML so far so maybe there's a
more appropriate function.)
> -Original Mes
[EMAIL PROTECTED] wrote:
> Hi,
>
> I'm drawing a very simple plot with both axes logarithmic (default base 10).
> Example:
> vec=c(1,10,100,1000,1,10,100,1000)
> plot(vec,vec,log="xy")
>
> The axes on the plot now show the technical notation like 1E+3 but I would
> prefer to have
dear R wizards: sorry, I am stumped. what is the parameter to just
move the location of the x-label and y-lable (not the labels on the
ticks of the x-axis and y-axis)? probably obvious, but not to me
right now... regards, /iaw
__
R-help@stat.math.eth
Does this do what you want:
x <- "Exonstart end
5'UTR 2254006022540121
1 2254012222540140
2 2254030322540493
3 2254155222541565
4 2254237322542519
5 2254426522544432
3'UTR 2254443322544856"
y
Hello there,
This question is relative to WindowsXP, using R 2.2.1:
I am looking for a function that allows a user to interactively choose a
directory so I can use list.files to process all the files in that
directory. I've looked at getwd, but this is not interactive. The
functions file.choo
Hi,
for testing if the slope of experimental data differs from a
given slope I'm using the function
"test_regression_against_slope" (see below).
I am now confronted with the problem that I have data which
requires a modelII regression (also called reduced major axes
regression (RMA) or geometric
choose.dir
it is documented on the same page as choose.files
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
On 7/10/2006 8:51 PM, Randy Zelick wrote:
> Hello there,
>
> This question is relative to WindowsXP, using R 2.2.1:
Time to upgrade. choose.dir (mentioned by Rich) was introduced in the
next release.
Duncan Murdoch
>
> I am looking for a function that allows a user to interactively choose a
dear R wizards: is it possible to instruct R to save or no-save from
inside R? or does this have to be given at invokation on the
command-line?The same question applies to "--no-restore-data",
although this presumably would have to be decided in a .First()
function or something like it.
on a
Thanks to all who responded. I upgraded to 2.3.1 straight away, and
choose.dir does exactly what I need.
Cheers,
=Randy=
R. Zelick email: [EMAIL PROTECTED]
Department of Biology voice: 503-725-3086
Portland State University fax: 5
ivo welch wrote:
> dear R wizards: is it possible to instruct R to save or no-save from
> inside R?
You can specify the action when you call q() or quit(), but you can't
change the default action.
> or does this have to be given at invokation on the
> command-line?The same question applies
Hi there,
I tried to install gstat package to R in Linux. I follow the
instruction to use command " R CMD INSTALL -1 lib pkgs". But R
shows "syntax error in "R CMD"".
Can someone help me with the installation of gstat package?
Thank you!
Jing
__
R-help@
Dear Prof. Kooperberg,
Thanks a lot for your clarification! That is very helpful. I will check
out the polynomial
spline functions in S-Plus that you mentioned.
Best regards,
-XG Su
Xiaogang Su, Assistant Professor
Department of Statistics and Actuarial Scienc
Hi,
How can I recover the variance-covariance matrix of the tobit model from
the variance-covariance of the survreg?
I first used to the survreg function and then I selected the variance
matrix. However, the last parameter is log(scale) and not the variance
of the standard deviation of the cens
On Mon, 10 Jul 2006, ivo welch wrote:
> dear R wizards: sorry, I am stumped. what is the parameter to just
> move the location of the x-label and y-lable (not the labels on the
> ticks of the x-axis and y-axis)? probably obvious, but not to me
> right now... regards, /iaw
They are all done by
It is -l (minus ell) not -1 (minus one)
Note that without actual cut-and-paste of the output, we cannot be sure if
this is the problem: this is why the posting guide asks for such details.
On Mon, 10 Jul 2006, Jing Gao wrote:
> Hi there,
> I tried to install gstat package to R in Linux. I foll
-
Find out what India is talking about on Yahoo! Answers India.
--- Begin Message ---
Subject: Other models of GARCH
To: R HELP
MIME-Version: 1.0
Content-Type: multipart/alternative; boundary="0-1656081585-1147428333=:42419"
Con
Hi R-Users,
First of all i apologize if this is too simple for you. I want to plot an index
of which i have daily data for 10 years. I want to plot specific dates in X
axis, like 20/05/91 20/12/92 20/07/94 etc. and the index value in
Y axis. Please suggest me.
Thanks,
Sumanta.
Geoffrey Matthews a écrit :
> I will be submitting my packages to CRAN soon. Uwe, perhaps the
> quad3d stuff should just be incorporated into klaR?
Hello Geoffrey,
I read your paper at
http://www.r-project.org/useR-2006/Presentations/
http://www.r-project.org/useR-2006/Slides/Matthews.pdf
I was
77 matches
Mail list logo