Displaying 20 results from an estimated 100 matches similar to: "effects in ANCOVA"
2011 May 04
3
Error in .Fortran Call
I have the following FORTRAN code converted to a DLL:
! my_xmean.f90
!
! FUNCTIONS/SUBROUTINES exported from my_function.dll:
! my_function - subroutine
!
subroutine my_xmean(X,N,XMEAN)
! Expose subroutine my_function to users of this DLL
!
!DEC$ ATTRIBUTES DLLEXPORT,C,REFERENCE,ALIAS:'my_xmean_'::my_xmean
! Body of my_function
DOUBLE PRECISION X(N)
XMEAN=0D0
DO J=1,N
2009 Feb 27
1
Sweave doesn't do csv.get()
Hi Everybody
I use R2.8.0 on Mac OS X. I set up LyX 1.6.1 to use Sweave today. I can
compile the test file I found on CRAN (
http://cran.r-project.org/contrib/extra/lyx/) without a problem and the
output looks very nice. In the test file the following R code is used.
<<myFirstChunkInLyX>>=
xObs <- 100; xMean <- 10; xVar <- 9
x <- rnorm(n=xObs, mean=xMean, sd=sqrt(xVar))
2012 May 02
1
returning value from source() in R .Net
Hello
I would like to source an R script from within a C# .Net application equivalent to:
source("my_r_code.r")
I can get this to run but am not sure how to retrieve R objects defined with script my_r_code.r at runtime.
For example, if "my_r_code.r" contains
#-- contents of my_r_code.r---------------------
x <- 1:10
xmean <- mean(x)
2006 Apr 10
1
Generic code for simulating from a distribution.
Hello all,
I have the code below to simulate samples of certain size from a
particular distribution (here,beta distribution) and compute some
statistics for the samples.
betasim2<-function(nsim,n,alpha,beta)
{
sim<-matrix(rbeta(nsim*n,alpha,beta),ncol=n)
xmean<-apply(sim,1,mean)
xvar<-apply(sim,1,var)
xmedian<-apply(sim,1,median)
2011 Aug 07
2
help annotating plot.xmean.ordinaly in rms package
Hello List:
Does anyone know if there's a way to annotate the plots produced by the
plot.xmean.ordinaly function in the rms package. I'd like to produce
publication quality figures, but first need to annotate the axis labels. Is
there a way to do that?
I appreciate any advice.
Chris
[[alternative HTML version deleted]]
2002 Jun 20
2
tabulate means & NA
Hi,
i have some problems to get means, but the data have got
NA's which should not be replaced or delete !
tapply(data,var,mean) # makes me problem,because when one value is
NA mean is NA,too.
Another attempt is this, because the summary does what i want , but
is not subsetable to summary[4] if i use tapply,table and the the
function below gives me
again all means as NA ?
xmean
2012 Mar 29
1
how to increase speed for function?/time efficiency of below function
i am using sarima() function as below
___________________________________________________________________________________________
sarima=function(data,p,d,q,P=0,D=0,Q=0,S=-1,tol=.001){
n=length(data)
constant=1:n
xmean=matrix(1,n,1)
if (d>0 & D>0)
fitit=arima(data, order=c(p,d,q), seasonal=list(order=c(P,D,Q),
period=S),
2004 Dec 02
3
R and Fortran in Windows
I just joined the list and appologize if this has been answered before
but I am trying to interface between R and the Compaq Visual Fortran
compiler version 6.6 for Windows.
I found the following instructions on the web -- and an example. When I
follow these directions exactly. R 2.0.0 crashes. Has anyone had any
experience with this?
Below are the instructions that I located:
Thanks
Dan
2013 Apr 19
2
NAMESPACE and imports
I am cleaning up the rms package to not export functions not to be called
directly by users. rms uses generic functions defined in other packages.
For example there is a latex method in the Hmisc package, and rms has a
latex method for objects of class "anova.rms" so there are anova.rms and
latex.anova.rms functions in rms. I use:
2006 Feb 28
2
Elegant way to express residual calculation in R?
Hi All,
I am illustrating a simple, two-way ANOVA using the following data and I'm
having difficulty in expressing the predicted values succinctly in R.
X<- data.frame(read.table(textConnection("
Machine.1 Machine.2 Machine.3
53 61 51
47 55 51
46 52 49
50
2009 Aug 13
2
Coding problem: How can I extract substring of function call within the function
In order to ease my students into the R environment I am developing a package which installs a variety of utility functions as well as slightly modified versions of some standard R functions -- e.g. mean, hist, barplot, .... In my versions of these standard R functions I either add options or alter some defaults that seem to create difficulties for most of my students -- for example, when they do
2005 Nov 09
2
help with legacy R code
Hi there,
Could somebody help me disect this legacy R script I inherited at work, I
have two questions:
1. I've tried to upgrade our R version from 1.6.2 (yeah, I know), to R 2.0,
but some of the lines in this script are not compatible with R 2.0, could
someone help me figure out where the problem is?
2. the jpeg generated (attached) seems to be off on some of the data, is
there a better way
2006 Aug 02
0
Trying to use segmented in a function
Hi folks
I wonder if anyone can help me. I want to run some simulations to see
how big a sample size might be necessary to distinguish a curved
bivariate relationship (e.g. something that might be best described
by a quadratic model) from a relationship that is two straight lines
with a sudden change in slope (e.g. something best described by a
breakpoint regression). I am using
2008 Oct 20
1
Calculate SPE in PLS package
Dear list,
I want to calculate SPE (squared prediction error) in x-space, can
someone help?
Here are my codes:
fit.pls<-
plsr(Y~X,data=DAT,ncomp=3,scale=T,method='oscorespls',validation="CV",x=
T)
actual<-fit.pls$model$X
pred<-fit.pls$scores %*% t(fit.pls$loadings)
SPE.x<-rowSums((actual-pred)^2)
Am I missing something here?
Thanks in advance.
Stella Sim
2007 May 21
1
PLS in R and SAS
Dear all:
I am comparing the PLS outputs of R and SAS for the following data set:
Y x1 x2 x3
3 6 2 2
3 1 5 5
4 7 4 1
5 6 5 6
2 4 3 2
8 5 0 9
where Y is the dependent variable and x1, x2, x3 are the independent variables. I found several PLS algorithms in R (NIPALS,SIMPLS,KERNEL PLS). SAS has SIMPLS and NIPALS.
The following are the NIPALS calculations of
2012 Nov 23
2
[LLVMdev] [cfe-dev] costing optimisations
On 23.11.2012, at 15:12, john skaller <skaller at users.sourceforge.net> wrote:
>
> On 23/11/2012, at 5:46 PM, Sean Silva wrote:
>
>> Adding LLVMdev, since this is intimately related to the optimization passes.
>>
>>> I think this is roughly because some function level optimisations are
>>> worse than O(N) in the number of instructions.
>>
2012 Jun 25
2
Fractional Factorial - Wrong values using lm-function
Hello.
I'm a new user of R, and I have a question regarding the use of aov and
lm-functions. I'm doing a fractional factorial experiment at our production
site, and I need to familiarize myself with the analysis before I conduct
the experiment. I've been working my way through the examples provided at
http://www.itl.nist.gov/div898/handbook/pri/section4/pri472.htm
2014 Sep 01
1
Correlation Matrix with a Covariate
R Help -
I'm trying to run a correlation matrix with a covariate of "age" and will
at some point will also want to covary other variables concurrently.
I'm using the "psych" package and have tried other methods such as writing
a loop to extract semi-partial correlations, but it does not seem to be
working. How can I accomplish this?
library(psych)
> set.cor(y =
2012 Mar 31
2
lm no calcula un coeficiente
Hola,
quiero hacer una regresión lineal
lm(y ~ x * grupo, data =datos)
y: numérica, x: numérica, grupo: factor con dos niveles (1 y 2)
pero no calcula el coeficiente de x:grupo2 a cuenta de una singularidad
Coefficients: (1 not defined because of singularities)
Estimate Std. Error t value Pr(>|t|)
(Intercept) -1.283e+06 2.276e+04 -56.359 < 2e-16 ***
x
2003 Jul 21
1
problem in pot with subset selection
Hi,
I try to make a selective plot, but it make an error. The subset only work
with length < 4, it is correct?
>
plot((ocorrencia/isca)~frag,subset=especieama==c("grupo1","grupo2","grupo3","AnoplotermesSp1","NeocapritermesOpacus"),pch=c(1,2,3,4,5),xlab="?rea
(ha)",ylab="Propor??o de iscas ocupadas por t?rmitas")