I have a fix in R2WinBUGS_2.1-15 which is on its way to CRAN and will
appear there within 2 days, probably.
Best wishes,
Uwe Ligges
Luwis Tapiwa Diya wrote:
Dear R users
I have been encountering an error (in the Winbugs window: Blackbox - index
out of range . ). So to me it seems
I have a fix in R2WinBUGS_2.1-15 which is on its way to CRAN and will
appear there within 2 days, probably.
Best wishes,
Uwe Ligges
Paul Heinrich Dietrich wrote:
This weekend I noticed that my R2WinBUGS connection was no longer working on
my Windows computer at work AND my Ubuntu linux
See the manual Writing R Extensions, Section 6.5.1: Printing from FORTRAN
Best wishes,
uwe Ligges
rkevinbur...@charter.net wrote:
Hello,
From the date time stamps I see that not only do I have the results of the
build in 'gnuwin32/front-ends' but I can also see the time stamps have
Kim Vanselow wrote:
Dear r-helpers!
Jan Thiele asked r-help-list in february, but I could not find an answer to his
question. I have the same problem:
While preparing figures of 'envfit' plots with vegan for publication, I ran
into a layout problem, that I found no solution for in the
What is the problem that you see?
(I think the trouble is that Laplace density is not smooth at x=0,
while a kernel density estimate returns a smooth density)
It would help to know what you are trying to achieve with this
exercise...
Giovanni
Date: Thu, 05 Nov 2009 14:36:44 +0100
From:
Hi Guys,
Of course, a backward, forward, or other methods can be used directly. But
concerning BMA, the model interpretation is far simple:
Bayesian Model Averaging accounts for the model uncertainty inherent in the
variable selection problem by averaging over the best models in the model
There is also a sparse PLS model in the spls package. It uses
lasso-like regularization to reduce the number of variables. I've had
a lot of success with it.
Max
2009/11/5 Ricardo Gonçalves Silva ricard...@terra.com.br:
Hi Guys,
Of course, a backward, forward, or other methods can be used
Dear All,
I'm trying to use the ols function in the Design library (version
2.1.1) of R to estimate parameters of a linear model, and then use the
contrast function in the same library to test various contrasts.
As a simple example, suppose I have three factors: feature (3
levels), group
Timothy Clough wrote:
Dear All,
I'm trying to use the ols function in the Design library (version
2.1.1) of R to estimate parameters of a linear model, and then use the
contrast function in the same library to test various contrasts.
As a simple example, suppose I have three factors:
try this:
x - read.table(textConnection(Year SpA
+ 2000 0
+ 2000 2
+ 2000 1
+ 2001 8
+ 2001 2
+ 2001 0
+ 2001 0
+ 2002 1
+ 2002 2), header=TRUE)
closeAllConnections()
# convoluted if you want a 'function' for each year
l.func - list('2000'=function(x) x / 146,
+
Dear R-users,
I've just installed the last version of R (2.10.0) in Windows and found a
completely new version of the help pages.
Basically, they now appear as a simple (even if slightly embellished) text file.
I tried to install it with a customized startup and realized that the old CHM
files
I took the subset of the data you send, did a 'write.csv' and then
imported it to Excel and the columns seem fine (see the .png file
attached).
So exactly what problems are you having?
On Thu, Nov 5, 2009 at 12:49 PM, Hayes, Daniel
d.j.ha...@liverpool.ac.uk wrote:
No luck with 'write.csv'
I need to use the Bhattacharyya distance metric to determine population
separation. Has anyone written a Bhattacharyya distance metric function in
R?
--
View this message in context:
http://old.nabble.com/Bhattacharyya-distance-metric-tp26221259p26221259.html
Sent from the R help mailing list
See the NEWS and CHANGES files. You can choose html help, if you prefer
that, or pdf. Support for compiled html help has been withdrawn.
Uwe Ligges
antonio.gasparr...@lshtm.ac.uk wrote:
Dear R-users,
I've just installed the last version of R (2.10.0) in Windows and found a completely new
Dear R,
Why rm(list-ls()) gives an error but rm(list=ls()) not? I remember the
operator â-â can be used anywhere...
Thanks!
Feng
--
Feng Li
Department of Statistics
Stockholm University
106 91 Stockholm, Sweden
http://feng.li/
[[alternative HTML version deleted]]
Antonio,
Starting from 2.10, R dropped its support for Windows CHM help. It has
been replaced with the web browser HTML based one. To activate this add
the following line to your Rprofile
options(help_type=html)
The pages come up in your browser. They all have the same links as
before.
- and = are not universally interchangable.
args(rm)
function (..., list = character(0L), pos = -1, envir = as.environment(pos),
inherits = FALSE)
The call
rm(list - ls())
assigns the result of ls() to the variable 'list' and passes that value as an
anonymous argument to rm()
Feng Li wrote:
Dear R,
Why rm(list-ls()) gives an error but rm(list=ls()) not? I remember the
operator ‘-’ can be used anywhere...
Yes, and it means that you make an assignment once passed to the first
argument ... in rm() and evaluated. Well, it is just never evaluated
since ...
I haven't tried this yet, but I should note that my Site Search page:
http://finzi.psych.upenn.edu,
and the associated R function RSiteSearch and related functions,
depend on compiled html help files.
The idea is to be able to search the help page of functions that you
don't have installed yet,
Jonathan Baron wrote:
I haven't tried this yet, but I should note that my Site Search page:
http://finzi.psych.upenn.edu,
and the associated R function RSiteSearch and related functions,
depend on compiled html help files.
On regular html, not compiled html, I believe.
The idea is to be
I have a column of dates in this format.. 19900501
I want to change
19900501 into 1990 05 01
then append a new column on the end (right hand side of spreadsheet)
:confused:
--
View this message in context:
http://old.nabble.com/19900501-into-1990-05-01-tp26217515p26217515.html
Sent from the
Hello,
Thank you very much for looking at this. I have a seasonal user for R. I
teach my undergrads and graduates students statistics using R and often find
myself trying to solve problems to process student collected data in an
efficient way.
In this case, I have a data.frame with multiple
Hello, I'm a totally newbie to R and I'm taking a class using S+.
In the class we use the multcomp command which takes a aov object and
calculates confidence intervals for all pairwise differences by the Fisher
least significant differences method.
How can I do this in R.
Thank you for taking
Hello all,
Thank you for you help with the first part, it worked perfectly. I now have
a dataset that is culled to only include rows with no NA's or only 1 NA.
What I want to achieve with the graph is:
Plot the new data matrix:
Alpha Beta Gamma Delta
A 12 34
B NA
Dear R,
Why rm(list-ls()) gives an error but rm(list=ls()) not?. I remember the
operator â-â can be used anywhere...
Thanks!
Feng
--
Feng Li
Department of Statistics
Stockholm University
106 91 Stockholm, Sweden
http://feng.li/
[[alternative HTML version deleted]]
Thank you all for your help.
None of the suggestions worked but this is because I believe I was opening the
created file into Excel the wrong way.
I managed to open the file in Word and then do indeed get a cvs format. This I
was then able to transform to a normal table.
Will now try to install
No luck with 'write.csv' either.
Thanks for helping though.
Cheers,
Daniel
-Original Message-
From: jim holtman [mailto:jholt...@gmail.com]
Sent: 05 November 2009 14:59
To: Hayes, Daniel
Cc: r-h...@lists.r-project.org
Subject: Re: [R] collumn error when exporting to Excel
try using
I am missing something fundamental. I ran the function nlm, but I don't
understand how to extract the optimal solution as a numeric vector. The
function produces it as one element of a list. I don't see anything in
the R documentation about converting such a list element to the vector it
Thanks. I misunderstood compiled. Sorry.
I need static html.
I cannot install all the packages (for several reasons). That is why
I wrote my script to extract the html pages only. The parts of R that
it uses are:
/usr/share/R/perl/build-help.pl
tools:::.writePkgIndices()
make.packages.html()
Hi:
Suppose that I have a data frame as below
x1 x2 x3 ... x10 wk1 wk2 ... Wk208 (these are the column names)
For each record, x1, x2, x3 ... x10 are attributes. and wk1, wk2, ..., wk208
are the sales recoreded for this attribute combination. Suppose that now, that
I want to do the following
Tena koe Satish
I'm not entirely sure what you want, but did you check aggregate()?
HTH
Peter Alspach
-Original Message-
From: r-help-boun...@r-project.org
[mailto:r-help-boun...@r-project.org] On Behalf Of Vadlamani,
Satish {FLNA}
Sent: Friday, 6 November 2009 10:16 a.m.
qwe84 wrote:
Hello, I'm a totally newbie to R and I'm taking a class using S+.
In the class we use the multcomp command which takes a aov object and
calculates confidence intervals for all pairwise differences by the Fisher
least significant differences method.
How can I do this in R.
Thank
Try something like this:
# get column numbers of 'wk' columns
col.num - grep(^wk, names(yourDF))
# split out only the 'wk' columns
result - lapply(split(yourDF[, col.num], paste(yourDF$x2,
paste(yourDF$x3))), function(.mrg){
colSums(.mrg)
})
On Thu, Nov 5, 2009 at 4:15 PM, Vadlamani, Satish
Something like:
format(as.Date(20090501,format=%Y%m%d), %Y %m %d)
r
-
Remko Duursma
Post-Doctoral Fellow
Centre for Plants and the Environment
University of Western Sydney
Hawkesbury Campus
Richmond NSW 2753
Dept of Biological Science
Macquarie
Join your colleagues for the 20th highly acclaimed international
forum to hear the latest implementations of biometrics...
The Winter 2010 BIOMETRICS SUMMIT:
Practical Implementation Strategies, Market Trends And Best Practices
In Government And Business
No luck with:
(col.names = NA and row.names = TRUE
Will check the other responses and else try to install writeXLS
Thanks for the help.
Daniel
-Original Message-
From: Leandro Marino [mailto:lean...@cesgranrio.org.br]
Sent: 05 November 2009 15:52
To: Hayes, Daniel;
Dear All,
I am trying to simulate data for a spline/piecewise regression model. I am
missing something fundamental in my simulation procedure because when I try to
fit my simulated data using the Gauss-Newton method in SAS, I am getting some
wacky parameter estimates. Can anyone please check
Gene,
You might want to look at function read.csv.ffdf from package ff which can read
large csv-files into a ffdf object. That's kind of data.frame which is stored
on disk resp. in the file-system-cache. Once you subscript part of it, you get
a regular data.frame.
Jens Oehlschlägel
--
If the headers all start with the same letter, A say, and the data
only contain numbers on their lines then just use
read.table(..., comment = A)
On Mon, Nov 2, 2009 at 2:03 PM, Gene Leynes gleyne...@gmail.com wrote:
I've been trying to figure out how to read in a large file for a few days
Hi Bob,
On Nov 5, 2009, at 3:04 PM, robertag...@discover.com robertag...@discover.com
wrote:
I am missing something fundamental. I ran the function nlm, but I
don't
understand how to extract the optimal solution as a numeric vector.
The
function produces it as one element of a list. I
Hi:
If I want to read a file with read.table. I want x1 and x2 to be read as
character and x3 as numeric. How to do this? Thanks.
Satish
x1 ,x2,x3
10,20,30
11 ,22,35
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
Thanks to everyone who replied.
I think the html based one works fine (even if I definitely preferred the old
chtml).
Only another comment: I chose html in my installation with Windows Vista, but
still the help pages are opened in text files in the newly installed version.
Of course, I can
See
?read.csv
and read the part about colClasses.
X - read.csv(textConnection(x1 ,x2,x3
+ 10,20,30
+ 11 ,22,35), colClasses = c(character, character, numeric)
+ )
str(X)
'data.frame': 2 obs. of 3 variables:
$ x1: chr 10 11
$ x2: chr 20 22
$ x3: num 30 35
-Ista
On Thu, Nov 5,
On 11/06/2009 02:24 AM, Zhijiang Wang wrote:
Dear All,
I have modified a .R code (/usr/local/lib/R/site-library/brainwaver/R), but
when I run it, why does it run the original code? Whatever I modified it,
the results did not change.
what happened?
Hi Zhijiang,
If you
I'm wondering if there is any option available in load() such that I
can specify which variable I want to load from an RData file. I don't
see such option in the help.
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
On 11/06/2009 03:21 AM, Joel Fürstenberg-Hägg wrote:
Hi all,
I'm making a PCA plot with eight variables (columns taken from a larger data frame
fieldTrial0809[idx,c(39:46)]). I then want the symbols in the plot to be colored as a
gradient from red to blue, depending on the value of another
Create an environment, load the RData file into the environment, copy
over the items of interest, delete the environment.
e1 - new.env()
load(gold.RData, e1)
ls()
[1] e1
ls(e1)
[1] a.mark aisleRange allSlots
balanced
[5] baseSlotsbox.offset
If you just want it to remain as character data:
sub(()(..)(..), \\1 \\2 \\3, 19000501)
[1] 1900 05 01
On Thu, Nov 5, 2009 at 2:04 PM, frenchcr frenc...@btinternet.com wrote:
I have a column of dates in this format.. 19900501
I want to change
19900501 into 1990 05 01
then append a
Hi folks, sorry for this beginner question but what means a p-value = NA on
a menn-whitney test?
v1 - c(0.022, 0.022, 0.022, 0.022, 0.022, 0.022)
v2 - c(0.022, 0.022, 0.022, 0.022, 0.022)
wilcox.test(v1, v2, PAIRED=FALSE)
W = 15, p-value = NA
I know that there's no statistical difference
Hi Ambar,
v1 and v2 are constants. they need to be variables.
try e.g.,
v1 - c(0.021, 0.022, 0.022, 0.022, 0.022, 0.022)
v2 - c(0.021, 0.022, 0.022, 0.022, 0.022)
wilcox.test(v1, v2, PAIRED=FALSE)
-Ista
On Thu, Nov 5, 2009 at 7:23 PM, Ambar Amarelo ambar.amar...@gmail.com wrote:
Hi folks,
Ambar Amarelo wrote:
Hi folks, sorry for this beginner question but what means a p-value = NA on
a menn-whitney test?
v1 - c(0.022, 0.022, 0.022, 0.022, 0.022, 0.022)
v2 - c(0.022, 0.022, 0.022, 0.022, 0.022)
wilcox.test(v1, v2, PAIRED=FALSE)
W = 15, p-value = NA
I know that there's no
Thank you. When the RData file is huge, it will take a long time to
load it even though I may only need a small variable. Could somebody
add an option to 'load()' to load only a few given variables?
On Thu, Nov 5, 2009 at 5:58 PM, jim holtman jholt...@gmail.com wrote:
Create an environment, load
Peng Yu wrote:
Thank you. When the RData file is huge, it will take a long time to
load it even though I may only need a small variable. Could somebody
add an option to 'load()' to load only a few given variables?
Check out the filehash package. It has an option that allows you to
Hi,
I am create a new DF that summarizes the mean angle per treatment/trial, of
the original DF (see below). I have had some success using:
(tapply(df$Angle, INDEX=interaction(df$State,
df$Trial), FUN=mean)); however, this gives the answer as a list, which
means I would have to split the name
I think the html based one works fine (even if I definitely preferred the
old chtml).
Only another comment: I chose html in my installation with Windows Vista,
but still the help pages are opened in text files in the newly installed
version.
The nice thing about the chtml that is missing from
Hi R users
I need the map of France’s « communes » (towns) to build a map
Is there a way to get it?
More generally:
How to do to get the map of a country and its different geographical levels?
Best regards
--
View this message in context:
On Thu, 2009-11-05 at 12:25 -0800, Diana Garrett wrote:
I need to use the Bhattacharyya distance metric to determine population
separation. Has anyone written a Bhattacharyya distance metric function in
R?
Something like this ( I think it's correct ):
--%---
DBhat - function(X1,X2){
# define
I recently started trying R and SSOAP and was able to successfully try a hello
world service. I am now trying to get a more complicated interface to work with
SSOAP and so far failed miserably at that and so need any help I can
get from here.
The service I am attaching is a prototype for a full
Tena koe Jim
?aggregate
HTH
Peter Alspach
-Original Message-
From: r-help-boun...@r-project.org
[mailto:r-help-boun...@r-project.org] On Behalf Of jimdare
Sent: Friday, 6 November 2009 2:54 p.m.
To: r-help@r-project.org
Subject: [R] Calculate Mean for each Treatment/Trial
On Thu, Nov 5, 2009 at 5:58 PM, jim holtman jholt...@gmail.com wrote:
Create an environment, load the RData file into the environment, copy
over the items of interest, delete the environment.
e1 - new.env()
load(gold.RData, e1)
ls()
[1] e1
ls(e1)
[1] a.mark aisleRange
Hi Koraleus,
Here is a suggestion to (almost) do what you want:
# data set
x - read.table(textConnection(Alpha Beta Gamma Delta
A 12 34
B NA 2 45
D 8 9 10 11
E 5 NA 7 13), header=TRUE)
closeAllConnections()
x
# plot
On Nov 5, 2009, at 7:53 PM, jimdare wrote:
I am create a new DF that summarizes the mean angle per treatment/
trial, of
the original DF (see below). I have had some success using:
(tapply(df$Angle, INDEX=interaction(df$State,
df$Trial), FUN=mean));
That's rather difficult to accept,
Peng Yu wrote:
How to remove the environment 'e1' after I get 'nfree'?
Did you try rm( e1 ) ?
-
Charlie Sharpsteen
Undergraduate
Environmental Resources Engineering
Humboldt State University
--
View this message in context:
Hello everyone,
I have tried to look for this everywhere and so far have no luck. I have a
plotMeans(DV, IV1, IV2, error.bars=se) graph that plots my data
(DV-continuous, IVs are factors, IV1 - two levels, IV2-four levels). I am
trying to increase a scale of my y-axis (to be consistent with my
Hi there,
data1-matrix(data=c(1,1.2,1.3,3/23/2004,1,1.5,2.3,3/22/2004,2,0.2,3.3,4/23/2004,3,1.5,1.3,5/22/2004),nrow=4,ncol=4,byrow=TRUE)
data1-data.frame(data1)
names(data1)-c(areaid,x,y,date)
data1
areaid x y date
1 1 1.2 1.3 3/23/2004
2 1 1.5 2.3 3/22/2004
3 2 0.2 3.3
Hi Sergey
I've attached a script with a ylim argument added to plotMeans.
If you are using Rcmdr, you can load this script via the Rcmdr File menu item
File
Open Script File
so save this plotMeans.R script somewhere,
and load it and run it in Rcmdr, then
do your plot command
plotMeans(DV,
You have to use R CMD REMOVE Package NAME to uninstall the package(from
Command line) .
Then R CMD build brainwaver.tar.gz
the R CMD INSTALL brainwaver.tar.gz to re install then the changes you made
you be effective
Zhijiang Wang wrote:
Dear All,
I have modified a .R code
I download a csv extract from a database and use read.csv to read it
from R and when there are large character fields with embedded blanks,
slashes etc. - R often sees one line as two lines (or more).
I verfied with readLines that an embedded blank in a character field
causes a spurious new line
I had to deal with this issue yesterday, and I'm guessing many who use lme4
have this same issue because of the lack of a predict.lmer function. It
seems to me that one would need to plug the X %*% b output from the
predict.lmerBin function below into the logistic distribution to get the
actual
hi all -
is there a package or library that contains a function for partitioning the
chi-square statistic of an I X J contingency table into its respective
independent parts?
i looked around for this, but i didn't find anything. perhaps there's
another name for this sort of analysis? i know it
Thanks, Steve. I believe that's all I needed. I couldn't find that
r$estimate syntax anywhere in the manual. I'm new to R, having used SAS
exclusively in the past. I was able to run the optimizing functions nlm
and optim successfully, but I couldn't figure out how to access the
estimates.
thank all friends to discuss this problem, my data is 54*160 matrix. PLS is
a good method, can it give a equation with y~selected little vriables? for
exmaple:
y Sv Sp Ms nCIRnAB nC nN nO nX
ZM1V
1 7.6231.45 33.44 2.372
I would appreciate your insights on the attached session.
Bob
Robert Agnew | Discover
Director Acquisition Analytics
Marketing – Analysis Pricing
2500 Lake Cook Road
Riverwoods, IL 60015
Tel 224-405-1425 Fax 224-405-4971
robertag...@discover.com
I'm using the plotmeans function to plot the means and 95% CI of my groups,
however I also want to get the actual numbers of the CI, i.e. the upper and
lower bounds. Is there anyway to get that information out of plotmeans? If
not, is there an easy way to calculate the CIs?
--
View this
Hello everyone,
I have tried to look for this everywhere and so far have no luck. I have a
plotMeans(DV, IV1, IV2, error.bars=se) graph that plots my data
(DV-continuous, IVs are factors, IV1 - two levels, IV2-four levels). I am
trying to increase a scale of my y-axis (to be consistent with my
For a code .R file, I add a function in it.
But I saw it was not modified in command window, while it was changed in
/usr/local/lib/R/site-library/brainwaver/R by openning it by editor.
So I doubt the two are different, but I just find one version this .R file,
which was just in
Hi,
I am trying to a strucutural change analysis on a certain data set. Attached
here, with variable name xcd
http://old.nabble.com/file/p26226190/xcd.rda xcd.rda
Am using the following command:
bp.inrz-breakpoints(INR~SPX+WPI,data=xcd,h=26)
But keep on getting this error whatever variables
You have to use R CMD REMOVE Package NAME to uninstall the package(from
Command line) .
Then R CMD build brainwaver.tar.gz
the R CMD INSTALL brainwaver.tar.gz to re install then the changes you made
you be effective
Zhijiang Wang wrote:
Dear All,
I have modified a .R code
Bert,
Thanks for the tip. Yes, strsplit works, and works fast! For me,
white-space tokenization means splitting at the white spaces, so the
^ and the outermost square brackets should/can be omitted.
Regards ... from Basel to South San Francisco,
Richard
On Nov 3, 2009, at 22:03 , Bert
Hello,
Can anyone recommend a good example of web implementation of R? Can't seem
to find anything on my own.
Thanks in advance!
Sam
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
On Thu, Nov 5, 2009 at 11:03 PM, antonio.gasparr...@lshtm.ac.uk wrote:
Thanks to everyone who replied.
I think the html based one works fine (even if I definitely preferred the old
chtml).
Only another comment: I chose html in my installation with Windows Vista, but
still the help pages
On Thu, Nov 5, 2009 at 10:54 PM, Sam Albers tonightstheni...@gmail.com wrote:
Hello,
Can anyone recommend a good example of web implementation of R? Can't seem
to find anything on my own.
Thanks in advance!
Sam
What do you mean by web implementation?
An interactive session?
A
82 matches
Mail list logo