Hi Rolf,
On Monday 09 June 2008 11:16:57 pm Rolf Turner wrote:
Your approach tacitly assumes --- as did the poster's question --- that
the probability of passing an item by one method is *independent* of
whether it is passed by the other method. Which makes the methods
effectively
Steven McKinney wrote:
You can either put code such as these first two lines
to grab the names of x and y before you alter x and y
in your function
badplot -
function(x, y, ylab, xlab, ...)
{
if (missing(ylab)) ylab - deparse(substitute(y))
if (missing(xlab)) xlab - deparse(substitute(x))
Dear experts,
I tried to put the two plots into one final PNG figure
with the following script.
However instead of giving 2 plots in one figure,
it only gives the the last plot in one figure.
What's wrong with my script below?
__BEGIN__
in_fname - paste(mydata.txt.,sep=)
out_fname -
Dennis Fisher wrote:
Colleagues
I have just encountered an interesting problem with readLines in R
2.7.0 in Windows Vista. I am trying to read a line that is created in
the following manner:
1. Intel Fortran (ifort) 10.1 creates two text files.
2. The OS concatenates
I am trying to check a package i made (SINFERS) and i get the following error
message in Command Prompt (everything before this is 'OK'):
**checking whether the package can be loadedERROR
Error in library.dynam(lib,package,package.lib):
shared library 'SINFERS' not found
Error:
Hi,
I have the following function that I want to apply to a list of 14
matrices (1536 x 170) of binary data:
DRes - function(x, nr = 1, metric = mixed, ...) {
require(analogue)
require(ade4)
m - c()
for (i in 1:nr) {
set.seed(i)
x1 - x[, sample(dimnames(x)[[2]],
On Mon, Jun 09, 2008 at 02:02:01PM -0700, Ken Spriggs wrote:
Hello,
This code works fine but is so fast I can't see anything but the last plot.
for (i in nrow(X)){
plot(as.numeric(d[i,])) }
I'd like to view a plot every 500 milliseconds, nrow(X) = 400. How?
Adding the line
Thank you very much!
That's exactly what I need as a starting point :-)
Antje
Greg Snow schrieb:
Here is an example that may get you started:
point1 - c(1.1, 1.7, 255)
point2 - c(2.2, 1.5, 180)
point3 - c(1.8, 2.2, 60)
mydf - as.data.frame( rbind(point1, point2, point3) )
names(mydf) -
first preallocate 'm' to the max (m - numeric(nr)) and then run Rprof to
see where time is being spent. Since there was not reproducible data
provided, it is hard to analyze beyond this point. Time is probably being
spent in one of the functions
On Tue, Jun 10, 2008 at 4:49 AM, Marc [EMAIL
Dear All,
I don't know whether this question is appropriate for this list or not. My
question is as follows:
I am fitting a non-linear mixed model using the SSlogis function. When the
parameters xmid and scal are both considered as random (let us say model M1),
the p-value for
1) cv.glm is not 'in R', it is part of contributed package 'boot'. Please
give credit where it is due.
2) There is nothing 'cross' about your 'home-made cross validation'.
cv.glm is support software for a book, so please consult it for the
definition used of cross-validation, or MASS (the
I presume that this is the function from package coda (not mentioned and
we have no sessionInfo() and the 'replicable code' does not work).
On Mon, 9 Jun 2008, Doran, Harold wrote:
Just out of curiosity, why might this be occuring:
class(x6)
[1] mcmc
crosscorr.plot(x6)
NULL
And what
Thanks for the messages. What I probably should have explained is that I have
a data.frame file with variables named A, B, etc.. I wanted to use paste
() to pull the names file$A, file$B, etc. and stick them in table() and so
get an output of table(file$A). For the record, someone read my
Dylan Beaudette wrote:
Hi,
Is anyone on the list familiar with an R implementation of Piper Diagrams?
Example:
http://faculty.uml.edu/nelson_eby/89.315/IMAGES/Figure%209-78.jpg
I am thinking that two calls to triax.plot (plotrix) along with some kind of
affine-transformed standard plot would
Hello
I don't know exactly what you want to do but:
-why do you use in your example h$counts and not h? Furthermore helpl
file says it should be a time series, why then rather not your time series?
-usually na.action will make the default action, which you can see by
getOptions(na.action)
Hello,
I am facing a problem in drawing heat map using R.
I have a 70X3 matrix and I want to draw a heat map with 1 coloumn on X axis
another on Y axis and want to show the value of 3rd coloumn using heat
colours.
Could you please help me with this .
Thank you .
Regards,
Sumit
PD == Peter Dunn [EMAIL PROTECTED]
on Tue, 10 Jun 2008 09:25:07 +1000 writes:
PD Hi all After posting what follows, Duncan Murdoch
PD suggested perhaps a bug in formatC, or an error on
PD documentation. Any comments?
It could be called a bug in your platform's
implementation
Hello everyone,
I am new to R, I have been using SAS for a while. Not surprisingly, I
find R much better in graphics, which is publication ready right away.
Recently, I have been trying to calculate some basic statistics using R.
I have a dataset of multiple rows per subject. For example:
Hi
I consider myself not a complete beginner in R, however an elegant
solution to this problem stumps me. I have a fairly long time series of
6 or so points, I need to gather the data to create a histogram of
the length of continuous zero periods in a data set.
So image the data looks like
Try using aggregate:
aggregate(tb[-1], list(tb$Subject), FUN=median, na.rm = T)
On Tue, Jun 10, 2008 at 10:18 AM, Michal Figurski
[EMAIL PROTECTED] wrote:
Hello everyone,
I am new to R, I have been using SAS for a while. Not surprisingly, I find
R much better in graphics, which is
Thank you Arne. Though I didn't detail all of my attempts, I did in
fact try a larger restriction matrix that accounted for the levels of
the factors.
However, I think my error in those attempts was in adding N columns to
the matrix for an N-level factor (rather than N-1 columns), because I
Hello,
I would be very thankful to hear how often you search for associative
rules in your daily work and what would you consider the most
important in your tasks:
- Data loading speed
- Data processing speed
- Advanced algorithms for rules post-processing
Thanks a lot for your time!
Best
On 6/10/2008 9:23 AM, Creighton, Sean wrote:
Hi
I consider myself not a complete beginner in R, however an elegant
solution to this problem stumps me. I have a fairly long time series of
6 or so points, I need to gather the data to create a histogram of
the length of continuous zero periods
I am tring to use the rollingRegression function from PerformanceAnalytics.
When I apply the function to R stock data (Boston), the function works fine.
But when I try to apply it to data that I imported from excel, using
xlsReadWrite, I recieve an error message:
Error in nrow(x): subscript
I have a data set something like this:
, Value
1972 , 117
1984 , 73
1969 , 92
1976 , 113
1999 , 80
1996 , 78
1976 , 98
1984 , 106
1976 , 99
it could be created with:
dafSamp -
Well,
Yes, that works just fine! Thank you!
--
Michal J. Figurski
HUP, Pathology Laboratory Medicine
Xenobiotics Toxicokinetics Research Laboratory
3400 Spruce St. 7 Maloney
Philadelphia, PA 19104
Henrique Dallazuanna wrote:
Try using aggregate:
aggregate(tb[-1], list(tb$Subject),
You could put the group averages back into dafSamp using ave():
dafSamp - data.frame(cbind(c(1972,1984,1969,1976,1999,1996,1976,1984,1976),
c(117,73,92,113,80,78,98,106,99)))
dafSamp$Ay - ave(dafSamp$X2, dafSamp$X1, FUN=mean)
dafSamp$vecAA - dafSamp$X2 * (dafSamp$Ay /
Hi
[EMAIL PROTECTED] napsal dne 10.06.2008 14:49:55:
Hello,
I am facing a problem in drawing heat map using R.
I have a 70X3 matrix and I want to draw a heat map with 1 coloumn on X
axis
another on Y axis and want to show the value of 3rd coloumn using heat
colours.
Could you please
Try
RSiteSearch(fit a sine)
On Tue, Jun 10, 2008 at 10:25 AM, Jon Loehrke [EMAIL PROTECTED] wrote:
I have been attempting to estimate the periodic contribution of an effect to
some data but have not been able to fit a sine wave within R.
It would be nice to start by being able to fit a sine
Look at the image function and the levelplot function in the lattice
package.
-Original Message-
From: sumit gupta [EMAIL PROTECTED]
To: r-help@r-project.org r-help@r-project.org
Sent: 6/10/08 6:55 AM
Subject: [R] Hello
Hello,
I am facing a problem in drawing heat map using R.
I have
I have been attempting to estimate the periodic contribution of an
effect to some data but have not been able to fit a sine wave within R.
It would be nice to start by being able to fit a sine wave with an
amplitude and frequency.
x-seq(0,20,by=0.5)
y-2*sin(2*pi*.5*x) #amplitude =2,
In MASS there is a nnet() function which returns, among other things, a
wts[] vector giving the weights in the neural network. However, in neither
the book nor in the online documentation for the routine can I find
documentation of which element of wts[] goes with which link in the
Actually, if you have the information on the individual units, then the
production line information is redundant and you can just use Mcnemar's test
comparing QA1 to QA2. Though as Rolf mentioned, this will not give the same
information as if you have true state and can compare false positives
The goal is to look at an evolution of a process over 400 steps. I don't
really need to page back and forth so much as I'd like to just watch it like
a flip book.
jholtman wrote:
Write it out to a PDF file and you can view at your leasure. On WIndows
you
have the option of paging back
So this works... Thanks. My mistake was leaving out seq_len()
for (i in seq_len(nrow(Volumes))){
plot(as.numeric(deltas[i,]))
Sys.sleep(0.5); }
Ken Spriggs wrote:
I appreciate the response but if this is what you have in mind I didn't
get anything different.
Any ideas why?
for (i
Could someone provide a link or examples of parsing XML document in R? Few
specific questions below:
For instance I can retrieve specific nodes using this:
node - xpathApply(xml, // %+% xtag, xmlValue)
1) I want to be able to retrieve parent node for this node, how can I do
this?
In addition to Gabor's suggestion, note the following warning from ?nls
Warning
Do not use nls on artificial zero-residual data.
The nls function uses a relative-offset convergence criterion that
compares the numerical imprecision at the current parameter estimates
to the residual
I believe read.xls has a colClasses argument. If you import using
read.xls(filename, colClasses='character')
everything will be imported as a string, and you can re-convert after
importing
Alberto Monteiro wrote:
I found a curious bug in read.xls. I don't know if it's reproducible.
It's like
Hello,
I have to read in a SPSS-file.
I use spss.get() from the Hmisc-package (version 3.4-3) and win-R(2.7.0)
My Problem now is, that the length of all variable names is cutted to 8
characters. Can I avoid this?
Thank you!
--
Ist Ihr Browser Vista-kompatibel? Jetzt die neuesten
Hello
Using R 2.7.0 on Windows.
I am running a linear discriminant analysis as follows
discrim1 - lda(normvar~ mafmahal+ mrfmahal+ mffmahal+ bafmahal+ brfmahal+
cofmahal+ bmfmahal+ cfmahal+ fractmahal+ antmahal+ absmifmahal+
absifmahal, subset = train)
prediction - predict(discrim1,
Alberto Monteiro wrote:
I found a curious bug in read.xls. I don't know if it's reproducible.
It's like this: suppose I do a read.xls in a spreadsheet. A column
begins with a number. Then, any strings below it will be rendered as NA.
If the column begins with a string, then it will be
Peter Flom wrote:
Hello
Using R 2.7.0 on Windows.
I am running a linear discriminant analysis as follows
discrim1 - lda(normvar~ mafmahal+ mrfmahal+ mffmahal+ bafmahal+ brfmahal+
cofmahal+ bmfmahal+ cfmahal+ fractmahal+ antmahal+ absmifmahal+ absifmahal, subset = train)
prediction
MM == Martin Maechler [EMAIL PROTECTED]
on Tue, 10 Jun 2008 15:05:05 +0200 writes:
PD == Peter Dunn [EMAIL PROTECTED]
on Tue, 10 Jun 2008 09:25:07 +1000 writes:
PD Hi all After posting what follows, Duncan Murdoch
PD suggested perhaps a bug in formatC, or an error on
PD
I found a curious bug in read.xls. I don't know if it's reproducible.
It's like this: suppose I do a read.xls in a spreadsheet. A column
begins with a number. Then, any strings below it will be rendered as NA.
If the column begins with a string, then it will be rendered correctly.
Alberto
I appreciate the response but if this is what you have in mind I didn't get
anything different.
Any ideas why?
for (i in nrow(Volumes)){
plot(as.numeric(deltas[i,]))
Sys.sleep(0.5); }
Thanks
Jan T. Kim wrote:
On Mon, Jun 09, 2008 at 02:02:01PM -0700, Ken Spriggs wrote:
Hello,
This
Hi,
I have the following problem while running a gam regression:
test.gam = gam(fmla,data=Final.Df)
where fmla is a regression formula that looks like: Var1 ~ poly(c(1:160),
degree = 5) + VarC1 + VarC2 + . VarC12
Within SPlus I used the following phrase to find the coefficients
Hello list!
I have a proble trying to perform a SAM analysis using the function samr
from the samr package. I have put the option *center.arrays=TRUE *in order
to scale all the experiments to median=0. I would like to retrieved the
scaled data but it seems that samr does not return it...Does
Hi Pratt --
ppatel3026 [EMAIL PROTECTED] writes:
Could someone provide a link or examples of parsing XML document in R? Few
specific questions below:
Always helpful to know what software you're using; here's mine
library(XML)
sessionInfo()
R version 2.8.0 Under development (unstable)
I wrote (in part)
I am running a linear discriminant analysis as follows
discrim1 - lda(normvar~ mafmahal+ mrfmahal+ mffmahal+ bafmahal+ brfmahal+
cofmahal+ bmfmahal+ cfmahal+ fractmahal+ antmahal+ absmifmahal+
absifmahal, subset = train)
prediction - predict(discrim1,
Hi,
How do I collapse (average in the simplest case) the values of those duplicated
ids (i.e., 2, 5, 6, 9) to give a table of unique ids ?
t - cbind(id=c(1:10, 2,5,6,9), value=rnorm(14))
_
[[alternative HTML
Dear experts,
I tried to put the two plots into one final PNG figure
with the following script.
However instead of giving 2 plots in one figure,
it only gives the the last plot in one figure.
What's wrong with my script below?
__BEGIN__
in_fname - paste(mydata.txt.,sep=)
out_fname -
I wrote
Hello
Using R 2.7.0 on Windows.
I am running a linear discriminant analysis as follows
discrim1 - lda(normvar~ mafmahal+ mrfmahal+ mffmahal+ bafmahal+ brfmahal+
cofmahal+ bmfmahal+ cfmahal+ fractmahal+ antmahal+ absmifmahal+
absifmahal, subset = train)
prediction -
Hi Jim,
This is genotype data of 170 samples. I selected subsets of SNP
optimized for different types of germplasm. So it is a matrix with 170
rows and 1536, 384 or 96 columns of binary data (0, 1). I have 14 of
such matrices in a list.
x - list()
for (i in 1:14) {
set.seed(i)
x[[i]] -
I from a first thought I would say that you are apply this wrong! The
fourier transform convolves a function (cos(x)+isin(x) (this may not be the
exact formula but I don't have my books near)) to the data and then
integrates over -1/2 to 1/2 takes the modulus and plots this- the
periodogram.
Hi All,
I have a problem of putting long titles on a graph:
for example,
x= seq(1:100)
y=seq(1:100)
plot(x,y,main=p=0.05:A-B=3,C-D=10,D-E=100,A-F=2,AFR-E=3,ACE-D=1,ADEF-M=0,AED-E=10,DE-F=3,AB-J=4,AC-J=10,ED-F=1,ED-B=4,AF-B=10,CD-S=10,AM-C=4)
R seems not able to print the whole title. The title
Hi,
I am trying to plot multiple lines on one plot such that all lines are of
the same color (black) and continuous. I need to distinguish these multiple
plots from each other by using different symbols (I am using pch=). However,
all my lines are broken on both sides of the symbols. This
Dear Hua,
Try using \n in the title or the cex.main argument. See ?par for
details.
# First attempt
x= seq(1:100)
y=seq(1:100)
plot(x,y,main=p=0.05:A-B=3,C-D=10,D-E=100,A-F=2,AFR-E=3,ACE-D=1,ADEF-M=0,\n
AED-E=10,DE-F=3,AB-J=4,AC-J=10,\n ED-F=1,ED-B=4,AF-B=10,CD-S=10,AM-C=4)
# Second attempt -
#try this cex.main=1.1 This is a par parameter ?par
plot(x,y,main=p=0.05:A-B=3,C-D=10,D-E=100,A-F=2,AFR-E=3,ACE-D=1,ADEF-M=0,AED-E=10,DE-F=3,AB-J=4,AC-J=10,ED-F=1,ED-B=4,AF-B=10,CD-S=10,AM-C=4,
cex.main=1.1)
On Tue, Jun 10, 2008 at 12:07 PM, Hua Li [EMAIL PROTECTED] wrote:
Hi All,
I have a
?tapply
?aggregate
--
Gregory (Greg) L. Snow Ph.D.
Statistical Data Center
Intermountain Healthcare
[EMAIL PROTECTED]
(801) 408-8111
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of Daren Tan
Sent: Tuesday, June 10, 2008 9:37 AM
To:
Look at ?print.trellis in the lattice package (specifically the split and more
argumnets).
--
Gregory (Greg) L. Snow Ph.D.
Statistical Data Center
Intermountain Healthcare
[EMAIL PROTECTED]
(801) 408-8111
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On
Dear Daren,
Try this:
tapply(t[,value],t[,id],mean)
HTH,
Jorge
On Tue, Jun 10, 2008 at 11:37 AM, Daren Tan [EMAIL PROTECTED] wrote:
Hi,
How do I collapse (average in the simplest case) the values of those
duplicated ids (i.e., 2, 5, 6, 9) to give a table of unique ids ?
t -
Hua Li wrote:
Hi All,
I have a problem of putting long titles on a graph:
for example,
x= seq(1:100)
y=seq(1:100)
plot(x,y,main=p=0.05:A-B=3,C-D=10,D-E=100,A-F=2,AFR-E=3,ACE-D=1,ADEF-M=0,AED-E=10,DE-F=3,AB-J=4,AC-J=10,ED-F=1,ED-B=4,AF-B=10,CD-S=10,AM-C=4)
R seems not able to print the whole
Dear Tariq,
A good starting point would be ?matplot. Also, you can try this:
set.seed(123)
X=cbind(1:10,matrix(rnorm(10*5),ncol=5))
matplot(X[,1],X[,2:5],type='o',lty=1,col=1,pch=1:5,xlab=Variable
X,ylab=Outcome)
and the add your any text you want around.
HTH,
Jorge
On Tue, Jun 10, 2008
In your plot and lines function call you specify type=b which says to use
both lines and points with gaps in the lines around the points. If you
change that to type=o then it will overplot the points and lines (no gaps).
The different types are documented in ?plot.default.
Hope this helps,
Here is one approach:
tmp -
p=0.05:A-B=3,C-D=10,D-E=100,A-F=2,AFR-E=3,ACE-D=1,ADEF-M=0,AED-E=10,DE-F=3,AB-J=4,AC-J=10,ED-F=1,ED-B=4,AF-B=10,CD-S=10,AM-C=4
tmp2 - gsub(',',', ',tmp)
tmp3 - strwrap(tmp2, 30)
par(mar=c(5,4,10,1)+.1)
plot(1:10, main=paste(tmp3, collapse=\n))
Hope this helps,
I have had good sucess with the par(mfrow=c(#,#)) for formating graphs and
they look good to me. I have seen a lot of use of the lattice package and
thought I would go fishing on the list for y'all's comments. Is there a
time when lattice would be easier more appropriate for certain graphics
Doran, Harold HDoran at air.org writes:
Just out of curiosity, why might this be occuring:
class(x6)
[1] mcmc
crosscorr.plot(x6)
NULL
# Replicable code
example(lmer)
x6 - mcmcsamp(fm1, n=1000)
crosscorr.plot(x6)
I cannot reproduce this with the Current Sourceforge version
Dieter
crosscorr.plot is in the coda package, which I should have noted. The
function wasn't producing a chart. It was only returning NULL. But, when
I updated packages, it is now.
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of Dieter Menne
Sent:
Instead of placing this information in the title we could
consider placing it in a legend.
s - a=1,b=2,c=3
plot(1:10)
legend(topleft, legend = strsplit(s, ,)[[1]], ncol = 1, cex = 0.5,
bty = n)
On Tue, Jun 10, 2008 at 12:07 PM, Hua Li [EMAIL PROTECTED] wrote:
Hi All,
I have a problem of
FWIW:
I use lattice exclusively. One difference from the mfrow approach is that
lattice automatically provides (if desired) the same scaling on all plots,
necessary to make visual comparisons. It also provides an easy way to
appropriately label and and arrange the plots, as well as to provide
mfrow is for classic graphics and when the graphs are unrelated.
lattice is for grid graphics and typically when the panels are related
such as being the same plot but conditioned on different factor levels.
These are typical uses only as one can write arbitary panels with lattice
and produce
Could someone please indicate if there is a package capable of computing a
hierarchical two-sided tobit model. Thanks in advance.
---
Anthony A. Pezzola
[EMAIL PROTECTED]
(02) 354-7823
Profesor de Ciencia Política
Instituto de Ciencia Política
Just wanted to say thanks to all suggestions.
My code works fine now *_0
Carl
__
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
Hi all,
I observed that the 'link' between R and BUGS (winBugs/linBugs) has
totally disappeared. I'm wondering what package can be use to run
bayesian models specified using the BUGS language in R (specifically
under Linux). Is there any other option besides JAGS ?
Thanks
Hi all,
I observed that the 'link' between R and BUGS (winBugs/linBugs) has
totally disappeared. I'm wondering what package can be use to run
bayesian models specified using the BUGS language in R (specifically
under Linux). Is there any other option besides JAGS ?
Thanks
DrakeGis wrote:
Hi all,
I observed that the 'link' between R and BUGS (winBugs/linBugs) has
totally disappeared. I'm wondering what package can be use to run
bayesian models specified using the BUGS language in R (specifically
under Linux). Is there any other option besides JAGS ?
Thanks
On 10/06/2008, at 6:08 PM, Ivan Adzhubey wrote:
Hi Rolf,
On Monday 09 June 2008 11:16:57 pm Rolf Turner wrote:
Your approach tacitly assumes --- as did the poster's question ---
that
the probability of passing an item by one method is *independent* of
whether it is passed by the other
Manisha -
Will this work instead?
if 'a' is your matrix, try
ifelse(a == A, 0, 1)
Best,
Erik Iverson
Manisha Brahmachary wrote:
Hello,
I will appreciate any suggestion for this simple problem. I have a matrix of
characters A, P, M.
I want to convert those characters to integers
If you make the x-values more dense, e.g. with x-seq(0,20,0.1) then your
example works. Plotting your data might give a hint to why your example
fails Btw: It is generally not a good idea to use F for a parameter because
F usually means FALSE.
Regards
Søren
Try this:
x - matrix(c(A, P, M, A, P, M), 2)
x[] - as.numeric(as.character(factor(x, labels = c(0,1,1
x
On 6/10/08, Manisha Brahmachary [EMAIL PROTECTED] wrote:
Hello,
I will appreciate any suggestion for this simple problem. I have a matrix
of
characters A, P, M.
I want to
Could you post your code so we can see what you are trying to do?
Thanks,
Patrick
From: [EMAIL PROTECTED] [EMAIL PROTECTED] On Behalf Of Eleni Christodoulou
[EMAIL PROTECTED]
Sent: Tuesday, June 10, 2008 11:20 AM
To: r-help@r-project.org
Subject: [R]
Hi,
I don't have too much programming experience - sorry in advance if this question
is very basic, and thanks for the help.
I have a vector of numbers and another vector of a subset of those numbers. How
can I create a third vector with all the numbers in the original list that
aren't in the
Thanks,
This works fine.
-Original Message-
From: Erik Iverson [mailto:[EMAIL PROTECTED]
Sent: Tuesday, June 10, 2008 4:19 PM
To: Manisha Brahmachary
Cc: r-help@r-project.org
Subject: Re: [R] convert characters into integers in a matrix
Manisha -
Will this work instead?
if 'a' is your
I have this frame in the long format which I want to convert to wide
format. There are three columns in the frame, the id for the rows in
the wide format, the name of the column where the final value (the third
column) is to be placed.
The complicating factor is that the long format is not
WinBugs doesn't run in LINUX. Neither OpenBugs.
DrakeGis wrote:
Hi all,
I observed that the 'link' between R and BUGS (winBugs/linBugs) has
totally disappeared. I'm wondering what package can be use to run
bayesian models specified using the BUGS language in R (specifically
under Linux).
I can not seem to get a ylim argument into plotMeans() is there anyway to do
this.
thanks in advance
Stephen
--
Let's not spend our time and resources thinking about things that are so
little or so large that all they really do for us is puff us up and make us
feel like gods. We are mammals,
The reshape function and even more powerful reshape package will probably do
what you want.
--
Gregory (Greg) L. Snow Ph.D.
Statistical Data Center
Intermountain Healthcare
[EMAIL PROTECTED]
(801) 408-8111
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On
Try plotmeans in the gplots package.
On Tue, Jun 10, 2008 at 5:14 PM, stephen sefick [EMAIL PROTECTED] wrote:
I can not seem to get a ylim argument into plotMeans() is there anyway to do
this.
thanks in advance
Stephen
--
Let's not spend our time and resources thinking about things that
I am using Sweave to generate a document in which ggplot2 figures are
embedded. I am using the following code in a particular plot within a
subfloat environment:
qplotARbosreg, echo = FALSE, fig = TRUE=
print(qplot(bosday2, Arday2, data= arbostimearonlyNArm,
colour=factor(bos2),
I'm tyring to fit a relatively simple nls model to some data, but keep coming
up against the same error (code follows):
Oto=nls(Otolith ~ Linf*(1-exp(-k(AGE-to))),
data = ages,
start = list(Linf=1000, k=0.1, to=0.1),
trace = TRUE)
The error message I keep getting is Error in
Hi
On 11/06/2008 at 10:56 a.m., Ranney, Steven [EMAIL PROTECTED]
wrote:
I'm tyring to fit a relatively simple nls model to some data, but keep coming
up against the same error (code follows):
Oto=nls(Otolith ~ Linf*(1-exp(-k(AGE-to))),
data = ages,
start = list(Linf=1000, k=0.1,
I am using Sweave to generate a document in which ggplot2 figures are
embedded. I am using the following code in a particular plot within a
subfloat environment:
qplotARbosreg, echo = FALSE, fig = TRUE=
print(qplot(bosday2, Arday2, data= arbostimearonlyNArm,
colour=factor(bos2),
Martin
I know you and Brian have been discussing this issue, and Brian
would officially declare it a 'buglet' if he used that word, but
for completeness (I should have done this anyway; my oversight):
What are the exact platform details?
sessionInfo()
R version 2.7.0 (2008-04-22)
Thanks to all for the advice. Took me a bit to get back to this, but
the following worked just fine for me with my 64-bit Ubuntu 8.04 OS:
R CMD INSTALL Rmpi_0.5-5.tar.gz --configure-args=--with-mpi=/usr/lib64/openmpi
On Thu, Jun 5, 2008 at 3:23 AM, Paul Hewson [EMAIL PROTECTED] wrote:
Or (more
Hi,
Are there any algorithms that handle numeric and factor variables
together in a cluster analysis?
Thank you,
Nagu
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide
If you can define a distance between two vectors (where each one has some
numerical and some categorical coordinates) then you can proceed with any
clustering algorithm.
One possibility to get such a distance is to use RandomForest which can produce
a proximity matrix which can be turned into
Hi, All:
Is there a way to get random effects for ARMA parameters?
Consider the following example from the 'corARMA' help page:
fm1Ovar.lme - lme(follicles ~ sin(2*pi*Time) + cos(2*pi*Time),
data = Ovary, random = pdDiag(~sin(2*pi*Time)))
fm5Ovar.lme -
Hi everyone.
I want to solve the following problem. I have a data.frame and I
create a dotplot using lattice.
Then I want to use the grid-package to create a combined graphic which
contains the dotplot as well as a textplot() (using package gplots) of
the data.frame next to the dotplot.
Example
DrakeGis wrote:
WinBugs doesn't run in LINUX. Neither OpenBugs.
There is a long thread ending on
http://tolstoy.newcastle.edu.au/R/e2/devel/07/08/4121.html
which documents how many people helped me get
OpenBUGS running under Linux (embedded in an
R package, of course).
HTH,
Tobias
99 matches
Mail list logo