search for: j's

Hello everyone, I want to calculate type I error rate for modified F statistic for one way robust anova. I need to find the j group trimmed mean and winsorized sum of squared deviations. Here I attached my code for j=2 to make it simple. Originally I have j=4. Hope you can help. I need to run it for 1000 times My problem is: i) the value of F-test obtain from my simulation below is in negative value..There might be som...

...for some reason, I can't seem to get it to work. I've included copies of my code below. The first set of code uses a loop (and it works fine), and the second set of code attempts to use mapply but I get a "subscript out of bounds" error. Any guidance would be greatly appreciated. Xj, Yj, and Wj are also lists, and s2, TAU, and GAMMA are scalars.   Thank You.   # THIS WORKS USING THE LOOP for (j in 1:J) { V.tilde.j <- solve((1/s2)*t(Xj[[j]])%*%Xj[[j]] + solve(TAU)) # Not singular case: if(round(det(t(Xj[[j]])%*%Xj[[j]]),8)!=0) {  Beta.hat.j <- solve(t(Xj[[j]])%*%Xj[[j]])%...

...y I want to create a probability matrix to be used in a trinomial tree going forward. This is the reason why I thought to build the matrix around 0 would be much more efficient. I need to loop through because the probabilities will depend on my node and is not always the same per row (e.g. if N> jmax, jmax being defined in another function) I found a workaround. Please see below. Thereafter I want to optimize this function. Hopefully it works. Many thanks, Eric HWMProb <- function(N) { ProbUP<- function( a, j, dt) 1 / 6 + ((j ^ 2 * Mfactor(a, dt) ^ 2 + j * Mfactor(a, dt)) / 2)...

I want to populate the matrix prb through the function HWMProb <- function (a,j,dt) that encapsulates different functions (please see code below), using j= 0:2 for each j. It only populates prb if I specify each function independently in the global environment and then run the loop with the iF statement, as per below. for (j in 0:2) { if (j==0) { prb["0",&quot...

<div> </div><div> </div><div>Hello,</div><div> </div><div>To me this is an issue in llvm loop vectorizer (if N is large enough to prevent complete unrolling of j-loop).</div><div> </div><div>Woud you mind to share stencil.ll than I would say more definitely what the issue is.</div><div> </div><div>Regards,</div><div>Serge.</div><div> </div><div> </div><div>22.10....

...eld so I cannot comment your approach. I wanted only to point out that building matrix your way is like scratching your left ear with right hand, especially in R. What if you want increase size of your matrix? E.g. you use function ProbUP once for row "0" and than for rows different from jmax (if I correctly understand your code). Use of any function depends on two parameters (row and column, let them call u and v) > ProbUP<- function( a, j, dt) 1 / 6 + ((j ^ 2 * Mfactor(a, dt) ^ 2 + j * Mfactor(a, > dt)) / 2) ######### Probability X going up > ProbMID<-...

...(30) - Linux box with Samba 2.0.6 as print / fax server Goal: ----- - To send FAXes from Office applications How we did it: -------------- in smb.conf: [fax] comment = FAX Service print ok = yes path = /ub/pkg/fax/spool/spool public = yes writable = no printer driver = HP LaserJet 4/4M PS print command = /bin/sh /ub/pkg/fax/spool/able/Wrap S %s %p %j %U lpq command = /bin/sh /ub/pkg/fax/spool/able/Wrap L %s %p %j %U lprm command = /bin/sh /ub/pkg/fax/spool/able/Wrap R %s %p %j %U so we call 'Wrap' with several parameters, %s to have the print file, %U to...

...t; stencil.ll is attached here. > > On Mon, Oct 23, 2017 at 11:37 AM, Serge Preis <spreis at yandex-team.ru> wrote: >> >> >> >> Hello, >> >> To me this is an issue in llvm loop vectorizer (if N is large enough to >> prevent complete unrolling of j-loop). >> >> Woud you mind to share stencil.ll than I would say more definitely what >> the issue is. >> >> Regards, >> Serge. >> >> >> 22.10.2017, 03:21, "hameeza ahmed" <hahmed2305 at gmail.com>: >> >> Hello, >&g...

...id(u=-2:2, v=-1:1) for each row you can than use one function with parameters u, v and a and dt. After you calculate vector of results, you can easily transform it to matrix by setting dim argument to it. Cheers Petr > -----Original Message----- > From: R-help [mailto:r-help-bounces at r-project.org] On Behalf Of > Eric.Pueyo at avivainvestors.com > Sent: Wednesday, November 1, 2017 1:31 PM > To: r-help at r-project.org > Subject: [R] repeat a function > > I want to populate the matrix prb through the function HWMProb <- function > (a,j,dt) that encapsulates diff...

...The existing code does the work with loops like the following: for (i = 0; i < len; i++) REAL(r)[i] = REAL(a)[(i / ncol) + (i % ncol) * nrow]; It seems a bit optimistic to expect a compiler to produce good code from this. I've re-written these loops as follows: for (i = 0, j = 0; i<len; i += 1, j += nrow) { if (j>=len) j -= (len-1); REAL(r)[i] = REAL(a)[j]; } The resulting improvement is sometimes dramatic. Here's a test program: M <- matrix(seq(0,1,12000),200,60) print(system.time({for (i in 1:1...

...forward -N ufw-after-logging-input -N ufw-after-logging-output -N ufw-after-output -N ufw-before-forward -N ufw-before-input -N ufw-before-logging-forward -N ufw-before-logging-input -N ufw-before-logging-output -N ufw-before-output -N ufw-logging-allow -N ufw-logging-deny -N ufw-not-local -N ufw-reject-forward -N ufw-reject-input -N ufw-reject-output -N ufw-skip-to-policy-forward -N ufw-skip-to-policy-input -N ufw-skip-to-policy-output -N ufw-track-forward -N ufw-track-input -N ufw-track-output -N ufw-user-forward -N ufw-user-input -N ufw-user-limit -N ufw-user-limit-accept -N ufw-user-logging...

...355 2329 890 372 2329 266 719 2329 324 761 2329 189 890 2329 355 913 2329 Below is what I have written: ordering <- vector("list", length(unique(mydata$group))) for (j in 1:length(unique(mydata$group))){ group.j <- mydata[mydata$group == unique(mydata$group)[j], ] ordering.j <- c() ordering.j[1] <- ifelse(group.j[1, ]$prev_event == 0, group.j[1, ]$event, group.j[1, ]$prev_event) for (i in 2:nrow(group.j)){ ordering.j[i] <- group.j[group.j$prev_event =...

Ok, before I definetly give up, and throw the laptop out of the window, or fill my data.frame manually, I'll ask for some help. I have a data.frame named MyTable with 3 columns, that looks like this: V1 V2 V3 red-j appearanceblood-n 105.032 red-j appearanceground-n 93.749 red-j appearancesea-n 102.167 red-j appearancesky-n 10.898 orange-j appearanceobject-n 109.354 orange-j as_adj_aspainting-n 93.248 orange-j colorbanknote-n 159.167 orange-j colorcar-n 117.985 yellow-j appearanceskin-n 109.5...

...39;'m looking for a way to route the packets to the right interface. Google gave my some solutions but no ones are working. Here are my iptable # Generated by iptables-save v1.2.11 on Fri Nov 25 12:21:59 2005 *filter :INPUT DROP [2:184] :FORWARD DROP [0:0] :OUTPUT DROP [3:188] -A INPUT -i lo -j ACCEPT -A INPUT -s 127.0.0.0/255.0.0.0 -i ! lo -j LOG -A INPUT -s 127.0.0.0/255.0.0.0 -i ! lo -j DROP -A INPUT -d 255.255.255.255 -i br0 -j ACCEPT -A INPUT -s 192.168.1.0/255.255.255.0 -i br0 -j ACCEPT -A INPUT -d 224.0.0.0/240.0.0.0 -i br0 -p ! tcp -j ACCEPT -A INPUT -s 192.168.1.0/255.255.255.0 -...

...using the iptables script provided by iptablesrocks.org. It's been quite useful, but I ran into a problem when I tried to connect samba. Without any iptables rules, I have no problem when connecting host os(ubuntu samba server) from guest os Windows XP. I referenced this article, http://troy.jdmz.net/samba/fw/, so I put the following sources in the middle of the source. -A INPUT -p udp -m udp --dport 137 -j ACCEPT -A INPUT -p udp -m udp --dport 138 -j ACCEPT -A INPUT -p tcp -m tcp --dport 139 -j ACCEPT -A INPUT -p tcp -m tcp --dport 445 -j ACCEPT Finally, it doesn't work. I fe...

...########################################## ############################################################################### # aplicando as regras contra ataques de negacao de servico # contra ping of death # iptables -t nat -A PREROUTING -i eth1 -p icmp --icmp-type echo-request -m limit --limit 1/s -j ACCEPT # contra ataques syn-flood # iptables -t nat -A PREROUTING -i eth1 -p tcp -m limit --limit 1/s -j ACCEPT # contra port scanners avancados (nmap) # iptables -t nat -A PREROUTING -i eth1 -p tcp --tcp-flags SYN,ACK FIN,RST -m limit --limit 1/s -j ACCEPT #########################################...

Does it happen due to loop carried dependence? if yes what is the solution to vectorize such codes? please reply. i m waiting. On Jul 1, 2017 12:30 PM, "hameeza ahmed" <hahmed2305 at gmail.com> wrote: > I even tried polly but still my llvm IR does not contain vector > instructions. i used the following command; > > clang -S -emit-llvm stencil.c -march=knl -O3 -mllvm -polly -mllvm > -polly-vector...

...on of R. The difference between them is that in the 'broken' version in line 42 I check if the items in the two-dimensional array are bigger than a certain value, and in the working one I do it in a separate loop. A 2.1.1 build for Solaris returns the following error Error in if ((X - (j - 1) * dS) > f[i, j]) f[i, j] = (X - (j - 1) * dS) : missing value where TRUE/FALSE needed On Windows both the stable 2.2.1 and 2.3.0rc gui versions will silently produce incorrect data. Nikolai.

...-f /var/lock/subsys/ipchains ] && /etc/init.d/ipchains stop) >/dev/null 2>&1 || true (rmmod ipchains) >/dev/null 2>&1 || true /sbin/iptables -F /sbin/iptables -X /sbin/iptables -Z /sbin/iptables -P INPUT DROP /sbin/iptables -A INPUT -m state --state ESTABLISHED,RELATED -j ACCEPT /sbin/iptables -A INPUT -p tcp ! --syn -j REJECT --reject-with tcp-reset /sbin/iptables -A INPUT -m state --state INVALID -j DROP /sbin/iptables -P OUTPUT DROP /sbin/iptables -A OUTPUT -m state --state ESTABLISHED,RELATED -j ACCEPT /sbin/iptables -A OUTPUT -p tcp ! --syn -j REJECT --reject-w...

...st with 12 matrices (representing one month). In each matrix there are as many rows as there are hours in the different months and 2 columns, since two meteorological parameters are measured. What I want to do is to calculate the hourly mean values for each month over the whole 13 year period for (j in 1:12){ (list[[1]][[j]] [1,1] + list[[2]][[j]] [1,1] + ...... list[[13]][[j]][1,1] )/13 } How can I do that nicely? Original code looks like this (zeilen[j] is the number of hours for j=1 = January, 2 = february and so on): mean_list<-list() for (j in 1:12){ mean_list[[j]]<-matri...