Dear Group,
Is there some way for me to package a few lines of R-code as exe and have it
running in the background? (unable to find info in the archives)
Even better if I can package it as an installation and send to my team who
do not have any background in programming / R? If they can install
Hi R-users,
Since it has been already asked for several times, I would like to share my
solution
that I have gathered from several sources on internet for portable
R+Emacs+MikTeX.
This may also help some future R-users.
My best,
Ozan
Installing R, MikTeX, Emacs (+auctex+aspell) on USB
Useful
Hi Assaf,
I ran the following code in R in order to get a smooth function.
IRfmonfd - with(DFt1, smooth.monotone(term,DFt1[,day], gr.fdPar))
However, I get an error notice. Do you have rough idea about what's going
on?
Error in eval(substitute(expr), data, enclos = parent.frame()) :
numeric
I have a data frame with three columns. The first one is filled with factors,
the second one with numeric values and the third one is an empty vector.
I need fill the third column with the results of the function
tapply(secondcolumn,firstcolumn,mean) so that I have for each row of the
third
Hi Michael,
I had already seen your solution (that does work).
If you're right about the issue in my function, then the error message
is confusing ('could not find the function get- '). Moreover, I
assign()ed and get() x from the .GlobalEnv, so there shouldn't be a
problem with scoping,
You could just create a new data frame with the result and cbind?
On Fri, Nov 12, 2010 at 2:11 PM, albechan alberto.case...@satt.biz wrote:
I have a data frame with three columns. The first one is filled with factors,
the second one with numeric values and the third one is an empty vector.
I
Hi
I am new on this forum. I am searching for a function in R which provides
all partitions of a set, say for the set
{1,2,3}
you get
{{1,2,3}}
{1,{2,3}}
{2,{1,3}}
{3,{1,2}}
{{1},{2},{3}}
. The number of partitions of a set is given by Bellsche`s number. The
number of possibilities of length
Hi everybody,
I have a problem when trying to do the quality control with the packages
simpleaffy and affyQCReport with the drosophila chip 2.0
At first I got the messeage, that the *.qcdef file is not there. I followed
the instructions in tha manual and created the file like that:
array
Hi Diana,
Have a look at the setparts function in the partitions package.
Michael
On 12 November 2010 20:03, Diana d.ti...@mathematik.uni-wuerzburg.de wrote:
Hi
I am new on this forum. I am searching for a function in R which provides
all partitions of a set, say for the set
{1,2,3}
you
Dear Scott,
The short answer is: No, this is not possible.
The longer answer is: In principle, this will be possible in the future. I am
working on functions for the metafor package for multivariate meta-analyses
(i.e., where we no longer assume that the observed outcomes or effect size
On 11/11/2010 10:29 PM, sachinthaka.abeyward...@allianz.com.au wrote:
seems like the texi2dvi doesnt exist on R2.12.0 anymore?
It's in the tools package. (The tools package is for tools used in R
like installing packages, showing help, etc. The utils package contains
tools meant for public
On 12/11/2010 3:23 AM, Santosh Srinivas wrote:
Dear Group,
Is there some way for me to package a few lines of R-code as exe and have it
running in the background? (unable to find info in the archives)
Even better if I can package it as an installation and send to my team who
do not have any
Hi Ivan,
I had already seen your solution (that does work).
If you're right about the issue in my function, then the error message is
confusing ('could not find the function get- '). Moreover, I assign()ed
and get() x from the .GlobalEnv, so there shouldn't be a problem with
scoping, right?
The first column has 331 rows, the number of factors in the second column is
20 repeated in order to have 331 rows too. What I need is to reapeat in the
third column each of the 20 results derived from
tapply(firstcolumn,secondcolumn,mean) as many times as each factor appears
in the second
Hi:
Look into the ave() function - here's a small demonstration:
d - data.frame(g = factor(rep(LETTERS[1:5], each = 5)),
x = rpois(25, 10))
d$mean - with(d, ave(x, g, FUN = mean))
head(d, 10)
head(d, 10)
g x mean
1 A 5 8.4
2 A 14 8.4
3 A 10 8.4
4 A 2 8.4
5 A 11
On Fri, Nov 12, 2010 at 7:05 AM, Tal Galili tal.gal...@gmail.com wrote:
Hello Barry,
With regards to the identical == FALSE, it didn't occur to me - great point,
thank you.
So question: how did you end up changing the -, so to enable the creation
of the .metadata object?
I haven't done that
On 11/12/2010 09:23 AM, Santosh Srinivas wrote:
Dear Group,
Is there some way for me to package a few lines of R-code as exe and have it
running in the background? (unable to find info in the archives)
Even better if I can package it as an installation and send to my team who
do not have any
Hi Diane,
Does this do what you want ?
listParts - function(n) {
# Generates a list of lists representing the partitions
# of an integer n
require(partitions)
x - 1:n
apply(setparts(n), 2, function(pp) tapply(x, pp, function(xx) list(xx)))
}
Michael
On 12 November 2010 21:15, Diana Tichy
Fancy that... vector spam :)
Michael
On 12 November 2010 20:30, Jeff Musgrave j...@teamteach2.com wrote:
Now you and your Vector Team can make more money.
Offer your current client base a chance to buy and sell
a product in high demand. An item that increases in value every day.
Are you
Dear list,
I trying to make a wind rose plot whit the command oz.windrose, from plotrix
package. My data, a matrix of percentages with the rows representing speed
ranges and the columns indicating wind directions was generated
using bin.wind.records command from same package:
[,1] [,2]
Hi!
I am trying to know which habitat variables most affect bird counts in a
radius of 100m. I obtained bird counts in 2751 spatial points, and measured
percentage of 21 habitat variables in these points.
I am applying a mixed model using the “lmer” function to the data, but I do
not know how
Hi,
It seems file.remove can't remove directories, regardless whether
it's an empty one, but how to delete a directory in R
file.remove(D:/ffdata)
[1] FALSE
Warning message:
In file.remove(D:/ffdata) :
cannot remove file 'D:/ffdata', reason 'Permission denied'
Thanks.
Xiaobo Gu
Date: Fri, 12 Nov 2010 11:42:01 +0100
From: rh...@eoos.dds.nl
To: r-help@r-project.org
Subject: Re: [R] Deploying code as exe
On 11/12/2010 09:23 AM, Santosh Srinivas wrote:
Dear Group,
Is there some way for me to package a few lines of
Hi!
I am trying to know which habitat variables most affect bird counts in a radius
of 100m. I obtained bird counts in 2751 spatial points, and measured percentage
of 21 habitat variables in these points.
I applied a mixed model using the lmer function to these data, but I do not
know
Thank you so much! That's exactly what I needed!
much love
--
View this message in context:
http://r.789695.n4.nabble.com/filling-a-vector-with-a-tapply-function-applied-to-another-vector-tp3039271p3039396.html
Sent from the R help mailing list archive at Nabble.com.
Hi,
This is a question regarding technique rather than an R specific issue. I
have been asked to evaluate a 30+ year long term continuous survey of bird
presence/absence data that has an associated ocular estimate of the
vegetation community percent coverage. The data are organized by
Dear R-list
This is just message to inform that the there is an issue with write.fwf in the
gdata library (from version 2.5.0 on). It does not seem to accept further
arguments to write.table like eol as the help file indicates as it stops when
executing tmp - lapply(x, format.info, ...).
Hello.
using the lm command, I always receive the error message couldn't find
the function is. Do you haven an idea, what that means?
Thanks for your help.
Best regards.
Kathrin L.
__
R-help@r-project.org mailing list
Dear Colleagues
I have been struggling with a problem and could not find a similar
previous posting in any of my searches, so I decided to write to the
list. I am trying to produce svg files from plots using RSVGTipsDevice
in R 2.12.0 in Ubuntu 10.04, but I am getting the following error
message:
Hi
Use ave instead of tapply. It is designed for such task
Regards
Petr
r-help-boun...@r-project.org napsal dne 12.11.2010 10:49:17:
albechan alberto.case...@satt.biz
Odeslal: r-help-boun...@r-project.org
12.11.2010 10:49
Re: [R] filling a vector with a tapply function applied to
Dear list,
I'm stuck with looking for a function of the *apply family, which I suppose
exists already – just I can't find it:
What I'm looking for is somewhere between sweep and mapply that does a
calculation vectorized over a matrix and a vector:
It should work complementary to sweep: for
On Nov 12, 2010, at 3:23 AM, Santosh Srinivas wrote:
Dear Group,
Is there some way for me to package a few lines of R-code as exe and
have it
running in the background? (unable to find info in the archives)
Even better if I can package it as an installation and send to my
team who
do
Skew as they are, your data certainly don't look normal. Try lognormal.
The chi-square test gives good results when all counts are 5 or more,
hence the warning.
At 12:25 AM 11/12/2010, Andrew Halford wrote:
Hi All,
I have a dataset consisting of abundance counts of a fish and I want to test
You need to reinstall these packages when you update R. Have you done
update.packages(checkBuilt=TRUE)
? It's the obvious way to find packages which need updating.
On Fri, 12 Nov 2010, leandro wrote:
Dear Colleagues
I have been struggling with a problem and could not find a similar
unlink(D:/ffdata, recursive = TRUE)
On Fri, Nov 12, 2010 at 6:29 AM, Xiaobo Gu guxiaobo1...@gmail.com wrote:
Hi,
It seems file.remove can't remove directories, regardless whether
it's an empty one, but how to delete a directory in R
file.remove(D:/ffdata)
[1] FALSE
Warning message:
In
-BEGIN PGP SIGNED MESSAGE-
Hash: SHA1
Hi
I know this has been asked before, but I can not find it - so my apologies.
I want to compile R myself, to optimize it for speed. There is a small
section (a few lines actually) in the R-admin manual (Compilation-flags)
which states that the
Jannis wrote:
Hi Sachin,
please read the posting guide and include a reproducible example of what
you want to do.
For your first question you should have a look at ?axis. Supplying the
'at' argument with the positions of the desired marks and the 'labels'
with text strings like '10.000$'
Steve:
Obviously, you would do better to consult a local statistician.
Presumably, you have none, so you would probably do better posting on
appropriate R SIGs -- e.g. R-sig-mixed-models, R-sig-ecology and the
liike. R-help is really not designed for this (although good souls
sometimes respond).
that was my bad. I let it in because I saw vector in the subject line, got
lazy and didn't check the contents.
mark
On Fri, Nov 12, 2010 at 5:45 AM, Michael Bedward
michael.bedw...@gmail.comwrote:
Fancy that... vector spam :)
Michael
On 12 November 2010 20:30, Jeff Musgrave
Bert Gunter gunter.berton at gene.com writes:
PLEASE follow the posting guide and give us the output of sessionInfo().
I was unable to duplicate your problem -- it worked fine for me.
sessionInfo()
R version 2.12.0 (2010-10-15)
Platform: i386-pc-mingw32/i386 (32-bit)
locale:
[1]
Dear Prof. Ripley
Thank you for your help. It works perfectly now. I did not realise that
just using update.packages() was not rebuilding the packages under the
new R version.
Best,
Leandro
On Fri, 2010-11-12 at 14:33 +, Prof Brian Ripley wrote:
You need to reinstall these packages when you
On 12/11/2010 8:58 AM, Jan Eitel wrote:
Bert Guntergunter.bertonat gene.com writes:
PLEASE follow the posting guide and give us the output of sessionInfo().
I was unable to duplicate your problem -- it worked fine for me.
sessionInfo()
R version 2.12.0 (2010-10-15)
Platform:
SILVIA DIAZ FERNANDEZ Silvia.Diaz at uclm.es writes:
Hi!
I am trying to know which habitat variables most affect
bird counts in a radius of 100m. I obtained bird
counts in 2751 spatial points, and measured percentage of
21 habitat variables in these points.
I applied a mixed model
-BEGIN PGP SIGNED MESSAGE-
Hash: SHA1
Hi Patrick,
thanks a lot. No, I am not using MKL, but I am going to try the flags
you mention in your blog.
Cheers,
Rainer
On 11/12/2010 03:58 PM, Richardson, Patrick wrote:
Rainer,
The attached link might give you some insight. It's for
In addition to Michael's answer, look at the ?par help page and look at the ask
argument, you can use this to have the graphics device wait for you to click
before starting a new plot.
--
Gregory (Greg) L. Snow Ph.D.
Statistical Data Center
Intermountain Healthcare
greg.s...@imail.org
The best approach if creating all the files using R is to change how you create
the graphs so that they all go to one file to begin with (as mentioned by
Joshua), but if some of the files are created differently (rgl, external
programs), then this is not an option.
One external program that
-Original Message-
From: r-help-boun...@r-project.org
[mailto:r-help-boun...@r-project.org] On Behalf Of Kathrin Linsler
Sent: Friday, November 12, 2010 5:22 AM
To: r-help@r-project.org
Subject: [R] error message: is
Hello.
using the lm command, I always receive the error
Hi Kathrin,
Can you provide:
1) the results of sessionInfo() on your system
2) a copy of the code from the R console leading up to and including
the error message? For example the exact lm() call, and if possible
some example data that you could use in the lm() call.
Sincerely,
Josh
On Fri,
Since I read the list in digest form (and was out ill yesterday) I'm
late to the discussion.
There are 3 steps for predicting survival, using a Cox model:
1. Fit the data
fit - coxph(Surv(time, status) ~ age + ph.ecog, data=lung)
The biggest question to answer here is what covariates you wish
With that limited information I fear that no one can really help you.
Please read the posting guide (link at the bottom of each mail) and
include your code with which we can reproduce your error. Only then we
can guess what causes this problem. For now I can only guess that you
spelled some
On Thu, Nov 11, 2010 at 11:45, Joon Yeong Kim joonyeen...@gmail.com wrote:
I tried using panel function but wasn't sure how to get the y value for the
peak of the densityplot.
I'm thinking there should be a way to retrieve the densityplot object so
that I can get the (x,y) values of the
Hello list,
here is what I stumbled upon:
1 test - data.frame(time=as.POSIXct((1:2),origin=2000-1-1))
1 test
time
1 2000-01-01 00:00:01
2 2000-01-01 00:00:02
1 rbind(test,b=1:2)
Fehler in as.POSIXct.numeric(value) : 'origin' muss angegeben werden
When I try to attach an
On 11/12/2010 9:57 AM, Alexander Senger wrote:
Hello list,
here is what I stumbled upon:
1 test- data.frame(time=as.POSIXct((1:2),origin=2000-1-1))
1 test
time
1 2000-01-01 00:00:01
2 2000-01-01 00:00:02
1 rbind(test,b=1:2)
Fehler in as.POSIXct.numeric(value) : 'origin'
Hello R-help,
I noticed the following surprising behavior when using %in% to find elements in
a vector generated using seq().
# weird!!!
c(7.7, 7.8, 7.9) %in% seq(4, 8, by=0.1)
[1] TRUE FALSE TRUE
# OK now
c(7.7, 7.8, 7.9) %in% round(seq(4, 8, by=0.1), 1)
[1] TRUE TRUE TRUE
# wait, how
I think it is likely I am missing something. Here is a very simple example:
R code:
mat - matrix(nrow = 10, ncol = 2, c(1,0,1,0,0,1,0,0,1,1),
c(5,4,1,6,3,6,5,3,7,9), dimnames = list(c(1,2,3,4,5,6,7,8,9,10),
c(column1,column2)))
g - glm(mat[1:10] ~ mat[11:20], family = binomial (link = logit))
Not weird at all!!! You need to understand how computers do
arithmetic. See R FAQ 7.31.
-- Bert
On Fri, Nov 12, 2010 at 10:55 AM, Vadim Patsalo patsa...@gmail.com wrote:
Hello R-help,
I noticed the following surprising behavior when using %in% to find elements
in a vector generated using
You are not creating your data set properly.
Your 'mat' is:
mat
column1 column2
11 0
21 0
30 1
40 0
51 1
61 0
71 0
80 1
90 0
10 1 1
What you really want
Terry,
My point was that if you are asking the question: What is the average time
to death based on a set of variables? The only logical approach for
calculating actual time to death is to use uncensored cases, because we do
not know the time to death for the censored cases and can only estimate
Dear list
I am calling the functiton bugs() provided by R2WinBugs to performs an IRT
analysis. The function returns a set of estimated parameters over n
replications/iterations. For each replication, two sets of person measures
(theta1 and theta2) and two sets of item difficulty parameters
I know such programs, however, for my specific problem I have an R
script that creates a report (which I have to create many times) and I
would like to append about 100 single paged post scripts at the end as
appendix. File names are controlled so it would be easy to detect
them; I just miss a
You install JGR and its dependencies as you would any other package from the
RGUI. Make sure you have a recent java installed and functional. You may
find jselect.exe on the RForge site useful for looking at this.
When I did the install, things seemed a little flaky, but eventually by
Hi Jan,
The issue isn't that the ... arguments aren't passed on. Rather, the
problem is that in the current implementation the ... arguments are passed
to format(), which doesn't understand the eol argument.
The solution is to modify write.fwf() to explicitly accept all of the
appropriate the
R 2.12 is not functioning for me On the MAC what the most painless way of
reverting
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide
Assuming I would go into the trouble of messing with the existing R
scripts that create mentioned postscripts/pdfs, how can I achieve that
an array of scripts append to a single ps/pdf? I would want the first
script to create the file if it does not yet exist and all other to
append to it with new
Hello!
I wonder how can one get upper and lower limits of a prediction interval
-- exact values? They are shown on caterpillar plot using ranef() with
argument postVar=TRUE, but I would like to know them. A while ago, some
discussions were opened on Confidence Intervals for Random Effect
On Fri, Nov 12, 2010 at 1:27 PM, Ralf B ralf.bie...@gmail.com wrote:
Assuming I would go into the trouble of messing with the existing R
scripts that create mentioned postscripts/pdfs, how can I achieve that
an array of scripts append to a single ps/pdf? I would want the first
script to create
Hi -
I have a dataframe of (x,y) values. I'd like to fit an exponential
curve to the data for further statistical analysis (pretty much the same
functionality provided by Excel's LOGEST worksheet array function). Can
someone point me to the (set of) functions/ package that is best suited to
Jim,
I respectfully disagree, and there is 5 decades of literature to back
me up. Berkson and Gage (1950) is in response to medical papers that
summarized surgical outcomes using only the observed deaths, and shows
important failings of the method. Ignoring the censored cases usually
gives
On Nov 12, 2010, at 4:22 PM, cran.30.miller_2...@spamgourmet.com wrote:
Hi -
I have a dataframe of (x,y) values. I'd like to fit an exponential
curve to the data for further statistical analysis (pretty much the
same
functionality provided by Excel's LOGEST worksheet array function).
On Nov 12, 2010, at 5:07 PM, David Winsemius wrote:
On Nov 12, 2010, at 4:22 PM, cran.30.miller_2...@spamgourmet.com
wrote:
Hi -
I have a dataframe of (x,y) values. I'd like to fit an exponential
curve to the data for further statistical analysis (pretty much the
same
functionality
Hello all,
I have a question about basic statistics. Given a PDF value of 0.328161,
how can I find out the value of -0.625 in R? It is like reversing the dnorm
function but I do not know how to do it in R.
pdf.xb - dnorm(-0.625)
pdf.xb
[1] 0.328161
qnorm(pdf.xb)
[1] -0.444997
On Nov 12, 2010, at 5:35 PM, Edwin Sun wrote:
Hello all,
I have a question about basic statistics. Given a PDF value of
0.328161,
how can I find out the value of -0.625 in R? It is like reversing
the dnorm
function but I do not know how to do it in R.
pdf.xb - dnorm(-0.625)
pdf.xb
Not sure if there's a pre-defined function for it, but use your basic
math skills: the normal distribution is
dnorm(x) = 1/(sqrt(2*pi)) * exp(-x^2/2),
so the inverse function (on the interval [0, infinity] is
f = function(x) {sqrt( -2*log(sqrt(2*pi) * x)) }
Since the dnorm function is not
There is no reversing the dnorm function -- dnorm is many to one in
general (in the normal case, 2 to 1 except for the mean). How would
you reverse dunif, for example?!
However, with that understanding you could do a simple one dimensional
search within the range you want via optimize(), as in
Hi:
This is kind of kludgy, but if the matrix and parallel vector are both
numeric, you could try something like
A - matrix(rnorm(12), nrow = 3)
v - 1:3
f - function(x) c(sum(x[-length(x)]^2), x[length(x)])
t(apply(cbind(A, v), 1, f))
v
[1,] 8.196513 1
[2,] 1.414914 2
[3,] 2.436660
Dear R People:
Is it possible to run R on an iPad, please?
For some reason, I'm thinking that you can't have Fortran, C, etc., so
you can't do it.
But I thought I would check anyway.
Thanks,
Erin
--
Erin Hodgess
Associate Professor
Department of Computer and Mathematical Sciences
University
Hi Erin,
This topic has been fairly well discussed a few months back in these
two threads:
http://www.mail-archive.com/r-help@r-project.org/msg97811.html
http://www.mail-archive.com/r-help@r-project.org/msg107623.html
That said, if R is made available on the cloud, it would of course be
Dear Group,
I'm having lots of problems getting RMySQL on a 64 bit machine. I followed
all instructions available but couldn't get it working yet! Please help.
See the output below.
I did a install of RMySQL binary from the revolution cran source. It seems
to have unpacked fine but gives this
I have been trying to update some R packages but I get the following error.
Can you advise how mow to get around this . I am using the R for 64 bit
windows.
--- Please select a CRAN mirror for use in this session ---
Warning in install.packages(update[instlib == l, Package], l, contriburl =
Hi all,
This could be very basic. I want to do exploratory factor analysis but I
don't have the data, rather I have the correlation matrix. How do I do this
with just the correlation matrix? I know for principal components, I can
just find the eigen values.
Best regards
--
View this message
See ?factanal, and ?ability.cov for a worked example.
On Fri, 12 Nov 2010, Brima wrote:
Hi all,
This could be very basic. I want to do exploratory factor analysis but I
don't have the data, rather I have the correlation matrix. How do I do this
with just the correlation matrix? I know for
Hello
On Sat, Nov 13, 2010 at 7:57 AM, Brima adamsteve2...@yahoo.com wrote:
Hi all,
This could be very basic. I want to do exploratory factor analysis but I
don't have the data, rather I have the correlation matrix. How do I do this
with just the correlation matrix? I know for principal
82 matches
Mail list logo