Hi Chris,
As Jeff mentioned, it is hard to tell what you want (correlations
between sets of coordinates could mean many things it seems like to
me), but here is something that perhaps helps:
## some data (usually nice if you provide this rather than us having
to make something up)
d1 <- cbind(x <
Isn't that one of the main things base R is used for?
Maybe your question is not specific enough. Can you provide some sample data
and your best estimate of what the result should look like?
The Posting Guide mentioned below could help you elicit more effective answers.
--
I have been having problems getting recent packages from Berkeley for Windows
and Linux, so I have had to change to UCLA anyway.
---
Jeff NewmillerThe . . Go Live...
DCN:Basic
I think you really have to show use your exact code that you did along
with an 'str' of each intermediate data structure since my quick test
does not bear out what you were saying:
> test.df <- data.frame(a1= 1:6, a2= 1:6, a3 = 1:6, a4 = 1:6, a5 = 1:6, a6 = 1:6
+ , a7=1:6, a8 = 1:6, a9 = 1:6,
almost certainly you restored an old workspace with a conflicting
definition of the c function.
type
c
to confirm,
also type
conflicts(detail=TRUE)
The repair is to remove the conflicting definition
rm(c)
The long term solution is to adopt the recommended practice of not saving
your worksp
Hello r world,
Does anyone know a function or package that can compute correlations between
sets of XY coordinates?
Thanks in advance for your help,
Chris
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read
Does "CNR" ring any bells as a department at Berkeley?The UCLA
mirror is hosted by the statistics department (cran.stat.ucla.edu) so
I am guessing your local mirror is hosted by the CNR department, and
they are probably the people to inform their system is running very
slowly. Maybe the server
My guess is that somewhere you've accidentally redefined "c" to be
another function. Try this:
Open the Terminal and type
R --vanilla
this will start R totally clean and then try
c(4,7,7)
c(1, 7:9)
I bet your issues go away. If that's the case, you can probably fix
things by clearing your old
On Sat, May 5, 2012 at 2:10 PM, Hillary Sardiñas wrote:
> I have been trying to download the latest version to my Macbook X version
> 10.6.8 from my institutions mirror http://cran.cnr.Berkeley.edu
> Usually is only takes minutes, but regardless of the browser or where I
> connect to the internet
Hi,
My best guess is that the c() you want to be calling is not the c()
you are calling. This can happen if it is masked by a function
definition in you global environment or in a package you load. What
happens when you type:
c
at the console? Also, from one of your sessions where it , c(), d
Hmmm, I know it might seem counterintuitive, but try a different
mirror -- I seem to remember some complaints about the Berkeley mirror
being slow a few weeks back (though I don't have physical evidence
thereof so I apologize for any accidental slander to the Berkeley
mirrorers) -- perhaps UCLA?
h
Ba -- far too much work to recreate (and I don't think you sent us
the file "act.lig"): here's a much better route:
Go to the step immediately before you're in trouble and use dput() on
your data. R will print out a nice plaintext representation that we
can copy and paste and reproduce *exactl
I recently changed to Mac OS X 10.7.3 from Windows and this simple function
is giving errors.
> c(4,7,7)
Error in c(4, 7, 7) : unused argument(s) (7)
> c(1,7:9)
Error in d[i, ] : incorrect number of dimensions
> c(1:5, 10.5, "next")
Error in c(1:5, 10.5, "next") : unused argument(s) ("next")
I would like to ask a general question about the randomForest predict
function and how it handles No Data values. I understand that you can omit
No Data values while developing the randomForest object, but how does it
handle No Data in the prediction phase? I would like the output to be NA
if any
I have been trying to download the latest version to my Macbook X version
10.6.8 from my institutions mirror http://cran.cnr.Berkeley.edu
Usually is only takes minutes, but regardless of the browser or where I
connect to the internet I have been unable to download because it would
take upwards 2 da
Hi,
I had a query regarding which object, a list or a dataframe, consumes more
R memory. Let me clarify this:
For example, I have a df of 6 rows and 12 columns, say 'test'. I do
object.size() and find it uses 3.3 KB of memory.
I run a loop and make a list, say 'testlist', of 6 elements, each ele
Dear all:
I have a question about how to get the optimal estimate of coefficients
using the penalized quantile regression (LASSO penalty in quantile
regression defined in Koenker 2005).
In R, I found both
rq(y ~ x, method="lasso",lambda = 30) and
rq.fit.lasso(x, y, tau = 0.5, lambda = 1, beta = .9
On May 5, 2012, at 2:42 PM, Josh Browning wrote:
Hello useRs!
So, I have a random question. I'm trying to build a character
string, then
evaluate it.
Actually you are trying to build a language object , a call or an
expression. You might have gotten further with:
do.call(ksvm, list(
The alarm() function just sends a ctrl-g (ASCII BEL character) to the console
terminal. If the terminal type recognizes this byte according to that
interpretation, and it is not locally configured to disable sound, then you
should hear an audible alert. Linux consoles typically do handle this. I
> I think that if you're not familiar with emacs or vim I suggest you to
use RStudio. If you're a programmer that knows how to use eclipse I think
eclipse with statet it will be the best solution. Another option that you
have is the notepad + + with NppToR plugin
>
> Best,
>
Att,
Leandro Marino
-
What's the "big picture" of what you're trying to do? eval(parse(text
= )) is often a less than optimal idea.
Some guesses:
Are you trying to construct a formula object (in the strict sense of
something that you pass to a modeling function)?
Maybe lazy evaluation of the deparse(substitute(x)) fl
Or more generally, see the overview page at:
http://www.sciviews.org/_rgui/projects/Editors.html
Uwe Ligges
On 05.05.2012 20:44, R. Michael Weylandt wrote:
There's a huge spectrum -- all the way from Vim / Emacs to notepad:
one that's designed with R in mind (though much heavier than notepad)
1. You need to call the function
pru()
after defining it.
2. Next excercise will be to get rid of all those loops.
3. Next excercide will be to reutn something useful from your function
(rather than printing).
4. Upgrade to a recent version of R.
Uwe Ligges
On 05.05.2012 17:31, Trying
There are many, off the top of my head:
--Emacs + ESS (what I use). You can get them separately or combined
from: http://vgoulet.act.ulaval.ca/en/emacs/
--RStudio
--Vim
--Eclipse + StatEt
Most of those work on different (all?) platforms.
Cheers,
Josh
On Sat, May 5, 2012 at 8:34 AM, Trying To
Hello useRs!
So, I have a random question. I'm trying to build a character string, then
evaluate it. I think an example would be the easiest way to explain:
kern.vec = c("rbfdot","polydot")
for( j in 1:length( kern.vec ) )
{
formula= paste("ksvm( ind ~ . ,
data=d.temp[,c(
There's a huge spectrum -- all the way from Vim / Emacs to notepad:
one that's designed with R in mind (though much heavier than notepad)
is RStudio, which is proving pretty popular in the R community:
http://rstudio.org/ If you get into more programming, it might be
worth it to adopt a general IDE
Dear users of metafor,
I am working on a meta-analysis using the metafor package. I have a
excel csv database that I am working with. I am interested in pooling
the effect measures for a particular subgroup (European women) in this
csv database. I am conducting both sub-group and meta-regression.
Hi all,
I´m using the windows writting pad (not Notepad the simplest version). I
think there should be greater and helping note pads to help programming.
Can you suggest one?
Many thanks
[[alternative HTML version deleted]]
__
R-help@r-pr
On Sat, May 5, 2012 at 8:34 PM, David Winsemius wrote:
> Is there a command you could call with system()?
>
Yes. Something like:
system("play /tmp/02Canon.mp3")
works jsut fine. But I was curious about alarm().
Liviu
__
R-help@r-project.org mailing li
I´m using this version of R
R version 2.13.2 (2011-09-30)
Copyright (C) 2011 The R Foundation for Statistical Computing
ISBN 3-900051-07-0
Platform: i386-pc-mingw32/i386 (32-bit)
2012/5/5 Trying To learn again
> Hi all,
>
> I´m re-starting (as my name indicates) my little knowlegde of CRAN R.
Hi all,
I´m re-starting (as my name indicates) my little knowlegde of CRAN R.
I made this function time before but I don´t know where is the error
because nothing appears as an error but the histogram plot doesn´t appear.
Should I install some special library to run sapply?
pru<-function(){
ra
Dear all,
the commented source code for the ROCR package
(http://cran.r-project.org/web/packages/ROCR) is now available on
github -- feel free to fork, add improvements, and contribute back!
https://github.com/ipa-tys/ROCR
Kind regards,
Tobias
___
R-
Hi,
library(pracma)
k=20
mu=4.5
casigma=17000
myint=function(j) {
quadinf(function(x)
(1/(1+exp(-x)))^j*(1-1/(1+exp(-x)))^(k-j)*dnorm(x,mu,casigma),-Inf,Inf)
}
sapply(0:k,myint)
##
Liviu;
Is there a command you could call with system()?
--
David.
On May 5, 2012, at 2:12 PM, Liviu Andronic wrote:
Dear all
I'd like to make a beeping sound in R, but alarm() doesn't beep? I
checked ?alarm but I couldn't find any pointers to system
configuration. Any ideas?
Regards
Liviu
It doesn't seem to be super robust and might depend on how you are
running R: for me (OS X 10.6), I get the beep at the terminal and in
the R.app GUI, but not within RStudio.
Michael
On Sat, May 5, 2012 at 2:12 PM, Liviu Andronic wrote:
> Dear all
> I'd like to make a beeping sound in R, but ala
Dear all
I'd like to make a beeping sound in R, but alarm() doesn't beep? I
checked ?alarm but I couldn't find any pointers to system
configuration. Any ideas?
Regards
Liviu
> sessionInfo()
R version 2.14.2 (2012-02-29)
Platform: x86_64-pc-linux-gnu (64-bit)
locale:
[1] LC_CTYPE=en_US.UTF-8
Hi,
I am a little confused at the output from predict() for a zeroinfl object.
Here's my confusion:
## From zeroinfl package
fm_zinb2 <- zeroinfl(art ~ . | ., data = bioChemists, dist = "negbin")
## The raw zero-inflated overdispersed data
> table(bioChemists$art)
0 1 2 3 4 5 6
On 03.05.2012 12:39, renu.s7 wrote:
Hi All,
I have a doubt. I used macros and i try to pass a value to a macro by
concatenating a bunch of strings. But it does not seem to work. Please help.
I have written down my code and the error message please tell me how to pass
the value that a string po
On 04.05.2012 14:34, Johannes Radinger wrote:
Hi,
I have two variables ranging both from 0 to 1 (n=500 each).
Now I am interested in plotting them both in one plot (using ggplot2).
So far I used ecdf() (from an example I found with google) to get
values for the cumulatice distribution functio
On 04.05.2012 13:54, kiinalist wrote:
Hi,
I tried to use naivebayes in package 'e1071'.
when I use following parameter, only one predictor, there is an error.
m<- naiveBayes(iris[,1], iris[,5])
table(predict(m, iris[,1]), iris[,5])
Error in log(sapply(attribs, function(v) { :
Non-numeric ar
On May 4, 2012, at 4:22 PM, Petr Savicky wrote:
On Fri, May 04, 2012 at 07:43:32PM +0200, Kehl Dániel wrote:
Dear Petr,
thank you for your input.
I tried to experiment with (probably somewhat biased) truncated means
like in the following code.
How I got the 225 as a truncation limit is a good
And you don't want to know about some of the other problems with the FBI’s
Uniform Crime Reports. IRRC, they are fine for what the FBI intended but a lot
of reseachers don't read the data descriptions as closely as they should.
John Kane
Kingston ON Canada
> -Original Message-
> From:
And you don't want to know about some of the other problems with the FBI’s
Uniform Crime Reports. IRRC, they are fine for what the FBI intended but a lot
of reseachers don't read the data descriptions as closely as they should.
John Kane
Kingston ON Canada
> -Original Message-
> From:
"Relative paths are not a mystery nor are they solely an aspect of linux.
> They are in fact the norm in DOS and Windows as well as in Linux or
> any other file system that use a "tree" structure."
Oh, is that all! I just never thought of that behaviour as 'relative paths'
just as something tha
On May 4, 2012, at 3:00 PM, Diego Rojas wrote:
Thanks, I know about it but i wat to find several local maxima, so
in other words I need a way to identify the places in the surface
where both slopes are equal to 0 and the second derivative is
negative.
There is no way that I know that wil
On Dec 10, 2010, at 9:27 AM, Scott Chamberlain wrote:
Dear R Community,
I recently switched to a Mac (10.6.5), and have installed
Textwrangler to
run code to R. However, I can't install the syntax highlighting file
because
I can't find the directory: "~Users/username/Library/Application
S
gaiarrido usal.es> writes:
> When I try to adjust a mixed model with random effects I can make this order
> without problem
>
> > lm.FA<-lme(absFA~trait*condition,random=~1|individual)
>
> But if I try to fit a model in which the response (absFA) is not the same in
> all individuals at differen
agent dunham hotmail.com> writes:
> I'd like to fix a mixed model. I have unbalance data, what should i use:
> lme in nlme package , or lmer in lme4.
>
> Thanks, user host.com as user host.com
>
More advanced mixed model questions belong on r-sig-mixed-models
r-project.org , but the answ
Perhaps you could contact the persons that supplied/created the file and
ask them what the format of the file exactly is. That is probably the
safest thing to do.
If you are sure that the lines containing only whitespace are
meaningless, then you could alter the previous code to make a copy
Hi All,
Can the MNP package available in R be used to analyze panel data as well?
*i.e., *if there are 3 observed discrete choices for three time periods for
the same individual , can i estimate a panel multinomial probit model which
allows correlated errors across time periods and individual het
Estaré ausente de la oficina desde el 04/05/2012 y no volveré hasta el
14/05/2012.
Responderé a su mensaje cuando regrese.
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http
Hi Johannes ,
Try this:
var1 <- rep(c("a","b"),c(100,100))
var2 <- runif(200,1,50)
df.test <- data.frame(var1,var2)
fn <- function(x){
x <- x$var2
x1 <- sort(x)
x2 <- seq(length(x))
x3 <- x2/max(x2)
df <- data.frame(x1,x2,x3)
df
}
require(plyr)
ddply(d
On Sat, May 05, 2012 at 09:21:10AM +0300, kiinalist wrote:
> Hi,
>
> I tried to use naivebayes in package 'e1071'.
> when I use following parameter, only one predictor, there is an error.
>
> > m <- naiveBayes(iris[,1], iris[,5])
> > table(predict(m, iris[,1]), iris[,5])
> Error in log(sapply(a
Hi,
I am really confused how ddply work, so maybe you can help me.
I created a function that sorts a vector etc.
fn <- function(x){
x1 <- sort(x)
x2 <- seq(length(x))
x3 <- x2/max(x2)
df <- data.frame(x1,x2,x3)
df
}
Probably this is not the best form of t
Given your starting values, there is nothing to optimize:
> f <- function(x, A=10, b=152, T=100, c=100) A*(1-exp(-b*(x-T) -
c*(sqrt(x) - sqrt(T
> f(time)
[1] -Inf -Inf -Inf -Inf -Inf -Inf -Inf -Inf -Inf -Inf -Inf -Inf -Inf
-Inf -Inf
[16] -Inf -Inf -Inf -Inf -Inf -Inf -Inf -Inf -Inf -Inf -Inf
Dear all,
I have a MA plot for "replicate" single channnel arrays, the loess curve is
aligned to M=0. However, I expect that M values be as close as possible to
M=0 for replicate arrays, which is not the case. How may I obtain a measure
that indicates the average distance of points to M=0?
Than
56 matches
Mail list logo
