Hello
On Sat, Nov 13, 2010 at 7:57 AM, 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 principal components, I can
> j
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 pri
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 in
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 =
cont
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 err
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
feasible
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:
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.4366
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
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 1-to-
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
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
> pno
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
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).
C
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 bias
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
p
On Fri, Nov 12, 2010 at 1:27 PM, Ralf B 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 the file if it do
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
BLUP
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
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 http://ww
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
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
caref
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 usef
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 (dif
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
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
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 wrote:
> Hello R-help,
>
> I noticed the following surprising behavior when using %in% to find elements
> in a vector generated using seq().
>
> # weird
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 = log
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
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) : 'origi
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 Thu, Nov 11, 2010 at 11:45, Joon Yeong Kim 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 curve directly from t
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 com
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
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,
> -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
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 i
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
801.408.8
-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 co
SILVIA DIAZ FERNANDEZ 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 using
On 12/11/2010 8:58 AM, Jan Eitel wrote:
Bert Gunter 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-
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
Bert Gunter 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] LC_C
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
wrote:
> Fancy that... vector spam :)
>
> Michael
>
> On 12 November 2010 20:30, Jeff Musgrave wrote:
> > Now you and your Vec
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). I
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$' sh
-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 LDFLA
unlink("D:/ffdata", recursive = TRUE)
On Fri, Nov 12, 2010 at 6:29 AM, Xiaobo Gu 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 file.remov
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
previous
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
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 not
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
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
> Odeslal: r-help-boun...@r-project.org
>
> 12.11.2010 10:49
>
> Re: [R] filling a vector with a tapply function applied to another
vector
>
>
>
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:
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
https://stat.ethz.ch/mailman/
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, ...).
G
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
subpopula
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!
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
> 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
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
_
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 to
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]
Fancy that... vector spam :)
Michael
On 12 November 2010 20:30, Jeff Musgrave 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 ready for this no-i
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 Tich
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 b
On Fri, Nov 12, 2010 at 7:05 AM, Tal Galili 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 - it's just
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 8
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 column.
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
> scopi
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 bac
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 u
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
estima
Hi Diana,
Have a look at the setparts function in the partitions package.
Michael
On 12 November 2010 20:03, Diana 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 get
> {{1,2,3}}
> {1,{2,3}}
>
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 drosoph
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 (thi
You could just create a new data frame with the result and cbind?
On Fri, Nov 12, 2010 at 2:11 PM, albechan 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 need fill the third c
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,
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 column
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
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 inf
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 a
82 matches
Mail list logo
