Displaying 20 results from an estimated 775 matches similar to: "Random Variable Generation"

2003 Dec 15

2

help in lme

To anyone who can help,
I have two stupid questions, and one fairly intelligent question
Stupid question (1): is there an R function to calculate a factorial of a number? That is...is there a function g(.) such that g(3) = 6, g(4) = 24, g(6) = 720, etc?
Stupid question (2): how do you extract the estimated covariance matrix of the random effects in an lme object?
Intelligent question

2003 Jan 14

1

MULTIPLE IPSEC TUNNELS

I am have a shorewall firewall and freeswan ipsec running on a redhat 8.0
Linux gateway machine. I have one working tunnel defined, all works well. I
am not clear how to define mutiple concurrent tunnels. I can not add further
interface entries as all the tunnels come in on ipsec0, do I still have
mutiple zone definitions? some of the tunnels will be dynamic roadwarriors
and as such would need a

2010 Feb 22

4

Siegel-Tukey test for equal variability (code)

Hi, I recently ran into the problem that I needed a Siegel-Tukey test for
equal variability based on ranks. Maybe there is a package that has it
implemented, but I could not find it. So I programmed an R function to do
it. The Siegel-Tukey test requires to recode the ranks so that they express
variability rather than ascending order. This is essentially what the code
further below does. After the

2012 Nov 23

3

Add column index number to str() output

I often find it would be very useful when inspecting a data structure using str() to know the column index number in order to rearrange the data in a manner amenable to my purpose. Is there a way to modify the display options of the str() function to add column index? I know this is really a trivial matter but it would increase the utility of the str() function for me and increase the speed with

2006 Apr 29

2

Creating your own params

I want to pass my own param called randid so I can do this
@rnd = params[:randid]
this will work if I can pass the :randid
--
Posted via http://www.ruby-forum.com/.

2013 Sep 24

19

[PATCH] curve25519-sha256@libssh.org key exchange proposal

Dear OpenSSH developers,
I''ve worked this week on an alternative key exchange mechanism, in
reaction to the whole NSA leaks and claims over cryptographic backdoors
and/or cracking advances. The key exchange is in my opinion the most
critical defense against passive eavesdropping attacks.
I believe Curve25519 from DJB can give users a secure alternative to
classical Diffie-Hellman

2008 Dec 05

2

Help with wavCWTPeaks

