similar to: Random Variable Generation

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

2003 Dec 15
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 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
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
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
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
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
2013 Sep 24
[PATCH] 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
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<`(`*tmp*`, value = c("1",
2011 Jan 19
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 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
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/ | head 38.77%
2008 Jan 04
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
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
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
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
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
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
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, ? ? ? ? ? ?+X.urban.suburban+X.CapWks,
2010 Dec 09
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
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
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
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> >To: "R-help at" <R-help at