Hi Spencer
Thank you for your reply. I tried as you shad suggested and it seems
to me that problem comes from this piece of code
contr <- object$contrasts
Browse[1]>
debug: for (i in names(dataMix)) {
if (inherits(dataMix[, i], "factor") && !is.null(contr[[i]])) {
levs <- levels(dat
You're right, that's exactly what I needed! Thanks!
Rehceb
On Thu, 2006-09-07 at 14:47 +0200, Dimitris Rizopoulos wrote:
> probably you're looking for a barplot, e.g.,
>
> v <- c(1, 1, 2, 2, 2, 3, 3, 4, 4, 4)
> plot(factor(v))
>
>
> I hope it helps.
>
> Best,
> Dimitris
>
>
> Dimitris R
Could somebody program this kind of plot type to R, if none exists,
based on mds or correlation tables or some more suitable method? What
do you think about idea? Does it work? None similar or better exists?
http://weightedassociationmap.blogspot.com/
Atte Tenkanen
University of Turku, Finla
If you don't need borders on the polygons then it can be simply done two points
at a time checking that neither point is an NA:
# data
x1 <- x2 <- 1:8
y1 <- c(1,5,6,1,4,5,5,5)
y2 <- c(0,3,3,NA,NA,1,3,4)
# plot
plot(x1,y1,type="l")
lines(x2,y2)
# fill in area between curves with green two points
Issue this command and then click anywhere on the plot.
loc <- locator(1); do.call(text, c(loc, "abc"))
On 9/7/06, Nair, Murlidharan T <[EMAIL PROTECTED]> wrote:
> I am trying to add text at specific location on my graph. I know this
> can be done in R but I can't recollect.
>
> I was trying usi
Hello, dear useR,
Is there anyways to get the posterior probabilites for each group by using
knn() instead of only get the proportions of winning class?
Example knn(train=Train[,-c(1:3)], test=Test, cl=group.id.train,k=K, prob=True)
will give you the proportions for winning votes, but I also wan
"MARK LEEDS" <[EMAIL PROTECTED]> writes:
> thanks. i assumed we we were talking about the standard textbook difference
> between the t test and pairwise t test.
> my bad.
Notice the difference between paired and pairwise...
> - Original Message -
> From: "Chuck Cleland" <[EMAIL PROTEC
thanks. i assumed we we were talking about the standard textbook difference
between the t test and pairwise t test.
my bad.
- Original Message -
From: "Chuck Cleland" <[EMAIL PROTECTED]>
To: "MARK LEEDS" <[EMAIL PROTECTED]>
Cc: "Li,Qinghong,ST.LOUIS,Molecular Biology" <[EMAIL PROTECTED
On Thu, 7 Sep 2006, Jonathan Boyd Thayn wrote:
> I am trying to install the rgdal package on my Mac OS X 3.9. I
> downloaded and installed the GDAL libraries from Fink and then tried to
> install rgdal and got the following message. I tried to determine if
> the GDAL libraries were in my path
MARK LEEDS wrote:
> no, because the formula for the test statistics ( even assuming that
> variances are equal ) of the two different tests are different. in the
> pairwise t test, the pairwise differences are
> viewed as one sample so it turns into a one sample test. any intro stat book
> will
no, because the formula for the test statistics ( even assuming that
variances are equal ) of the two different tests are different. in the
pairwise t test, the pairwise differences are
viewed as one sample so it turns into a one sample test. any intro stat book
will have the formulas.
All,
I'm trying to create an augPred plot in the nlme library, similar to the
plot on
p.43 of Pinheiro & Bates (Mixed Effects Models in S and S-Plus) for
their Pixel data.
My data structure is the same as the example but I still get the error
msg below.
> comp.adj.UKV <- groupedData(adj.UKV
I am trying to add text at specific location on my graph. I know this
can be done in R but I can't recollect.
I was trying using locator() to identify the position and use identify()
but I can get it to work. Can someone jog my memory?
Thanks ../Murli
[[alternative HTML version dele
Are there functions to read image files in jpg, gif or even a pdf file?
Thanks ../Murli
[[alternative HTML version deleted]]
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the post
Hi,
If I set the p.adjust="none", does it meant that the output p values from the
pairwise.t.test will be the same as those from individual t.tests (set
var.equal=T, alternative="t")?
I actually got different p values from the two tests. See below. Is it supposed
to be this way?
Thanks
Johnny
try something like the following:
lis1 <- c(lapply(1:10, rnorm, n = 10))
lis2 <- c(lapply(1:10, rnorm, n = 10))
lis1[[5]] <- lis2[[8]] <- numeric(0)
ind <- sapply(lis1, length) > 0 & sapply(lis2, length) > 0
lis1 <- lis1[ind]
lis2 <- lis2[ind]
mapply(f
Hi Andrew,
Thanks a lot, That would give me what I want.
But using my own data and models resulted in this:
> plot(fitted(tcos31.c.cp, level=1), FCR.c$g.cp)
Error in xy.coords(x, y, xlabel, ylabel, log) :
'x' and 'y' lengths differ
This is quite correct, as there are some missing values
Dear all,
I'm a newbie to R and I would really apperciate any help with the following:
I have two lists, l1 and l2:
l1:
$"A*0101"
[1] 0.076 0.109 0.155 0.077 0.09 0 0 0.073
[9] 0.33 0.0034 0.0053
$"A*0247"
[1] 0 0 0.5 .004 0 0 0
$"A*0248"
[1] 0 0 0.3 0 0.06
l2:
$"A*1101"
[1] 0.17
Using R in batch mode should work on both Windows and Linux:
R CMD BATCH (assuming that R.exe is in your path)
Even without R's location in your path, you could issue the following
command at the prompt (in windows):
"c:\Program Files\R\R-2.3.1\bin\R.exe" CMD BATCH --vanilla --slave
i:\R_HOME\bat
On 7 September 2006 at 14:39, Zodet, Marc W. (AHRQ) wrote:
| Is there any way to run an R script file (i.e., *.R) from the command
| prompt in the console window. Ultimately, I'm looking to put such code
| in a script file so that it can set off other R scripts/programs as
| needed.
Which platfo
One obvious alternative is an SQL join, which you could do directly in
a DBMS, or from R via RMySQL / RSQLite /... Keep in mind that creating
indexes on user/userid before the join may save a lot of time.
> -Original Message-
> From: [EMAIL PROTECTED]
> [mailto:[EMAIL PROTECTED] On Behalf
XLSolutions Corporation (www.xlsolutions-corp.com) is proud to
announce our 2-day October 2006 "R/S-plus Fundamentals and Programming
Techniques" : www.xlsolutions-corp.com/Rfund.htm
*** Washington DC / October 12-13, 2006
*** Seattle Wa / October 19-20
*** San Francisco / October 26-27
Re
On Thu, 7 Sep 2006, Prof Brian Ripley wrote:
> Which version of R?
Previously, 2.3.1.
> Please try 2.4.0 alpha, as it has a different and more efficient
> algorithm for the case of 1-1 matches.
I downloaded and installed R-latest, but got the same error message:
Error: cannot allocate vector
All:
Is there any way to run an R script file (i.e., *.R) from the command
prompt in the console window. Ultimately, I'm looking to put such code
in a script file so that it can set off other R scripts/programs as
needed.
Thanks.
Marc
[[alternative HTML version deleted]]
I'm working with an historical image that may be (one of?) the first
uses of gray-scale shading to show the pattern of values in a
matrix/table, later used by Bertin in his 'reorderable matrix'
and sometimes called a scalogram.
The image is at
http://euclid.psych.yorku.ca/SCS/Gallery/images/Priv
On Thu, 2006-09-07 at 11:57 -0500, Erin Hodgess wrote:
> Dear R People:
>
> Way Off Topic:
>
> Is anyone aware of a website that contains important dates
> in statistics history, please?
>
> Maybe a sort of "This Day in Statistics", please?
>
> I thought that my students might get a kick out of
I am trying to install the rgdal package on my Mac OS X 3.9. I
downloaded and installed the GDAL libraries from Fink and then tried to
install rgdal and got the following message. I tried to determine if
the GDAL libraries were in my path but I'm not sure how to do that.
Any ideas? Thanks.
On Thu, 2006-09-07 at 12:14 -0500, Hao Chen wrote:
> Hello Marc Schwartz
>
> On Thu, Sep 07, 2006 at 07:54:05AM -0500, Marc Schwartz wrote:
> > On Thu, 2006-09-07 at 06:18 -0500, Hao Chen wrote:
> > > Hi,
> > >
> > > I am using barplot and would like to know if it is possible to have bars
> > > f
Hello Marc Schwartz
On Thu, Sep 07, 2006 at 07:54:05AM -0500, Marc Schwartz wrote:
> On Thu, 2006-09-07 at 06:18 -0500, Hao Chen wrote:
> > Hi,
> >
> > I am using barplot and would like to know if it is possible to have bars
> > filled with one color while use a different color for the shading li
Dear R People:
Way Off Topic:
Is anyone aware of a website that contains important dates
in statistics history, please?
Maybe a sort of "This Day in Statistics", please?
I thought that my students might get a kick out of that.
(actually I will probably enjoy it more than them!)
Thanks for any
Hello,
I want to colorize the area between two curves, but one of these
curves isn't continuous.
The best solution I found is the 2nd example in the help of polygon,
but how can I get no area filling for the missing data in the 2nd curve.
example:
x1 = c(1:8)
x2 = c(1:8)
y1 = c(1,5,6,1,4,5,5,5)
[EMAIL PROTECTED] asked:
>I am trying to divide the columns of a matrix by the first row in the
>matrix.
Dividing columns of a matrix by a vector is a pretty fundamental
operation, and the query resulted in a large number of suggestions:
x/matrix(v, nrow(x), ncol(x), byrow = TRUE))
sweep(x, 2, v
Hi, Patrick: Thanks very much. I'll try it. Spencer Graves
Patrick Burns wrote:
> This is a very common computation in finance.
>
> On the public domain page of the Burns Statistics website
> in the financial part is the code and R help file for
> 'factor.model.stat'. Most of the complication
On 07 Sep 2006 17:20:29 +0200, Peter Dalgaard <[EMAIL PROTECTED]> wrote:
> Martin Maechler <[EMAIL PROTECTED]> writes:
>
> > > "DB" == Douglas Bates <[EMAIL PROTECTED]>
> > > on Thu, 7 Sep 2006 07:59:58 -0500 writes:
> >
> > DB> Thanks for your summary, Hank.
> > DB> On 9/7/06,
Martin Maechler <[EMAIL PROTECTED]> writes:
> > "DB" == Douglas Bates <[EMAIL PROTECTED]>
> > on Thu, 7 Sep 2006 07:59:58 -0500 writes:
>
> DB> Thanks for your summary, Hank.
> DB> On 9/7/06, Martin Henry H. Stevens <[EMAIL PROTECTED]> wrote:
> >> Dear lmer-ers,
> >> M
On 9/7/06, Martin Maechler <[EMAIL PROTECTED]> wrote:
> > "DB" == Douglas Bates <[EMAIL PROTECTED]>
> > on Thu, 7 Sep 2006 07:59:58 -0500 writes:
>
> DB> Thanks for your summary, Hank.
> DB> On 9/7/06, Martin Henry H. Stevens <[EMAIL PROTECTED]> wrote:
> >> Dear lmer-ers,
>
> "DB" == Douglas Bates <[EMAIL PROTECTED]>
> on Thu, 7 Sep 2006 07:59:58 -0500 writes:
DB> Thanks for your summary, Hank.
DB> On 9/7/06, Martin Henry H. Stevens <[EMAIL PROTECTED]> wrote:
>> Dear lmer-ers,
>> My thanks for all of you who are sharing your trials and tri
On Thu, 7 Sep 2006, alex lam (RI) wrote:
> Dear list,
>
> I have been trying to run the function "qvalue" under the package qvalue
> on a vector with about 20 million values.
>
> > asso_p.qvalue<-qvalue(asso_p.vector)
> Error: cannot allocate vector of size 156513 Kb
> > sessionInfo()
> Version
Dear Dr Bates,
Many thanks for such a useful response to my problem.
Regarding Variance Components .
The VarCorr function runs fine for lmer objects once the nlme package is
removed.
Regarding the format of the nested random effects for an lmer object, you
said:
In recent versi
Dear list,
I have been trying to run the function "qvalue" under the package qvalue
on a vector with about 20 million values.
> asso_p.qvalue<-qvalue(asso_p.vector)
Error: cannot allocate vector of size 156513 Kb
> sessionInfo()
Version 2.3.1 (2006-06-01)
i686-pc-linux-gnu
attached base packages
Joris De Wolf wrote:
> Are your sure your second solution does not work? Try again...
>
>
Turns out the second approach did work - but only once I stopped
cutting-and-pasting between two different operating systems (Linux and
Windows under Linux). Apparently, some of the cut-and-paste things
Dear Douglas,
> I would be happy to re-institute p-values for fixed effects in the
> summary and anova methods for lmer objects using a denominator
> degrees of freedom based on the trace of the hat matrix or the rank
> of Z:X
Please do!
> if others will volunteer to respond to the "these answer
Hi,
I am currently analysising a counting process form of a cox model allowing for
the inclusion of time dependent covariates. An example model I have fitted is
modlqol<-coxph(Surv(Tstart,Tstop,cens.time)~tmt.first+risk
+lqol+cluster(id),data=cat)
summary(modlqol)
My question is quick. I am
> I would like to join repeated measures for patients across two visits using
> a line. The program below uses symbols to represent each patient. Basically,
> I would like to join each pair of symbols.
This is easy in ggplot:
install.packages("ggplot")
library(ggplot)
qplot(visit, var, id=patien
Thanks for your summary, Hank.
On 9/7/06, Martin Henry H. Stevens <[EMAIL PROTECTED]> wrote:
> Dear lmer-ers,
> My thanks for all of you who are sharing your trials and tribulations
> publicly.
> I was hoping to elicit some feedback on my thoughts on denominator
> degrees of freedom for F ratios
On Thu, 2006-09-07 at 06:18 -0500, Hao Chen wrote:
> Hi,
>
> I am using barplot and would like to know if it is possible to have bars
> filled with one color while use a different color for the shading lines.
>
> The following code colors the shading lines, leaving the bars in white:
>
> barpl
On Thu, 2006-09-07 at 14:35 +0200, Rotkiv, Rehceb wrote:
> Hello everyone,
>
> I would be glad if you could help out an R-beginner here... I have a
> vector of categorial data like this
>
> > v <- c(1, 1, 2, 2, 2, 3, 3, 4, 4, 4)
>
> When I do
>
> > hist(v)
>
> I get the x-axis of the histogram
Some of the coordinates might not match and also I do
not have the same number of observations in every
table but I want to get only the common ones back.
This is where it gets tricky!I have tried merge, scan
and every joining function I could find but nothing
seems to do what I want.
the R-sig-geo
probably you're looking for a barplot, e.g.,
v <- c(1, 1, 2, 2, 2, 3, 3, 4, 4, 4)
plot(factor(v))
I hope it helps.
Best,
Dimitris
Dimitris Rizopoulos
Ph.D. Student
Biostatistical Centre
School of Public Health
Catholic University of Leuven
Address: Kapucijnenvoer 35, Leuven, Belgium
Tel:
Hello everyone,
I would be glad if you could help out an R-beginner here... I have a
vector of categorial data like this
> v <- c(1, 1, 2, 2, 2, 3, 3, 4, 4, 4)
When I do
> hist(v)
I get the x-axis of the histogram with floating point labels: 1.0, 1.5,
2.0, etc. Is it possible to tell R that th
You can use the 'zoo' or 'its' packages. For 'zoo' see the
documents listed at the end of:
http://cran.r-project.org/src/contrib/Descriptions/zoo.html
On 9/6/06, [EMAIL PROTECTED] <[EMAIL PROTECTED]> wrote:
>
> hi all
>
> i have date and the return series like below, but the dates are not in
> u
On Thu, 7 Sep 2006, isidora k wrote:
> Hi everyone!
> I have 100 tables of the form:
> XCOORD,YCOORD,OBSERVATION
> 27.47500,42.52641,177
> 27.48788,42.52641,177
> 27.50075,42.52641,179
> 27.51362,42.52641,178
> 27.52650,42.52641,180
> 27.53937,42.52641,178
> 27.55225,42.52641,181
> 27.56512,42.526
Hi everyone!
I have 100 tables of the form:
XCOORD,YCOORD,OBSERVATION
27.47500,42.52641,177
27.48788,42.52641,177
27.50075,42.52641,179
27.51362,42.52641,178
27.52650,42.52641,180
27.53937,42.52641,178
27.55225,42.52641,181
27.56512,42.52641,177
27.57800,42.52641,181
27.59087,42.52641,181
27.60375,
try this:
do.call(cbind, l)
Best,
Dimitris
Dimitris Rizopoulos
Ph.D. Student
Biostatistical Centre
School of Public Health
Catholic University of Leuven
Address: Kapucijnenvoer 35, Leuven, Belgium
Tel: +32/(0)16/336899
Fax: +32/(0)16/337015
Web: http://med.kuleuven.be/biostat/
http:/
Dear list,
I have a list of data.frames (generated by "by"), that I want to "stack" into a
single data.frame.
I can do this by cbind, but only by subsetting the list explicitly like this:
cbind(l[[1]],l[[2]],l[[3]],l[[4]])
I find this ugly and not very general.
I tried
cbind(l
Dear lmer-ers,
My thanks for all of you who are sharing your trials and tribulations
publicly.
I was hoping to elicit some feedback on my thoughts on denominator
degrees of freedom for F ratios in mixed models. These thoughts and
practices result from my reading of previous postings by Doug
Ernst O Ahlberg Helgee wrote:
> Hi!
> Im sorry to bother you but I cant fix this.
> I use the lattice function levelplot and I want the colorkey at the
> bottom, how do I get it there? I have tried changing colorkey.space and
> changing in legend but I cant get it right, plz help
>
> btw I'd l
Hi,
I am using barplot and would like to know if it is possible to have bars
filled with one color while use a different color for the shading lines.
The following code colors the shading lines, leaving the bars in white:
barplot(1:5, col=c(1:5), density=c(1:5)*5)
while the colors are applied
Hi!
Im sorry to bother you but I cant fix this.
I use the lattice function levelplot and I want the colorkey at the
bottom, how do I get it there? I have tried changing colorkey.space and
changing in legend but I cant get it right, plz help
btw I'd like to speceify strings to appear at the tick
Hi CG,
I think that the best pair of summary plots are
1) the fitted values without random effects against the observed
response variable, and
2) fitted values with random effects against the observed response
variable.
The first plot gives a summary of the overall quality of the fixed
e
Dear all.
R 2.3.1, W2k.
I am working with a field trial series where, for the moment, I do
regressions using more than one covariate to explain the protein levels
in malting barley.
To do this I use lme() and a mixed call, structured by both experiment
(trial) and repetition in each experimen
This is a very common computation in finance.
On the public domain page of the Burns Statistics website
in the financial part is the code and R help file for
'factor.model.stat'. Most of the complication of the code
is to deal with missing values.
Patrick Burns
[EMAIL PROTECTED]
+44 (0)20 8525 0
Which version of R?
Please try 2.4.0 alpha, as it has a different and more efficient
algorithm for the case of 1-1 matches.
On Wed, 6 Sep 2006, Adam D. I. Kramer wrote:
> Hello,
>
> I am trying to merge two very large data sets, via
>
> pubbounds.prof <-
> merge(x=pubbounds,y=prof,by.x="user"
Version 0.4.4 of odfWeave is available from CRAN. A Windows binary
should be available shortly.
This version requires base R version 2.3.1 or greater.
Changes from the last version include
- Non-English character sets are handled better. For example, Chinese
characters can be included in R cod
hi all
i have date and the return series like below, but the dates are not in
uniform intervals. Please show me the way how to create a time series in
'R' so that dates are also associated with the returns.
thanks in advance
Sayonara With Smile & With Warm Regards :-)
G a u r a v Y a
65 matches
Mail list logo