I cannot understand the following error printed out when I try to get the extrema of my time series.
I would appreciate some suggestion as I really cannot interpret the error. I might not be using a proper
set of parameters in calling such functions. I am learning by doing ...
> aa.peak <- wavCWTPeaks (aa.tree)
Error in `row.names<-.data.frame`(`*tmp*`, value = c("1",

2011 Jan 19

4

Using subset to filter data table

I am having difficulty understanding how I would constrain a data set by
filtering out ''records'' based on certain criteria.
Using SQL I could query using ''select * from my.data where LithClass in
(''sand'', ''clay'')'' or some such.
Using subset, there seem to be ghosts left behind (that is, all of the
LithClass *.Labels* remain

2017 Dec 01

3

gnu X sysv hash performance

I got curious how the lld produced gnu hash tables compared to gold. To
test that I timed "perf record ninja check-llvm" (just the lit run) in a
BUILD_SHARED_LIBS build.
The performance was almost identical, so I decided to try sysv versus
gnu (both produced by lld). The results are interesting:
% grep -v '^#' perf-gnu/perf.report-by-dso-sym | head
38.77% ld-2.24.so

2008 Jan 04

3

GLMMs fitted with lmer (R) & glimmix (SAS)

I''m fitting generalized linear mixed models to using several fixed effects (main effects and a couple of interactions) and a grouping factor (site) to explain the variation in a dichotomous response variable (family=binomial). I wanted to compare the output I obtained using PROC GLIMMIX in SAS with that obtained using lmer in R (version 2.6.1 in Windows). When using lmer I''m

2007 Apr 12

0

LME: incompatible formulas for groups

Dear R-Users,
I am currently working with LME to analyse repeated measures data. I encounter a problem when including both a random effect and a correlation structure with different grouping levels into the LME model. The error message is:
Error in lme.formula(diameter ~ flowers*timef + competition*timef + population*timef, :
Incompatible formulas for groups in "random" and

2012 Dec 31

5

cut ()

Hello List,
My goal is to create a 5 category variable (p1_st_data$ob_mrj_cat), based on the p1_st_data$obt_mrj_p variable, using the following code for 50 States and District of Columbia (N=51).
p1_st_data$ob_mrj_cat <- cut (p1_st_data$obt_mrj_p, quantile (p1_st_data$obt_mrj_p, (0:5/5), include.lowest=TRUE))
The issue is that, for Utah, I am getting an <NA> instead of (42,48.7]

2009 Feb 26

11

OpenSSL::Cipher decrypt returns 'wrong final block length'

I just want to encrypt a string submitted through a form before saving
it to the DB. And then decrypt it again when I need to retrieve and
use it.
Im trying to use the OpenSSL::Cipher library. I have the following
module for encryption/decryption
[code]
require ''openssl''
module AESCrypt
# Decrypts a block of data (encrypted_data) given an encryption key
# and an

2007 Dec 08

0

help for segmented package

Hi,
I am trying to find m breakpoints of a linear regression model. I
used the segmented package. It works fine for small number of
predicators and breakpoints.(3 r.v. 3 points). However, my model has
14 variables it even would not work even for just one breakpoints!.
The error message is always estimated breakpoints are out of range.
Since my problem is time related problem. So I

2017 Jun 12

1

plotting gamm results in lattice

Hi Maria
If you have problems just start with a small model with predictions and then plot with xyplot
the same applies to xyplot
Try
library(gamm4)
spring <- dget(file = "G:/1/example.txt")
str(spring)
'data.frame': 11744 obs. of 11 variables:
$ WATERBODY_ID : Factor w/ 1994 levels "GB102021072830",..: 1 1 2 2 2 3 3 3 4 4 ...
$ SITE_ID

2017 Jun 12

1

plotting gamm results in lattice

Dear all,?
I hope that you can help me on this. I have been struggling to figure this out but I haven't found any solution.
I am running a generalised mixed effect model, gamm4, for an ecology project. Below is the code for the model:
model<-gamm4(LIFE.OE_spring~s(Q95, by=super.end.group)+Year+Hms_Rsctned+Hms_Poaching+X.broadleaved_woodland? ? ? ? ? ? ?+X.urban.suburban+X.CapWks,

2010 Dec 09

5

Extremely poor write performance, but read appears to be okay

Hello,
I''m writing from the otherside of the world from where my systems are,
so details are coming in slow. We have a 6TB OCFS2 volume across 20 or
so nodes all running OEL5.4 running ocfs2-1.4.4. The system has worked
fairly well for the last 6-8 months. Something has happened over the
last few weeks which has driven write performance nearly to a halt.
I''m not sure how to

2010 Jul 22

0

Please advise acf and pacf in order to determine order of Arima

I have data as below.Please let me know how the ACF and Pacf used to
determine the order od arima model.
Is there any rules need to be followed to determine order.Please advise
> turkey.price.ts
Jan Feb Mar Apr May Jun Jul Aug Sep Oct Nov Dec
2001 1.58 1.75 1.63 1.45 1.56 2.07 1.81 1.74 1.54 1.45 0.57 1.15
2002 1.50 1.66 1.34 1.67 1.81 1.60 1.70 1.87 1.47 1.59 0.74 0.82

2004 May 26

0

2 way repeated measures ANOVA using R: syntax and reporting question

Dear Friends
I have a technical question about conducting 2 way repeated measures ANOVA analysis using R.
1. Data set: repeated measurement of activity over night (2 hr. intervals)
repeated (within subject)factor: Hours
Between subject: Species, Sex
Dependent variables: specimens
Here is how the data arranged for the analysis:
Subject Replicate Hour Sp. Specimens
1 1 1 a
2 1 1 a

2004 May 26

0

2 way repeated measures ANOVA using R: syntax and reportingquestion

Dear Gideon
Unless you are willing to accept uncorrelated measures in time and
homogeneity of variance I would go for a linear mixed model instead. You
might be surprised on how diffrent can be your results. Try;
>?lme
I hope that this helps
Francisco
>From: GIDEON WASSERBERG <wasserberg at wisc.edu>
>To: "R-help at lists.R-project.org" <R-help at