what is wrong?
Your data is ihmo not correct imported
myArray - read.table(c:/mydata.csv,sep=,,na.strings=)
...and it seems you have using quotes in your csv file so you should
play with the
additional parameter quotes in read.table.
Type ?read.table
Then check your data i.e. with
See ?xyplot and its argument legend.
Uwe Ligges
Ravi Varadhan wrote:
Hi,
How can I add legends in the xyplot function, in the lattice library?
Here is a simulation example:
x - runif(90)
z - sample(1:3, 90, rep=T)
y - rnorm(90, mean = x^2 + z, sd=1)
Hi,
I'm looking for software that can perform kriging on systems with
dimensionality higher than 3, say d=5.
Are anyone aware of packages in R that can do this?
Thanks,
Eivind Smørgrav
---
The information contained in this
Hi,
Fields might be what you are looking for.
Christian
2006/1/17, Eivind Smørgrav [EMAIL PROTECTED]:
Hi,
I'm looking for software that can perform kriging on systems with
dimensionality higher than 3, say d=5.
Are anyone aware of packages in R that can do this?
Thanks,
Eivind
Hi
the BACCO bundle includes package emulator, which is effectively
kriging in arbitrary dimensions.
best wishes
Robin
On 17 Jan 2006, at 09:45, Eivind Smørgrav wrote:
Hi,
I'm looking for software that can perform kriging on systems with
dimensionality higher than 3, say d=5.
Are
John Logsdon j.logsdon at quantex-research.com writes:
I am trying to extract the solution from a simple lme calculation.
y-c(-1.118,-.5,.5,1.118,10)
gp-factor(c(rep('one',4),'two'))
res-lme(y~1,rand=~1|gp)
...
Random effects:
Formula: ~1 | gp
(Intercept) Residual
Thank you very much for your help but I think there is an error for the
answer to the first problem I spent time on searching the solution but
I failed to find it. I tried to put which.max instead of which.min
but it doesn't work. I tried to do my best but i didn't have any idea to
solve this
Hi all,
In R, it is not possible to set the font size of axis labels directly
(AFAIK). Instead, scaling factors for the font chosen by the graphics
device can be supplied. It appears that there is no constant font size
for axis labels. My impression is that the axis label font size is
scaled
I'd like to fit a GLM to each of a number of subsets of some data. The `family'
argument to `lmList' (in lme4) has given me cause for optimism, but so far I've
only been able to achieve linear model fits. For example
df - data.frame(gp = gp.temp - factor(rep(1:3, each = 100)),
x = x.temp -
sorry, i let toto$Date in the function within lapply() instead of x$Date !
now, it works :
toto
Num Date Place X Y
1 1 1/1/04 0:48x1 1 1
2 2 1/1/04 8:02x1 NA NA
3 4 1/1/04 1:55x4 3 7
4 3 1/1/04 2:14x3 2 9
5 4 1/1/04 1:19x4 3 7
6 4
Hello all,
I am using RSPerl package to display plots in R . I want to be able to compare
some data in boxplot and in histogram and in probplot ( from e1071 package)
- ( cumulative probability plot - normal)
I have a nested arrray of data ... from perl like the following
data[0][0] =
From: Wensui Liu [EMAIL PROTECTED]
Just curious how xlispstat is used in the industry and what's it strengthen
compared with other computing languages such as R or matlab?
Almost not at all, though there are a few holdouts.
On a related note, I've been doing some interesting things with a
On Tue, 17 Jan 2006, Hilmar Berger wrote:
In R, it is not possible to set the font size of axis labels directly
(AFAIK).
Nor of anything else, since you are limited to the fonts available on the
output device.
Instead, scaling factors for the font chosen by the graphics
device can be
Why does this happen when I do a sqlQuery?
sessionInfo()
Version 2.3.0 Under development (unstable) (2006-01-04 r36984)
powerpc-apple-darwin8.3.0
attached base packages:
[1] methods stats graphics grDevices utils datasets
[7] base
other attached packages:
RODBC
1.1-5
b -
I have dowloaded the Source Code of R,and I want to know the process of
chi-sqared test,but how can I found it?
[[alternative HTML version deleted]]
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE
On 1/17/2006 8:33 AM, [EMAIL PROTECTED] wrote:
I have dowloaded the Source Code of R,and I want to know the process of
chi-sqared test,but how can I found it?
Which function are you using? Let's assume chisq.test. Then the first
thing to do is to type the function name:
chisq.test
Hello,
The follwing code generates 1000 dendograms from 1000 input binary matrices.
I dont't know how to generate a consensus dendogram from the 1000.
Can you help me ?
#code
library(ade4)
library(cluster)
library(stats)
for (i in 1:1000) {
#read each file
Thank you Andy. It worked.
Farrel Buchinsky, MD --- Mobile (412) 779-1073
Pediatric Otolaryngologist
Allegheny General Hospital
Pittsburgh, PA
-Original Message-
From: Liaw, Andy [EMAIL PROTECTED]
Sent: Monday, January 16, 2006 20:47
To: r-help@stat.math.ethz.ch; Farrel Buchinsky
Dear R useRs,
I have a vector with positive and negative numbers:
A=c(0,1,2,3,0,4,5)
Now if i-th element in vector A is 0, then i-th element in vector B
is a+1
else i-th element in vector b=a (or 0)
vector A: 0 1 2 3 0 4 5
vector B: 0 2 3 4 0 5 6
What's the right way to do this. I still
B - ifelse(A 0, A + 1, A)
?ifelse
Andrej Kastrin wrote:
Dear R useRs,
I have a vector with positive and negative numbers:
A=c(0,1,2,3,0,4,5)
Now if i-th element in vector A is 0, then i-th element in vector B
is a+1
else i-th element in vector b=a (or 0)
vector A: 0 1 2 3 0 4 5
one way is:
a - c(0,1,2,3,0,4,5)
b - ifelse(a 0, a + 1, a)
b
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: +32/(0)16/336899
Fax:
Andrej Kastrin wrote:
vector A: 0 1 2 3 0 4 5
vector B: 0 2 3 4 0 5 6
What's the right way to do this. I still have some problems with for and
if statements...
?ifelse perhaps...
A
[1] 0 1 2 3 0 4 5
B=ifelse(A0,A+1,0)
B
[1] 0 2 3 4 0 5 6
does a sort of element-wise
On 1/12/06, Seth Falcon [EMAIL PROTECTED] wrote:
On 12 Jan 2006, [EMAIL PROTECTED] wrote:
I'm using emacs-21.4 on debian unstable, together with the latest
ESS implementation. I try to change indentation to 4 by following
the advise in R-exts: It results in the following lines in my
On Mon, 16 Jan 2006, Dimitri Joe wrote:
(i) I get a big R file (for example, a 15Mb Stata file became a 42Mb R
file; after cleanup.import() from the Hmisc package, it drooped to 35Mb,
but that's still more than 2x the original Stata file) which, in turn, I
suspect is due the fact that
(ii)
Dear useRs!
I athought this was a trival question, however I could not fin dan answer in
the help files for print, format or formatC. I would like to print a
numerical matrix so that all cells (all rows/colums) are printed:
a) with the same number of decimal places (numbers after the decimal
If addition to the ifelse solution already posted one could do this
since a logical expression used in a numeric context is regarded
as 1 for TRUE and 0 for FALSE.
B - A + (A0)
On 1/17/06, Andrej Kastrin [EMAIL PROTECTED] wrote:
Dear R useRs,
I have a vector with positive and negative
On Tue, 17 Jan 2006, Daniel Farewell wrote:
I'd like to fit a GLM to each of a number of subsets of some data. The
`family' argument to `lmList' (in lme4) has given me cause for optimism,
but so far I've only been able to achieve linear model fits. For example
df - data.frame(gp = gp.temp
a-c(0,1,2,3,0,4,5)
b-vector(length=length(a))
b[a0]-a[a0]+1
b[a=0]-a[a=0]
b
[1] 0 2 3 4 0 5 6
On Tue, 2006-17-01 at 15:30 +0100, Andrej Kastrin wrote:
Dear R useRs,
I have a vector with positive and negative numbers:
A=c(0,1,2,3,0,4,5)
Now if i-th element in vector A is 0, then
To add onto an already clear explanation (a comment on precision in Stata).
Indeed Stata stores all numbers as floats
(also known as single precision or 4-byte reals). One way you could check this
is to save a small subset of your data
with all numbers as doubles in stata and see how that size
Hi Uwe,
I am aware of the legend option in xyplot, but I can't figure out how to
make it work. In particular, I do not understand how to use legend along
with the key argument. For my example, I would like to have a legend box
inside the plotting frame to indicate the 3 different smoothed
Many thanks! That'll work great. It's always good to discover a new, general,
function like by().
I would still be interested to know how the family argument to lmList() should
be used.
Daniel
Thomas Lumley [EMAIL PROTECTED] 01/17/06 3:15 pm
On Tue, 17 Jan 2006, Daniel Farewell wrote:
I'd
Here are a couple of alternatives to try:
noquote(format(round(M,3)))
noquote(apply(round(M,3), 2, format))
On 1/17/06, Aleš Žiberna [EMAIL PROTECTED] wrote:
Dear useRs!
I athought this was a trival question, however I could not fin dan answer in
the help files for print, format or
Hello,
probably you should have a look to the RandomFields package.
Regards. Olivier
--
Message: 59
Date: Tue, 17 Jan 2006 10:45:46 +0100
From: Eivind Sm?rgrav [EMAIL PROTECTED]
Subject: [R] Kriging for d3
To: r-help@stat.math.ethz.ch
Message-ID:
[EMAIL
Hello,
I want to test if the Weibull distribution is appropriate for the
failure time. When trying to reproduce an example from MASS (the book,
Ch. 13.2), I type
library(survival)
library(MASS)
leuk.wei - survreg( Surv(time)~ag+log(wbc),data=leuk)
ntimes -
On 1/17/06, Ravi Varadhan [EMAIL PROTECTED] wrote:
Hi Uwe,
I am aware of the legend option in xyplot, but I can't figure out how to
make it work. In particular, I do not understand how to use legend along
with the key argument. For my example, I would like to have a legend box
inside the
You could also write your numbers using SI suffixes.
Below is a function I use for this purpose.
The option near allows you to require that numbers between 0.001 and 1000 are
written as decimal; i.e. 0.010 appears as 0.010 instead of 10.0m.
##
## num2SI converts numbers to SI suffixed numbers
##
Thanks, Uwe
that clears up why I can't make R2WinBUGs work with OpenBUGS and WinBUGS1.5 :)
Both work pretty good with Wine in a GUI. I noticed that when I tried
rbugs, it does succeed in starting WinBUGS GUI, but then nothing
happens. I'll get WinBUGS1.4 and see what happens.
In the meanwhile,
Thank you!
That solves my a) problem!
Best,
Ales
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of Gabor Grothendieck
Sent: Tuesday, January 17, 2006 4:55 PM
To: Aleš Žiberna
Cc: r-help@stat.math.ethz.ch
Subject: Re: [R] Printing numerical matrices
Hi,
I'm looking for a more elegant (and faster) solution to my current
problem than the code at the end. I'm sure there is one, but can't think
where to look - any pointers would be very welcome. The problem is one
of resampling within an array. This array consists of 0s, 1s and NAs.
For each
On 1/17/06, Daniel Farewell [EMAIL PROTECTED] wrote:
Many thanks! That'll work great. It's always good to discover a new, general,
function like by().
I would still be interested to know how the family argument to lmList()
should be used.
As you have discovered the family argument has no
Dear mailing group,
I have loaded an Excel file into R by calling it .csv and using the
read.csv function in R. However then I want to use the (limma package
specific, I believe) function hclust, which clusters data in a tree
dendrogram, by similarity. However, I receive the errors msg.s: 1)
Hello,
The following code does'nt work for me. The last command reports an error. I
have created a consensus tree using the consensus comand from phylo but
cannot manipulate the phylo object afterwards to create a dendogram , by
transforming the phylo object into a hclust object and then into a
I am trying to use the friedman.test() on a data frame, d, but
I am receiving the following error message:
d
AW HS IAC WA
1 6 8 3 5
2 2 2 3 6
3 7 7 8 3
4 8 5 4 5
20 2 5 2 7
21 7 7 6 7
22 7 8 6 8
23 6 8 4 5
24 5 7 5 2
friedman.test(d)
Hi,
I am despairing of getting the indices for a vector:
First, I have a table from which I kick out a number of rows with na.omit.
Next, I use this table for clustering with kmeans and cl$cluster contains my
clusters. The cl$cluster is a vector which still contains my original indices
Dear R users,
There is a method called style analysis where you make a regression being
Y=fund yield and X=benchmarks yield, where we have the restrictions to
calculatethe linear regression:
1. The regression must don have the intercept term.
2. The coefficient sum must be one.
3. All
Hi All,
I am wondering if there is any literature or any prior implementations
of cluster analysis for only nominal (categorical) variables for a
large dataset, apprx 20,000 rows with 15 variables.
I came across one or two such implementations, but they seem to assume
certain data distributions.
Dear All:
I have a question on using coxph for multiple genes:
I have written code to loop through all 22283 genes in the Hgu-133A and
apply coxph on survival data.
However, I don't know how to work with the result for each gene:
Dear All:
I have a question on using coxph for multiple genes:
I have written code to loop through all 22283 genes in the Hgu-133A and
apply coxph on survival data.
However, I don't know how to work with the result for each gene:
I have been asked to analyse the results of (what is to me) a very
complicated experiment.
The dependent measure is the estimated distance, which is measured as a
function of the actual distance. There are also several other IVs.
The plot of log estimated distance as a function of log distance
On 1/16/06, Daniel A. Powers [EMAIL PROTECTED] wrote:
-- R-List
Can someone tell me how to get fitted values etc. after fitting lmer?
for example, from lme, I can fit mod.1 - lme() and get fitted values,
coefficients, etc. in this way
mod.1$fitted[,1] or mod.1$fitted[,2] etc.
It
Hao,
I'm not sure but you have specified your modell as tautology. The
formula below should be enough:
survtest - coxph(Surv(fup_interval, endpoint) ~ geneid, data =
pcc.primary.stg.3.cox)
HTH
Fredrik
- Original Message -
From: Hao Liu [EMAIL PROTECTED]
To: r-help@stat.math.ethz.ch
I have been using the ECDF function in the Hmisc package to produce
cumulative distribution function plots. The problem is that for small
datasets the steps look bad (not my characterization but from the
client). Is there a way to get the same information but smoothed? I have
tried the densityplot
Thanks for your reply, but I am not sure if I understand your reply right:
are you saying the geneid in your post is a vector of all 22283 genes? I
want to do cox fit to each and everyone of the 22283 genes and need a way to
summarize it -- to do that, I need to know how to parse the result for
Hi all,
I want to create an array of datetime.
If I have a datetime object dt
dt - strptime(10Jan2006 00:00:15, %d%b%Y %H:%M:%S)
dt
[1]2006-01-10 00:00:15
I want to make an array of dt, say 100 size. I got those error.
[1] 2006-01-10 00:00:15
dtarray-array(dt, dim=c(100));
Error in
On 1/17/06, Mike Bock [EMAIL PROTECTED] wrote:
I have been using the ECDF function in the Hmisc package to produce
cumulative distribution function plots. The problem is that for small
datasets the steps look bad (not my characterization but from the
client). Is there a way to get the same
The plotmeans() function uses the order of the factor levels. To change the
order in the plot, change the order of the factor levels.
For example:
data(state)
plotmeans(state.area ~ state.region)
Plots the groups in the order
levels(state.region)
levels(state.region)
[1] Northeast
Dear all,
I have run a hierarchical log-linear analysis using loglin {stats}
and came up with a specific model. loglin returns me the parameter
estimates giving me an idea in which direction the observed
frequencies deviate from the expected ones for my different factors.
To assess the
Thanks, I can do calculatation for each gene, however what I want to do is
to fit the model on each and every gene, store their result and then be able
to access all of them, then filter on them, to do that:
1. I need to know how to access each part of the coxph result, like its
coefficient,
I am using the intervals function with a gnls (example below). How can
I use some of the coefficients below in another calculation? Or asked
another way, how do I index these values?
Thanks.
Rick
intervals(yy$fit[[1]])
Approximate 95% confidence intervals
Coefficients:
Norman,
You're missing a step. You need to convert the data file into a
'dist' object, which is either a distance or dissimilarity matrix. This
is typically done by function dist(), but may also be done by other
functions which produce dist objects, like daisy() in package cluster,
Chang Shen Chang_Shen at progressive.com writes:
I am new to R. I try to search the web but could not find the answer so I
post it here asking for help.
Broad web searches may not be as helpful as following the advice in the posting
guide:
http://www.R-project.org/posting-guide.html
Hi
I have a linear model and I want to tests whether the model terms are
significant
I used anova but F and P value depend on the order of the terms in the
model.
I have repeated the same analysis in another stat software and F and P value
did not differ with order of terms in the model.
Can
I have two variables X and Y that have a bivariate standard normal
distribution. From the data I can obtain the probability that
P(X = a, Y= b) where a and b are cutoff values. Is there a way to
calculate the Corr(X,Y) using this information?
Thank you very much,
Marina
Hi:
I am using step.glm() as follows:
form1 - as.formula(haspdata ~ 1)
lg.mod1 - glm ( formula=form1, data=st.mtx,
family=binomial , na.action=na.omit )
upper - as.formula( haspdata ~ (
c5+childnm+educ1+incpov1+marital+msa+racekid+racemom+sex+shotcard )^2)
As you apparently haven't received any answer yet ...
Assuming your table is a data frame (str() will tell you), I believe you
are confusing the row names with indices. Please read the docs on data frame
for details (?data.frame). Row names are character strings that can be
indexed as such:
Douglas Bates [EMAIL PROTECTED] writes:
On 1/16/06, Daniel A. Powers [EMAIL PROTECTED] wrote:
-- R-List
Can someone tell me how to get fitted values etc. after fitting lmer?
for example, from lme, I can fit mod.1 - lme() and get fitted values,
coefficients, etc. in this way
Paul Johnson wrote:
Thanks, Uwe
that clears up why I can't make R2WinBUGs work with OpenBUGS and WinBUGS1.5 :)
Both work pretty good with Wine in a GUI. I noticed that when I tried
rbugs, it does succeed in starting WinBUGS GUI, but then nothing
happens. I'll get WinBUGS1.4 and see what
I suggest you consult a local statistician or read up on linear
models,perhaps by reading the relevant section of VR's MASS. Also search
the R-Help archives for Type III SS. Further comments below.
-- Bert Gunter
Genentech Non-Clinical Statistics
South San Francisco, CA
The business of the
Hi
I do not get why the symbols function produces warnings when axes=F is
added. The following example illustrate this
symbols(0,10, inches=T, circles=1, axes=F, xlab=, ylab=)
Warning message:
parameter axes could not be set in high-level plot() function
I augmented symbols and added the
On 17 Jan 2006 22:00:38 +0100, Peter Dalgaard [EMAIL PROTECTED] wrote:
Douglas Bates [EMAIL PROTECTED] writes:
On 1/16/06, Daniel A. Powers [EMAIL PROTECTED] wrote:
-- R-List
Can someone tell me how to get fitted values etc. after fitting lmer?
for example, from lme, I can fit
On Tue, 17 Jan 2006, Jean Eid wrote:
Hi
I do not get why the symbols function produces warnings when axes=F is
added. The following example illustrate this
symbols(0,10, inches=T, circles=1, axes=F, xlab=, ylab=)
Warning message:
parameter axes could not be set in high-level plot()
Mike Bock wrote:
I have been using the ECDF function in the Hmisc package to produce
cumulative distribution function plots. The problem is that for small
datasets the steps look bad (not my characterization but from the
client). Is there a way to get the same information but smoothed? I have
Hi, Doug:
I think it would help me to have an option like you just suggested --
with your comment below included in the help file for that option.
Thanks for all your hard work and creativity in this.
Best Wishes,
spencer graves
Douglas Bates wrote:
Dear Martin,
The anova() function computes sequential (sometimes called type-I) tests;
I'm not sure what Jmp does by default, but either type-II or type-III
tests would produce these results for a model of this structure. See the
Anova() function in the car package for type-II and III tests.
I
On 1/16/06, Antje Schüle [EMAIL PROTECTED] wrote:
Dear group,
First I provide you with an example, I found in the newsgroup. Then I'd like
to explain my problem to you by means of the output.
enviro -
data.frame(Year = rep(2001:2002, each = 365),
Day = rep(1:365, 2),
Dear sir or ma'am,
I have difficulty in installing packages. When I click Install package(s)...
in the Packages menu, it takes forever and it finally shows the following
message:
--- Please select a CRAN mirror for use in this session ---
Warning: unable to access index for repository
Please use POSIXct and not POSIXlt objects, which are lists.
On Tue, 17 Jan 2006, Chang Shen wrote:
Hi all,
I want to create an array of datetime.
If I have a datetime object dt
dt - strptime(10Jan2006 00:00:15, %d%b%Y %H:%M:%S)
dt
[1]2006-01-10 00:00:15
I want to make an array of dt,
On Tue, 17 Jan 2006, Thomas Lumley wrote:
On Tue, 17 Jan 2006, Jean Eid wrote:
Hi
I do not get why the symbols function produces warnings when axes=F is
added. The following example illustrate this
symbols(0,10, inches=T, circles=1, axes=F, xlab=, ylab=)
Warning message:
parameter axes
Hello!
Re R packages:
- R2WinBUGS is compatible with WinBUGS-1.4.x only, its newest version
can speak with WinBUGS under wine thanks to user contributions. But it
still depends on WinBUGS-1.4.x, hence Windows only (considering wine as
Windows).
However Andrew Gelman, has added also
You seem to be using Windows. I do suggest you ask local advice, as the
problem is almost certainly local your site.
You may need to set a proxy -- see the rw-FAQ.
You may need to set a policy to allow R through a firewall.
You may need permission to download at all, and indeed the trouble you
Use multinom or glm instead. The IFP algorithm used in loglin is not
designed to find parameter estimates let alone standard errors.
Fitting log-linear models is discussed in all good books on R/S, e.g.
MASS (see the FAQ for full details).
On Tue, 17 Jan 2006, Becker Cordula wrote:
Dear
There is no step.glm in R!
I suggest that you forget about step.glm, and start afresh with step in R
(which is a cut-down version of MASS's stepAIC rather than emulating
step.glm in S). In particular, it avoids many of the pitfalls of step.glm.
Please do study the posting guide, and tell us
?friedman.test
Says:
Description:
Performs a Friedman rank sum test with unreplicated blocked data.
Usage:
friedman.test(y, ...)
y: either a numeric vector of data values, or a data matrix.
So assuming your data, d, is unreplicated blocked data, perhaps
d=as.matrix(d)
friedman.test(d)
or
I mentioned this on the ESS list a little while ago, I made an installer
for ESS and XEmacs on Windows.
It worked for me but given my minimal knowledege of ESS and XEmacs, it
might not be the right way to do it and may or may not work for you.
I hosted it and the source for the inno setup
I suggest you not worry about the loops. A decade ago (e.g., with
S-Plus 3 or 3.1), loops were a major problem. Releases of S-Plus and R
since then have made substantial improvements in loop computations.
My preferred solution to your problem, as I understand it, is as
Thank you Gabor, the 'aggregate' function
operating on zoo(WS) does the job beautifully!
'zoo' has a lot of other goodies too, great
package (thanks for that too).
Cheers,
Augusto
Augusto Sanabria. MSc, PhD.
Mathematical Modeller
Risk Research
When I try to install RMySQL from source on Mac OS X 10.4.4 (iMac G5
2.1 Ghz), I get this...
[tcom122-223-dhcp:2.1.1/Resources/library] root# setenv PKG_CPPFLAGS
/usr/local/mysql/lib
[tcom122-223-dhcp:2.1.1/Resources/library] root# setenv PKG_LIBS
/usr/local/mysql/include
Thanks! Let me ask this question again. Clearly, if you experts
can't make R talk to WinBUGS, then I can't either. So, If bugs()
doesn't work, What is the next best thing?
With wine, I can run OpenBUGS and WinBUGS, but I cannot send jobs from
R to a BUGS program (still trying, some people say
Hello,
I am new to using R and I am having problems get boot() to work properly.
Here is what I am trying to do:
I have statistic called cs. cs takes a data matrix (154 x 5) and calculates
12 different scores for me. cs outputs the data as a vector (12 x 1). cs
doesn't really use
The first thing you are doing wrong is that you are not including a
copy of cs for us to see ;).
Based on what you have written, I speculate that cs does not use the
index correctly. if so then a simple, although inefficient,
workaround is to rewrite cs:
cs - function(dataframe, index) {
Before I forget, found the working recipe for rbugs. My mistake
before was not
realzing that the n.iter value is total iterations, including burnin,
and so by setting
n.iter=1000 and n.burnin=1000, I was leaving 0 iterations for the updates.
### Paul Johnson 2006-01-18. This does work!
### Works
Paul Johnson wrote:
Do you mean to say that you have actually made OpenBUGS run with
R2WinBUGS in Linux?
No, I did not say this.
Gelman's page seems to state that OpenBUGS support is brought in from
BRugs, which is still Windows-only.
Well, Gelman changed his site a bit. Few days (weeks?)
Hi, R users:
I have a data.frame (not a matrix), I got a vector with the same length
as the
number of records (rows) of the data frame, and each element of
that vector is the column number (in a specific range of columns) of the
corresponding
record that I must set to zero.
How can I do
93 matches
Mail list logo