Displaying 20 results from an estimated 9000 matches similar to: "L(x,y,t)?"
2024 Aug 09
3
If loop
Can someone help me with the if loop below? In the subroutine, I
initialize all of (joint12,marg1,marg2,cond12,cond21) as FALSE, and call
with only one of them being TRUE:
,...,joint12=FALSE,marg1=FALSE,marg2=FALSE,cond12=FALSE,cond21=FALSE,,,,
joint12 seems to always kick in, even though I call with, e.g., marg1
being TRUE and everything else being FALSE. My attempts with if... else
if were
2024 Aug 09
3
If loop
"Or use <<- assignment I think. (I usually return, but return can only
return one object and I think you want two or more"
You can return any number of objects by putting them in a list and
returning the list.
Use of "<<-" is rarely a good idea in R.
-- Bert
On Fri, Aug 9, 2024 at 1:53?AM CALUM POLWART <polc1410 at gmail.com> wrote:
>
> OK. The fact
2024 Aug 09
1
If loop
The following (using if else) did not help. Seemed like joint12 always
kicked in.
??? me1<-me0<-NULL.
??? if(joint12){
????? {me1<-cbind(me1,v1$p12);? me0<-cbind(me0,v0$p12)}
??? } else if(marg1) {
????? {me1<-cbind(me1,v1$p1);?? me0<-cbind(me0,v0$p1)}
??? } else if(marg2) {
????? {me1<-cbind(me1,v1$p2);?? me0<-cbind(me0,v0$p2)}
??? } else if(cond12){
?????
2024 Aug 09
1
If loop
Is something wrong in the initialisation part that we don't see?
joint12 <- marg1 <-F
marg1 <-T
if (joint12) {
print ("joint 12")
cat (joint12)
}
if (marg1) {
print("marg 1")
cat(marg1)
}
Would probably be my diagnostic approach
On Fri, 9 Aug 2024, 04:45 Steven Yen, <styen at ntu.edu.tw> wrote:
> Can someone help me with the if loop below?
2024 Aug 09
1
If loop
Thanks. Hmm. The loop is doing what it is supposed to do.
> try1<-function(joint12=FALSE,marg1=FALSE,marg2=FALSE,
+??????????????? cond12=FALSE,cond21=FALSE){
+ # ***************************************************
+ # Testing if loop
+ # ***************************************************
+ if(joint12){
+?? {print ("joint12"); cat(joint12,"\n")}
+?? {print
2024 Aug 09
2
If loop
OK. The fact it's in a function is making things clearer.
Are you trying to update the values of an object from within the function,
and have them available outside the function. I don't speak functional
programming articulately enough but basically
v <- 1
funA <- function() {
v <- v+1
}
funA()
cat (v)
# 1
You either return the v from the function so
funB <- function() {
2004 May 14
5
Tagging identical rows of a matrix
I would like to generate a vector having the same length
as the number of rows in a matrix. The vector should contain
an integer indicating the "group" of the row, where identical
matrix rows are in a group, and a unique row has a unique integer.
Thus, for
a <- c(1,2)
b <- c(1,3)
c <- c(1,2)
d <- c(1,2)
e <- c(1,3)
f <- c(2,1)
mat <- rbind(a,b,c,d,e,f)
I would
2005 Oct 10
1
using innov in arima.sim
Hello,
I have used the arima.sim function to generate a lot of time series, but to day I got som results that I didn't quite understand. Generating two time series z0 and z1 as
eps <- rnorm(n, sd=0.03)
z0 <- arima.sim(list(ar=c(0.9)), n=n, innov=eps)
and
z1 <- arima.sim(list(ar=c(0.9)), n=n, sd=0.03),
I would expect z0 and z1 to be qualitatively similar. However, with n=10 the
2003 Jul 14
2
Subsetting a matrix
I'd welcome some comments or advice regarding the situation described
below.
The following illustrates what seems to me to be an inconsistency
in the behaviour of matrix subsetting:
> Z<-matrix(c(1.1,2.1,3.1,1.2,2.2,3.2,1.3,2.3,3.3),nrow=3)
> Z
[,1] [,2] [,3]
[1,] 1.1 1.2 1.3
[2,] 2.1 2.2 2.3
[3,] 3.1 3.2 3.3
> dim(Z)
[1] 3 3
>
2006 Sep 01
1
integration problem with gamma function
Dear R-list members,
I have a problem with translating a mathematica script into R. The whole
script is at the end of the email (with initial values for easy
reproduction) and can be pasted directly into R. The problematic part
(which is included below of course) is
<--- Original Mathematica --->
(* p_svbar *)
UiA = Ni (Dsi - 2Di A + A^2)/2;
UiiA = Nii (Dsii - 2Dii A + A^2)/2;
psvbar =
2011 Sep 07
1
Very slow assignments
I'm creating an object of a S4 class that has two slots: ListExamples, which
is a list, and idx, which is an integer (as the code below).
Then, I read a data.frame file with 10000 (ten thousands) of lines and 10
columns, do some pre-processing and, basically, I store each line as an
element of a list in the slot ListExamples of the S4 object. However, any
kind of assignment operation (<-)
2011 Jul 20
3
select element from each row of the matrix
I have a 5 column matrix like
12 10 8 6 3
10 9 8 7 5
14 NA 4 NA NA NA
15 NA 10 NA 5
...
I want to select the position of the first entry for each row <=5
for example, for the first row, I want to select the last element and return
its position as 5;
for th e third row, I want to select the third element and return its
position as 3;
similarly for the 4th row, I want to select the fifth
2019 Dec 11
3
Weird update_llc_test_checks behavior?
I've run update_llc_test_checks on a set of tests and am seeing some
weird behavior. The CHECK lines appear directly after the function's
first line, even if there are multiple arguments. E.g.:
define <vscale x 4 x i32> @sel_nxv4i32(<vscale x 4 x i1> %p,
; CHECK-LABEL: sel_nxv4i32:
; CHECK: // %bb.0:
; CHECK-NEXT: mov z0.s, p0/m, z1.s
; CHECK-NEXT: ret
2004 Feb 10
2
How to compute the minimal distanct between a point and curve in N-dim space
Dear All,
In the N-dimensional space, give a data point A and a curve f,
how to write the explicit expression for calculating the
minimal distance between A and f?
Or have to use some nonlinear optimization method to calcualte it?
Thanks for your point.
Fred
[[alternative HTML version deleted]]
2010 Nov 29
2
Significance of the difference between two correlation coefficients
Hi,
based on the sample size I want to calculate whether to correlation
coefficients are significantly different or not. I know that as a first step
both coefficients
have to be converted to z values using fisher's z transformation. I have
done this already but I dont know how to further proceed from there.
unlike for correlation coefficients I know that the difference for z values
is
2014 Jul 09
3
error com un archivo
Hola a todos,
Me gustaria pedir vuestra ayuda a encontrar el error que no consigo
encontrar en este archivo. He revisado todo mil veces y probado y no doy
con ello.Adjunto el archivo con Google drive porque es muy grande.
?
monicap_50.csv
<https://docs.google.com/file/d/0B8o2KrPEgG7ATlBMc19lTVk1d3M/edit?usp=drive_web>
?
Este es el script, y lo que no entiendo que pasa es que tengo 592044
2004 Oct 11
1
Bug#275946: Acknowledgement (newline not recognized when logcheck sends emails)
I upgraded to 1.2.28, same results.
Here are the rules I added.
^\w{3} [ :0-9]{11} [._[:alnum:]-]+ perdition\[[0-9]+\]: Connect:
^\w{3} [ :0-9]{11} [._[:alnum:]-]+ pure-ftpd: [^[:space:]]+ \[NOTICE\]
^\w{3} [ :0-9]{11} [._[:alnum:]-]+ pure-ftpd: [^[:space:]]+ \[INFO\]
^\w{3} [ :0-9]{11} [._[:alnum:]-]+ exact\[[0-9]+\]:
^\w{3} [ :0-9]{11} [._[:alnum:]-]+ slapd\[[0-9]+\]:
^\w{3} [ :0-9]{11}
2012 Jun 12
4
How do I connect dots in the RGL package?
Hi all,
I just started using the rgl package in R, and I need to know how I connect
the dots I created. A simple line from one co?rdinate (x0,y0,z0) to another
(x1,y1,z1) will do. I tried using the code rgl.lines, but it doesn't accept
lines that don't stem from O (0,0,0). Please help!
Thanks!
--
View this message in context:
2005 Dec 12
2
More floating point errors
Hi!
It seems that using a different compiler (I assume the
development is done with MSC) uncovers various issues
with the floating point.
Here is the second one I'm seeing - DOMAIN error from
sqrt().
In the open_loop_nbest_pitch(), this line:
g = DIV32(corr[i-start],
10+SHR32(MULT16_16(spx_sqrt(e0),spx_sqrt(energy[i-start])
),6));
The values are:
e0 is 7.2
i is 142,
start is 17
2012 Apr 03
1
object of type 'S4' is not subsettable
hey there!
The object 'cit' contains:
> cit
#####################################################
# Johansen-Procedure Unit Root / Cointegration Test #
#####################################################
The value of the test statistic is: 5.3484 9.0681 10.6433
---------------
I want R to save the value 5.3484 in a new object. I am used to use the
command x=cit[a] where a