On Sun, 25 Jan 2009, Patrick Burns wrote:
Michael Kubovy wrote:
Dear R-helpers,
[...]
(2) If I remember dnorm() and want to be reminded of the call, I also get a
list of pages.
It sounds to me like here you want:
args(dnorm)
or, for functions hidden in a namespace, argsAnywhere().
On Sun, 25 Jan 2009, Jonathan Greenberg wrote:
Ah, perfect -- so would the "ideal" R_LIBS_USER setting (to more or less
guarantee the libraries will work on every possible computer) be something
along the lines of:
~/myRlibraries/%V%p%o%a
Or is this overkill?
%V is overkill. On some OSes %
On Sun, 25 Jan 2009, Karun Gahlawat wrote:
Prof Ripely,
Thanks. SUN gnu-iconv package should overwrite the sun version and so
it does. Apparently, it does not work.
As I said, you also need to pick up the correct header.
I built this library from gnu source with gcc and now it configures
and
Hi Jim,
I run the following code
*ds <- read.csv(file="D:/Shreyasee laptop data/ASC Dataset/Subset of the ASC
Dataset.csv", header=TRUE)
> attach(ds)
> str(dos)*
I am getting the following message:
*Factor w/ 12 levels "-00-00","6-Aug",..: 6 6 6 6 6 6 6 6 6 6 ...*
Thanks,
Shreyasee
On
do:
str(dos)
str(patientinformation1)
They must be the same length for the command to work: must be a one to
one match of the data.
On Sun, Jan 25, 2009 at 10:23 PM, Shreyasee wrote:
> Hi Jim,
>
> I tried the code which u provided.
> In place of "dos" in command "pat1 <- rbinom(length(dos), 1,
Probably not intentional, but there doesn't appear to be a link to R or any R
related material on the site. Ha. Found that interesting. Still a good
list...
--- On Sat, 1/24/09, David C. Howell wrote:
From: David C. Howell
Subject: Re: [R] Stat textbook recommendation
To: r-help@r-project.
Hi: if i understand, i think
newx<-x[ x %in% key]
should give you what you want.
On Sun, Jan 25, 2009 at 10:27 PM, Akshaya Jha wrote:
Hi,
I have the following datasets:
x=data I am looking through
key=a set of data with the codes I want
I have the following issue:
I want the subset of x
Hi Akshaka,
Take a look at ?"%in%". Here is an example for the help:
> x<-1:10
> key<-c(1,3,5,9)
> x %in% key
[1] TRUE FALSE TRUE FALSE TRUE FALSE FALSE FALSE TRUE FALSE
> x[ x %in% key ]
[1] 1 3 5 9
HTH,
Jorge
On Sun, Jan 25, 2009 at 10:27 PM, Akshaya Jha wrote:
> Hi,
>
> I have the f
Hi,
I have the following datasets:
x=data I am looking through
key=a set of data with the codes I want
I have the following issue:
I want the subset of x which has a code contained in the key dataset. That is,
if x[i] is contained in the key dataset, I want to keep it. Note that x may
contain
Hi Jim,
I tried the code which u provided.
In place of "dos" in command "pat1 <- rbinom(length(dos), 1, .5) # generate
some data"
I added "patientinformation1" variable and then I gave the command for
"tapply" but its giving me the following error:
*Error in tapply(pat1, format(dos, "%Y%m"), fun
YOu can save the output of the tapply and then replicate it for each
of the variables. The data can be used to plot the graphs.
On Sun, Jan 25, 2009 at 9:38 PM, Shreyasee wrote:
> Hi Jim,
>
> I need to calculate the missing values in variable "patientinformation1" for
> the period of May 2006 to
Thanks John--just needed to rtfm a little farther down :)
Anthony
John Fox wrote:
Dear Anthony,
From ?sem:
"If given as NA, the program will compute a start value, by a slight
modification of the method described by McDonald and Hartmann (1992). Note:
In some circumstances, some start values
Hi Jim,
I need to calculate the missing values in variable "patientinformation1" for
the period of May 2006 to March 2007 and then plot the graph of the
percentage of the missing values over these months.
This has to be done for each variable.
The code which you have provided, calculates the missi
Dear Anthony,
>From ?sem:
"If given as NA, the program will compute a start value, by a slight
modification of the method described by McDonald and Hartmann (1992). Note:
In some circumstances, some start values are selected randomly; this might
produce small differences in the parameter estimate
Here is an example of how you might approach it:
> dos <- seq(as.Date('2006-05-01'), as.Date('2007-03-31'), by='1 day')
> pat1 <- rbinom(length(dos), 1, .5) # generate some data
> # partition by month and then list out the number of zero values (missing)
> tapply(pat1, format(dos, "%Y%m"), functi
Hello-
If I input a variance-covariance matrix and specify NA for start values,
how does sem determine the start value? Is there a default?
Anthony
--
Anthony Steven Dick, Ph.D.
Post-Doctoral Fellow
Human Neuroscience Laboratory
Department of Neurology
The University of Chicago
5841 S. Maryla
Dear helpers,
As the title says, my question is not directly related to R.
I find, however, that there are many people who are both knowledgeable and
kind in this email list, and so decided to give it a try.
I do stochastic simulations. Parameter values used in simulation often come
from the obs
Hi Jim,
The dataset has 4 variables (dos, patientinformation1, patientinformation2,
patientinformation3).
In dos variable ther are months (May 2006 to March 2007) when the surgeries
were formed.
I need to calculate the percentage of missing values for each variable
(patientinformation1, patientinf
What does you data look like? You could use 'split' and then examine
the data in each range to count the number missing. Would have to
have some actual data to suggest a solution.
On Sun, Jan 25, 2009 at 8:30 PM, Shreyasee wrote:
> Hi,
>
> I have imported one dataset in R.
> I want to calculate
Hi,
I have imported one dataset in R.
I want to calculate the percentage of missing values for each month (May
2006 to March 2007) for each variable.
Just to begin with I tried the following code :
*for(i in 1:length(dos))
for(j in 1:length(patientinformation1)
if(dos[i]=="May-06" && patientinfor
sorry,
there is a comma missing;
fish[fish$GeoArea == 1 & fish$Month == 10, ]
Am 25.01.2009 um 23:33 schrieb Jörg Groß:
fish[fish$GeoArea == 1 & fish$Month == 10]
Am 25.01.2009 um 23:06 schrieb pfc_ivan:
I am a beginner using this R software and have a quick question.
I added a file i
Hello!
I am VERY new to genetic optimization and have a question about
rgenoud package (http://sekhon.berkeley.edu/rgenoud/):
Is it possible to modify the settings of rgenoud so that one could see
not only the "winner" solution but also the "runner ups"? By runner
ups I mean at least several othe
It is not that you are out of memory; one of your two 's2yg' objects
is not large enough (improperly dimensioned?) so you get the subscript
error. You can put the following in your script to catch the error
and then examine the values:
options(error=utils::recover)
On Sun, Jan 25, 2009 at 6:25
It looks in data and if not found there in environment(formula)
so try this:
mylm <- function(model, wghts) {
lm(model, data.frame(wghts), weights = wghts)
}
On Sun, Jan 25, 2009 at 4:20 PM, Wacek Kusnierczyk
wrote:
> dear list,
>
> below is an edited version of my response to an r user
I am writing a Gibbs sampler. I think it is outputting some of what I want,
in that I am getting vector of several thousand values (but not 10,000) in a
txt file at the end.
My question is, is the error message (see below) telling me that it can't
output 10,000 values (draws) because of a limitat
Hi Jörg,
If I understood,
apply(yourdataframe,2,function(x) x[diff(which(x==1 | x==3))])
should do what you want.
HTH,
Jorge
On Sun, Jan 25, 2009 at 5:53 PM, Jörg Groß wrote:
> is there a way to do that without generating a data.frame?
>
> In my real data, I have a big data.frame and I have
How about:
a <- c(1,2,3,3,2,1,6,3,2)
b <- c(NA,a[-length(a)])
c <- c(a[-1],NA)
a[b==1 & c==3]
[1] 2 6
hth,
Kingsford Jones
On Sun, Jan 25, 2009 at 3:02 PM, Jörg Groß wrote:
> Hi,
>
> I have a quit abstract problem, hope someone can help me here.
>
> I have a vector like this:
>
>
> x <- c(1
The data frame is not essential. I was just trying to keep things tidy.
Try this:
nxt <- c(tail(x, -1), NA)
prv <- c(NA, head(x, -1))
x[nxt == 3 & prv == 1]
On Sun, Jan 25, 2009 at 5:53 PM, Jörg Groß wrote:
> is there a way to do that without generating a data.frame?
>
> In my real data, I hav
on 01/25/2009 04:29 PM Jörg Groß wrote:
> Hi,
>
> I have a quit abstract problem, hope someone can help me here.
>
> I have a vector like this:
>
>
> x <- c(1,2,3,4,5,2,6)
> x
>
> [1] 1 2 3 4 5 2 6
>
> now I want to get the number where the previous number is 1 and the next
> number is 3
> (t
Prof Ripely,
Thanks. SUN gnu-iconv package should overwrite the sun version and so
it does. Apparently, it does not work. I built this library from gnu
source with gcc and now it configures and builds but fails on make
check for regressions.
This topic is touched in the manual with some 'blas' and
Try this:
> DF <- data.frame(x, nxt = c(tail(x, -1), NA), prv = c(NA, head(x, -1)))
> DF
x nxt prv
1 1 2 NA
2 2 3 1
3 3 4 2
4 4 5 3
5 5 2 4
6 2 6 5
7 6 NA 2
> subset(DF, nxt == 3 & prv == 1)$x
[1] 2
On Sun, Jan 25, 2009 at 5:29 PM, Jörg Groß wrote:
> Hi,
>
> I have a
Hi,
fish.new <- fish[fish$GeoArea==1 & fish$Month==10,]
HTH,
Stephan
pfc_ivan schrieb:
I am a beginner using this R software and have a quick question.
I added a file into the R called fish.txt using this line.
fish<-read.table("fish.txt", head=T, fill=T)
The .txt file looks like this.
fish[fish$GeoArea == 1 & fish$Month == 10]
Am 25.01.2009 um 23:06 schrieb pfc_ivan:
I am a beginner using this R software and have a quick question.
I added a file into the R called fish.txt using this line.
fish<-read.table("fish.txt", head=T, fill=T)
The .txt file looks like this. Since
Hi,
I have a quit abstract problem, hope someone can help me here.
I have a vector like this:
x <- c(1,2,3,4,5,2,6)
x
[1] 1 2 3 4 5 2 6
now I want to get the number where the previous number is 1 and the
next number is 3
(that is the 2 at the second place)
I tried something with tail(x
Jonathan Greenberg wrote:
Our lab has a lot of different unix boxes, with different hardware, and
I'm assuming (perhaps wrongly) that by setting a per-user package
installation directory, the packages will only work on one type of
hardware. Our systems are all set up to share the same home dir
On 25-Jan-09 22:29:25, Jörg Groß wrote:
> Hi,
> I have a quit abstract problem, hope someone can help me here.
> I have a vector like this:
>
> x <- c(1,2,3,4,5,2,6)
> x
> [1] 1 2 3 4 5 2 6
>
> now I want to get the number where the previous number is 1 and the
> next number is 3
> (that is the
dear list,
below is an edited version of my response to an r user asking me for
explaining some issues related to r's evaluation rules. i find the
problem interesting enough to be forwarded to the list, hopefully for
comments from whoever may want to extend or correct my explanations.
(i'd like
is there a way to do that without generating a data.frame?
In my real data, I have a big data.frame and I have to compare over
different columns...
Am 25.01.2009 um 23:42 schrieb Gabor Grothendieck:
Try this:
DF <- data.frame(x, nxt = c(tail(x, -1), NA), prv = c(NA, head(x,
-1)))
DF
Ah, perfect -- so would the "ideal" R_LIBS_USER setting (to more or less
guarantee the libraries will work on every possible computer) be
something along the lines of:
~/myRlibraries/%V%p%o%a
Or is this overkill?
--j
Prof Brian Ripley wrote:
On Sun, 25 Jan 2009, Henrik Bengtsson wrote:
Th
see
?subset
Or use indexing, which is covered in section 2.7 of an introduction to
R (but note that a data frame has 2 dimensions)
hth,
Kingsford Jones
On Sun, Jan 25, 2009 at 3:06 PM, pfc_ivan wrote:
>
> I am a beginner using this R software and have a quick question.
>
> I added a file int
Hi,
I have a quit abstract problem, hope someone can help me here.
I have a vector like this:
x <- c(1,2,3,4,5,2,6)
x
[1] 1 2 3 4 5 2 6
now I want to get the number where the previous number is 1 and the
next number is 3
(that is the 2 at the second place)
I tried something with tail(x,
I am a beginner using this R software and have a quick question.
I added a file into the R called fish.txt using this line.
fish<-read.table("fish.txt", head=T, fill=T)
The .txt file looks like this. Since it contains like 30 lines of data I
will copy/paste first 5 lines.
Year GeoArea
I am interested in familiarizing with the functions belonging to paclage
clusterSim.
This package on-line documentation indicates a file pathname which is clear to
me on Linux systems where $VAR
identifies VAR as an environment variable. But what if I am using a Windows
system ?
I do not know
On Sun, 25 Jan 2009, Henrik Bengtsson wrote:
The script .Rprofile evaluates R code on startup. You could use that
to test for various environment variables. Alternatively, use Unix
shell scripts to set system environment variables to be used in a
generic .Renviron. See help(Startup) for more
The script .Rprofile evaluates R code on startup. You could use that
to test for various environment variables. Alternatively, use Unix
shell scripts to set system environment variables to be used in a
generic .Renviron. See help(Startup) for more details.
/Henrik
On Sun, Jan 25, 2009 at 11:22
Our lab has a lot of different unix boxes, with different hardware, and
I'm assuming (perhaps wrongly) that by setting a per-user package
installation directory, the packages will only work on one type of
hardware. Our systems are all set up to share the same home directory
(and, thus, the sam
Does something like this help:
> x <- matrix(runif(25,-2,2), 5)
> x
[,1][,2] [,3] [,4][,5]
[1,] 0.6188957 -1.14716746 1.9046828 -1.9476897 1.96735448
[2,] -0.5872109 -1.48251061 0.9271700 0.8622643 -0.01762569
[3,] -0.9189594 -0.08752786 -0.5730924 -1.58
On Sun, 25 Jan 2009, Karun Gahlawat wrote:
Uwe,
Sorry I missed it. I do have gnu iconv..
SUNWgnu-libiconv
ls -lra /usr/lib/*iconv* | more
lrwxrwxrwx 1 root root 14 Jan 23 21:23 /usr/lib/libiconv.so -> li
bgnuiconv.so
lrwxrwxrwx 1 root root 22 Jan 23 21:23 /usr/lib/
Dear Users R,
If you use the library agricolae, I would like to have a review to improve my
library.
Please, You might fill the satisfaction survey and send to email.
http://tarwi.lamolina.edu.pe/~fmendiburu/survey.htm
Thanks for your response.
Felipe de Mendiburu
http://tarwi.lamolina.edu.p
On 25 January 2009 at 04:39, new ruser wrote:
| Can anyone please refer me to all firms that offer and/or are developing a
| commercially supported version of R for 64 -bit Windows? - Thanks
Try contacting Revolution-Computing.com --- to the best of my knowledge they
expect to have such a produc
On the topic of visualizing correlation, see also
Murdoch, D.J. and Chow, E.D. (1996). A graphical display of large
correlation matrices.
The American Statistician 50, 178-180.
with examples here:
# install.packages('ellipse')
example(plotcorr, package='ellipse')
On Sat, Mar 8, 2008 at 3:01 A
This script is correct,
use library agricolae, graph.freq() is similar hist(), aditional parameters
size<- c(0,10,20,50,100)
f<-c(15,25,10,5)
library(agricolae)
h<-graph.freq(size,counts=f,axes=F)
axis(1,size)
axis(2,seq(0,30,5))
#
# Other function:
# is necesary histogram h with hist() or graph.
Ok,
use library agricolae, graph.freq() is similar hist(), aditional parameters
size<- c(0,10,20,50,100)
f<-c(15,25,10,5)
library(agricolae)
h<-graph.freq(size,counts=f,axes=F)
axis(1,x)
axis(2,seq(0,30,5))
Other function:
# is necesary histogram h with hist() or graph.freq()
h<-graph.freq(x,cou
At 5:42 AM -0800 1/23/09, Josh B wrote:
Hello,
I was wondering if someone could tell me how to output, to file, the
residuals from a REML model-fit. The type of residuals I am
interested in are the simple "original raw values - model fit" type.
?residuals
To find out how to get residuals f
I am beta test XLSolutions commercially supported version of R for 64-bit
windows.
www.xlsolutions-corp.com
--- On Sun, 1/25/09, stephen sefick wrote:
> From: stephen sefick
> Subject: Re: [R] commercially supported version of R for 64 -bit Windows?
> To: newru...@yahoo.com
> Cc: r-help@r-p
I am beta test XLSolutions commercially supported version of R for 64-bit
windows.
www.xlsolutions-corp.com
--- On Sun, 1/25/09, stephen sefick wrote:
> From: stephen sefick
> Subject: Re: [R] commercially supported version of R for 64 -bit Windows?
> To: newru...@yahoo.com
> Cc: r-help@r-p
Many thanks to both Drs. Bates and Fox for the help!
I also figured out yesterday what Dr. Fox just said regarding the
interpretations of those coefficients for a balanced design. Thanks
Dr. Bates for the suggestion of using solve(cbind(1, contr.sum(4))) to
sort out the factor level effects. Model
Karun Gahlawat wrote:
Uwe,
Sorry I missed it. I do have gnu iconv..
SUNWgnu-libiconv
ls -lra /usr/lib/*iconv* | more
lrwxrwxrwx 1 root root 14 Jan 23 21:23 /usr/lib/libiconv.so -> li
bgnuiconv.so
lrwxrwxrwx 1 root root 22 Jan 23 21:23 /usr/lib/libgnuiconv.so ->
.
Uwe,
Sorry I missed it. I do have gnu iconv..
SUNWgnu-libiconv
ls -lra /usr/lib/*iconv* | more
lrwxrwxrwx 1 root root 14 Jan 23 21:23 /usr/lib/libiconv.so -> li
bgnuiconv.so
lrwxrwxrwx 1 root root 22 Jan 23 21:23 /usr/lib/libgnuiconv.so ->
../gnu/lib/libiconv.so
And
Why? Revolution Computing may do what you want.
On Sun, Jan 25, 2009 at 7:39 AM, new ruser wrote:
>
>
>
>
> Can anyone please refer me to all firms that offer and/or are developing a
> commercially supported version of R for 64 -bit Windows? - Thanks
>
>
>
>
>
>[[alternative HTML versio
Dear Doug and Gang Chen,
With balanced data and sum-to-zero contrasts, the intercept is indeed the
general mean of the response; the coefficient of a1 is the mean of the
response in category a1 minus the general mean; the coefficient of a1:b1 is
the mean of the response in cell a1, b1 minus the ge
Although the message is rather unreadable, I guess you want to look at
?reshape.
Uwe Ligges
oscar linares wrote:
Dear Rxperts,
I would like to convert the following:
StudyStudy.NameParameterDestSrcFormValueMin
MaxFSD
1NT_1-0BFK(03)03A0.
On Fri, Jan 23, 2009 at 4:58 PM, Gang Chen wrote:
> With the following example using contr.sum for both factors,
>
>> dd <- data.frame(a = gl(3,4), b = gl(4,1,12)) # balanced 2-way
>> model.matrix(~ a * b, dd, contrasts = list(a="contr.sum", b="contr.sum"))
>
> (Intercept) a1 a2 b1 b2 b3 a1:
Dear Rxperts,
I would like to convert the following:
StudyStudy.NameParameterDestSrcFormValueMin
MaxFSD
1NT_1-0BFK(03)03A0.128510.0E+001.
0.41670E-01
1NT_1-0BFL(00,03)0003D0.36577
1NT_1-0BFL(00
Thanks for the suggestion replacing the leading 0 with a space instead
of nothing to preserve the layout, and for explaining why there is no
option for this.
Yes, I see why this sounds like a bad idea. The reason I asked is that
I use Sweave to write statistical reports, and I like to get the
form
Dear Friends,
Thanks to Rolf Turner, Brian Ripley and Patrick Burns for their
answers.They don't quite resolve the problem, which I now realize is
due to non-standard behavior of JGR, at least on my machine (I
verified that Mac GUI works entirely as expected):
My installation**
Thanks to Ken and Bernardo for their attempts to answer my question,
but I was apparently unclear as to what I meant by "computational
neuroscience".
The tools Ken and Bernardo suggest provide means to analyze data from
neuroscience research, but I'm actually looking for means to simulate
biologic
Can anyone please refer me to all firms that offer and/or are developing a
commercially supported version of R for 64 -bit Windows? - Thanks
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.et
At 14:55 23/01/2009, David Freedman wrote:
Hi - wouldn't it be possible to bootstrap the difference between the fit of
the 2 models? For example, if one had a *linear* regression problem, the
following script could be used (although I'm sure that it could be
improved):
There are a number of m
Mike Lawrence thatmike.com> writes:
> I've noticed that many computational neuroscience research groups use
> MATLAB. While it's possible that MATLAB may have some features
> unavailable in R, I suspect that this may instead simply be a case of
> costly tradition, where researchers were taught MAT
Check out the source to helloJavaWorld package or one of the other
packages that uses rJava. Some of them are:
CADStat Containers JGR RFreak RJDBC RKEA RLadyBug
RWeka gWidgetsrJava helloJavaWorld iplots openNLP rSymPy
rcdk rcdklibs wordnet
On Sun, Jan 25, 2009 at 5:48 AM, cameron.bracken
wrote:
Have you tried c() from the latticeExtra package?
It worked for me (see below)
library(grid)
library(lattice)
x <- seq(0, 10, length=100)
y <- sin(x)
y2 <- 10*sin(x)
f <- rep(c("1", "2"), each=50)
p1 <- xyplot(y~x,groups=f, ylab="BIG LABEL",
# auto.key=list(space="right"),
par.settings = list
If your purpose is simply to represent a correlation matrix it in a more
compact way see ?symnum, the corrgram package and an example in the
book Multivariate Data Visualization (regarding which gives a lattice
implementation).
On Fri, Mar 7, 2008 at 2:15 PM, Martin Kaffanke
wrote:
> Thank you, t
See:
https://stat.ethz.ch/pipermail/r-help/2008-September/173522.html
On Sat, Jan 24, 2009 at 11:01 PM, pluribus wrote:
> I need to create a vector of dates, weekdays only for a function I am
> working on. Thanks to the chron library, I have managed to accomplish
> this, but is there is a better
I want to call the jar file from R. I want to be able to do this without
using system(). Normally a command line call would look like
java -jar eps2pgf.jar -m directcopy myfile.eps
My question is: Can I call this program using the rJava package or any other
(command line options and all)? I
On Fri, 2009-01-23 at 08:53 -0400, Mike Lawrence wrote:
> Hi all,
>
> I've noticed that many computational neuroscience research groups use
> MATLAB. While it's possible that MATLAB may have some features
> unavailable in R, I suspect that this may instead simply be a case of
> costly tradition, w
Dear R-help,
I am creating a two lattice plots (a densityplot() and xyplot()) that
have the same x-axes and then 'printing' them onto the same page, one
above the other (see end of email for an example to generate the graph).
With different labels on the y-axis for each plot the left spacing i
I need to create a vector of dates, weekdays only for a function I am
working on. Thanks to the chron library, I have managed to accomplish
this, but is there is a better / easier way. This is what I have thus
far.
range.dates <- seq.dates('02/02/2009', '03/13/2009', by =
'days')
Michael Kubovy wrote:
Dear R-helpers,
[...]
(2) If I remember dnorm() and want to be reminded of the call, I also
get a list of pages.
It sounds to me like here you want:
args(dnorm)
Patrick Burns
patr...@burns-stat.com
+44 (0)20 8525 0696
http://www.burns-stat.com
(home of "The R Infe
79 matches
Mail list logo