Displaying 20 results from an estimated 2000 matches similar to: "R-SIG-Debian Digest, Vol 35, Issue 1"
2008 Oct 29
2
Help using tapply with multiple variables
Dear list,
I have the function (as a simple example, which is actually part of a
larger function)
pres.test<-function(N0=N0, N1=N1)
{
dt<-5
r<-log(N1/N0)/dt
r
}
which calculates the annual growth rates in a population
Where N0 is the population classified into age intervals, say 5
years, at time=1995, and N1 is the population by 5 year age classes
at time=2000.
2010 Jan 05
1
Multivariate Poisson GLM??
Dear R Users,
I'm working on a problem where I have a multivariate response vector of
counts and a continuous predictor.
I've thought about doing this the same way you would do a Multvariate
regression model with normally distributed data, but since these data are
counts, they are probably better modeled with a Poisson distribution.
For example
y1<-rpois(100,3.5)
y2<-rpois(100,1.5)
2010 Mar 10
1
Strange result in survey package: svyvar
Hi R users,
I'm using the survey package to calculate summary statistics for a large
health survey (the Demographic and Health Survey for Honduras, 2006), and
when I try to calculate the variances for several variables, I get negative
numbers. I thought it may be my data, so I ran the example on the help
page:
data(api)
## one-stage cluster sample
dclus1<-svydesign(id=~dnum, weights=~pw,
2011 Jan 13
1
Openbugs and rbugs on mac with wine
Hello list,
I?ve been trying to get OpenBUGS running on my mac using the wine
emulator. I can run Openbugs just fine by doing:
wine ~/OpenBUGS312/OpenBUGS.exe
In the terminal, so OpenBUGS works. When I try to run the schools
example using rbugs(), the OpenBUGS process starts in wine, but it
just sits there, no log, no script, no output of any sort. The rbugs
() call makes the init,
2008 Jul 16
4
Likelihood ratio test between glm and glmer fits
Dear list,
I am fitting a logistic multi-level regression model and need to test the difference between the ordinary logistic regression from a glm() fit and the mixed effects fit from glmer(), basically I want to do a likelihood ratio test between the two fits.
The data are like this:
My outcome is a (1,0) for health status, I have several (1,0) dummy variables RURAL, SMOKE, DRINK, EMPLOYED,
2007 Sep 27
1
Reading SAS SD2 Data file
I am trying to read a SAS dataset into R that is in a .SD2 format, which
is the Version 6 standard format from SAS.
I see the routines that read the SAS XPORT format (foreign, Hmisc), but
is there any way to read this one?
Any help would be greatly appreciated.
Corey Sparks
Corey Sparks
Assistant Professor
Department of Demography and Organization Studies
University of Texas-San Antonio
One UTSA
2011 Jun 24
3
Fwd: Extract element of a list based on an index value
> Dear list,
>
> I have some data on a geneaology, here is a subset:
> warmerge[1:11,c(1,6,25)]
> Warrior SibID birth.year
> 1100 3793 2013 1926
> 4 2013 2024 1934
> 1094 3769 2024 1918
> 632 2747 2037 1928
> 176 2083 2039 1944
> 187 2085 2039 1949
> 192 2086 2039 NA
> 495
2009 Feb 25
1
Using gregexpr with multiple search elements
Dear list,
I am trying to use gregexpr to see if entries in a dataframe have
either of two possible values for a string.
here's an example
text<-c("fat", "rat", "cat", "dog", "log", "fish")
If I just wanted to find if any one of the elements in text match the
pattern "at" I would do
gregexpr("\\at", text)
2008 Dec 18
1
svyglm fit method?
Dear list,
I am using svyglm in the survey library to fit a binomial logistic
regression accounting for sample design. The documentation says the
models are not fit by maximum likelihood, so my question is what is
the fitting method? Pseudo likelihood? Generalized least squares?
Thanks for any insight.
Best,
Corey
Corey Sparks
Assistant Professor
Department of Demography and
2007 Dec 05
4
coxme frailty model standard errors?
Hello,
I am running R 2.6.1 on windows xp
I am trying to fit a cox proportional hazard model with a shared
Gaussian frailty term using coxme
My model is specified as:
nofit1<-coxme(Surv(Age,cen1new)~ Sex+bo2+bo3,random=~1|isl,data=mydat)
With x1-x3 being dummy variables, and isl being the community level
variable with 4 levels.
Does anyone know if there is a way to get the standard error
2011 Jan 09
1
question about the chow test of poolability
Good day R-listers,
My question is more a statistical question than an R related question,
so please bear with me
i'm currently applying the chow test of poolability
in fact i'm working with panel N=17 T=5 , and my model looks like
this : Yit= a0+B1X1+B2X2+B3X3+B4X4+eit
My question is the following when i'm Testing for the equality of
the coefficients of the unpooled data (the
2010 Mar 04
1
logistic regression by group?
Hi,
Looking for a function in R that can help me calculate a parameter that
maximizes the likelihood over groups of observations.
The general formula is:
p = exp(xb) / sum(exp(xb))
So, according to the formulas I've seen published, to do this "by group" is
product(p = exp(x_i * b_i) / sum(exp(x_i b_i)))
Where i represents a "group" and we iterate through each group.
2010 Mar 10
1
Displaying political boundaries
?
Dear all,?
?
?
I am new in the R environment. How can I enable political boundaries in R?
I want plot a variable (e.g., temperature) only for Mozambique.
?
?Sincerely,
?
Arlindo
?
?
?
Arlindo
____________________________________________________________________________________
Veja quais s?o os assuntos do momento no Yahoo! +Buscados
http://br.maisbuscados.yahoo.com
2010 Jul 18
2
Neural Network
Hi all,
I am working for my master's thesis and I need to do a neural network to
forecast stock market price, with also external inputs like technical
indicators.
I would like to know which function and package of R are more suitable for
this study.
Thanks a lot for your response,
Arnaud TREBAOL.
--
Arnaud Trébaol
T.I.M.E. Student
Ecole Centrale de Lille (09)
Politecnico di Milano (10)
2010 Jul 25
2
R equivalent of SAS proc freq
Dear R-users,
I am looking for a R function that would be the equivalent of the SAS proc
freq (
http://support.sas.com/documentation/cdl/en/procstat/63104/HTML/default/viewer.htm#/documentation/cdl/en/procstat/63104/HTML/default/procstat_freq_sect006.htm).
The table, ftable, xtabs functions are close but do not quite offer the same
capabilities (e.g. they just return counts and no %ages as far as
2011 Jul 08
1
Visualizing a dissimilarity matrix in Euclidean space
Hi,
I have a set of nodes and a dissimilarity matrix for them, as well as a csv
file in which the diss matrix has been converted to [node_1, node_2,
dissimilarity] format. I would like to visualize this as a graph in
Euclidean space (that is, similar nodes clumped together in clusters),
rather than the seriation visualization given by dissplot(). I am using
Network WorkBench for my
2010 Jul 22
2
Multilevel survival model
* Please cc me if you reply as I am a digest subscriber *
Hi,
I am wondering how I can run a multilevel survival model in R? Below is
some of my data.
> head(bi0.test)
childid famid lifedxm sex age delta
1 22.02 22 CONTROL MALES 21.36893 0
2 13.02 13 MAJOR MALES 21.18001 0
3 64.02 64 CONTROL MALES 20.09377 0
4 5.02 5 CONTROL FEMALES
2013 Jan 15
0
smbldap-tool scripts not executing
I'm running a samba pdc+ldap and am having issue with a windows 7 client.
Checking the logs, I see nothing about smbldap-useradd. If I manually run
smbldap-useradd -w hostname$, I'm able to add the client to the domain, but
then run into the trust issue preventing anyone from logging in.
The PDC was running Samba 3.6 from Sernet, but downgraded to 3.5 as part of
CentOS 6.
As a test, I put
2010 Apr 23
3
Event History Data Recoding
Dear R list,
I have an event history data set that is structured like this:
Legislative act Discussion Agreement Time Event
Act1 2006-05-30 2006-06-19 20 1
Act2 2004-03-01 2004-06-14 105 1
.
.
.
I have information on the meetings in the legislature between adoption
periods in a separate variable (the
2009 Sep 22
2
Subsetting dataframes based on column names
Dear R users,
I am interested in taking the columns from multiple dataframes, the
problem is that the different dataframes have different combinations
of the same variable names, here's a simple example:
a<-rep(1:10)
b<-rep(1:10)
c<-rep(21:30)
d<-rep(31:40)
dat.a<-data.frame(a,b,c,d)
names(dat.a)<-c("a", "b", "c", "d")