Displaying 20 results from an estimated 1000 matches similar to: "Power calculation with measurement error"
2007 Sep 17
1
Create correlated data with skew
Hi all,
I understand that it is simple to create data with a specific
correlation (say, .5) using mvrnorm from the MASS library:
> library(MASS)
> set.seed(1)
>
> a=mvrnorm(
+ n=10
+ ,mu=rep(0,2)
+ ,Sigma=matrix(c(1,.5,.5,1),2,2)
+ ,empirical=T
+ )
> a
[,1] [,2]
[1,] -1.0008380 -1.233467875
[2,] -0.1588633 -0.003410001
[3,] 1.2054727 -0.620558768
2007 Jul 13
2
Suggestion to extend aggregate() to return multiple and/or named values
Hi all,
This is my first post to the developers list. As I understand it,
aggregate() currently repeats a function across cells in a dataframe
but is only able to handle functions with single value returns.
Aggregate() also lacks the ability to retain the names given to the
returned value. I've created an agg() function (pasted below) that is
apparently backwards compatible (i.e.
2007 Sep 27
3
Aggregate factor names
Hi all,
A suggestion derived from discussions amongst a number of R users in
my research group: set the default column names produced by aggregate
() equal to the names of the objects in the list passed to the 'by'
object.
ex. it is annoying to type
with(
my.data
,aggregate(
my.dv
,list(
one.iv = one.iv
,another.iv = another.iv
,yet.another.iv = yet.another.iv
)
2007 Oct 01
3
optimize() stuck in local plateau ?
Hi all,
Consider the following function:
####
my.func = function(x){
y=ifelse(x>-.5,0,ifelse(x< -.8,abs(x)/2,abs(x)))
print(c(x,y)) #print what was tested and what the result is
return(y)
}
curve(my.func,from=-1,1)
####
When I attempt to find the maximum of this function, which should be
-.8, I find that optimize gets stuck in the plateau area and doesn't
bother testing the
2008 May 09
1
lme() with two random effects
Hi all,
I have collected response time data from 178 participants ('sub') for
each combination of 4 within-Ss factors ('con','int','tone','cue').
Additionally, I have recorded the gender of each participant, so this
forms a between-Ss factor ('gender'). Normally this would be analyzed
using aov:
2007 Aug 08
3
SWF animation method
Hi all,
Just thought I'd share something I discovered last night. I was
interested in creating animations consisting of a series of plots and
after finding very little in the usual sources regarding animation in
R directly, and disliking the imagemagick method described here
(http://tolstoy.newcastle.edu.au/R/help/05/10/13297.html), I
discovered that if one exports the plots to a
2007 May 24
2
Calculation of ratio distribution properties
Hi all,
Looking to calculate the expected mean and variance of a ratio
distribution where the source distributions are gaussian with known
parameters and sample values are correlated. I see (from wikipedia:
http://en.wikipedia.org/wiki/
Ratio_distribution#Gaussian_ratio_distribution) that this calculation
is quite involved, so I'm hoping that someone has already coded a
function to
2007 Jun 16
0
Fwd: How to set degrees of freedom in cor.test?
You could calculate the confidence interval of the correlation at
your desired df: http://davidmlane.com/hyperstat/B8544.html
The below code takes as arguments the observed correlation, N, and
alpha, calculates the confidence interval and checks whether this
includes 0.
cor.test2=function(r,n,a=.05){
phi=function(x){
log((1+x)/(1-x))/2
}
inv.phi=function(x){
2009 Sep 30
3
Choose IAX or SIP trunking?
Asterisk 1.6 behind a firewall and NAT. We are terminating multiple DID
calls, originating and transferring.
A provider offers both SIP and IAX trunking. Cateris paribus, what is
the preferred solution to choose? What points to consider?
I can name the provider if this is not against this list policy--is it?
Thanks,
-kkm
2007 Aug 30
1
xyplot() groups scope issue
Hi all,
Tinkering with a wrapper for xyplot that will help me plot a bunch of
plots in a data analysis I'm doing and I ran into an odd error that
I'm guessing is a scope issue. Here's a very simple version of the code:
###############
#load lattice
library(lattice)
#create the wrapper function
do.xyplot = function( plot.formula, data.frame, plot.groups){
print(plot.groups)
2007 Oct 06
1
Tricky vectorization problem
Hi all,
I'm using the code below within a loop that I run thousands of times
and even with the super-computing resources at my disposal this is
just too slow. The snippet below takes about 10s on my machines,
which is an order of magnitude or two slower than would be
preferable; in the end I'd like to set the number of monte carlo
experiments to 1e4 or even 1e5 to ensure stable
2007 Jun 20
2
how to create cumulative histogram from two independent variables?
Hi all,
I am extremely newbie to R. Can anybody jump-start me with any clues as to
how do I get a cumulative histogram from two independent variables,
cumhist(X,Y) ?
-jose
[[alternative HTML version deleted]]
2007 Jun 22
1
connecting to running process possible?
Hello,
i'm trying to find a more modern system to reproduce the functionality that
was available through the Histoscope program (from Fermilab). Namely, the
capability of connecting to a running process and having plots update in
realtime in response to new data. Is this possible with R? Thank you,
Charles Cosse
[[alternative HTML version deleted]]
2008 Jul 10
2
Lattice: merged strips?
Hi all,
By default a call to xyplot from the Lattice package when using 2
factors [eg xyplot( dv~iv | XY * AB ) ] yields the following shingle
structure:
|_A_|_A_|_B_|_B_|
|_X_|_Y_|_X_|_Y_|
However, I'm wondering if it is possible to merge the upper shingle
within levels of that factor, as in:
|___A___|___B___|
|_X_|_Y_|_X_|_Y_|
Mike
--
Mike Lawrence
Graduate Student, Department of
2007 Jun 26
3
surprising difference in log()
Hello everybody
My collegue and I noticed a strange behaviour of R on different
platforms. It's a simple computation, but results are rather different.
On Windows XP:
> floor(log(8,2))
[1] 3
which is what one should expect.
Here's instead the result with Mac OS X (same version, 2.5.0
(2007-04-23))
> floor(log(8,2))
[1] 2
Is it a "bug" in R or in the operating
2007 Aug 06
4
Marking and remarking of incoming traffic
I can use DSMARK to mark on the Egress side. Is there a way to
mark/change the DSCP value of an incoming packet on the ingress side?
Thanks.
Jon Flechsenhaar
Boeing WNW Team
Network Services
(714)-762-1231
202-E7
2015 Feb 06
0
PSJIP Leak handle
I have an Asterisk 13 that only processes app Transfer with PJSIP, to the
tune of 60 per second. No voice calls.
After like 2 hours, I can no longer get into Asterisk. This command,
asterisk -r, fails, and also "asterisk -rx core show channels", etc. I am
returned to the bash prompt. I checked the handles and
lsof | grep asterisk |wc -l
7098126
I think there is a kind of handle leak
2007 Jul 26
2
logistic regression
Greetings,
I am working on a logistic regression model in R and I am struggling with the code, as it is a relatively new program for me. In searching Google for 'logistic regression diagnostics' I came Elizabeth Brown's Lecture 14 from her Winter 2004 Biostatistics 515 course (http://courses.washington.edu/b515/l14.pdf) . I found most of the code to be very helpful, but I am
2007 Jul 22
4
using R for a reaction-time experiment
I want to use R to run a reaction-time experiment: Something appears on the
screen, I respond by typing something (one keystroke), the system measures
the speed of my response. R would be great for this if only I didn't have to
hit Enter to enter that keystroke. I am doing such experiments now but they
require two actions per trial: hit keystroke, hit Enter.
Is there some way that R can be
2008 Jul 17
2
Sampling distribution (PDF & CDF) of correlation
Hi all,
I'm looking for an analytic method to obtain the PDF & CDF of the
sampling distribution of a given correlation (rho) at a given sample
size (N).
I've attached code describing a monte carlo method of achieving this,
and while it is relatively fast, an analytic solution would obviously
be optimal.
get.cors <- function(i, x, y, N){
end=i*N