How to apply a function to each entry of a matrix in R -


i have store log likelihood of densities in matrix, follows:

matrix.loglik              [,1]      [,2]      [,3]     [,4] [,5] [1,]     0.00000    0.0000    0.0000   0.0000    0 [2,]   -34.41018    0.0000    0.0000   0.0000    0 [3,] -2275.14867 -765.8642    0.0000   0.0000    0 [4,]    64.96982  264.7709 -256.1461   0.0000    0 [5,]   358.17822  260.1582  427.3490 363.2247    0

i apply aic function. 

aic.log <- function(x,y=2){   -2*x+2*y }

where x log likelihood value (the entries of matrix.loglik) , y number of parameters. then, store result in lower triangular matrix similar matrix.loglik.

thank answer. however, matrix must stay lower triangular matrix. zeros entries in matrix must still zeros.

i try both answer , got this:

              [,1]      [,2]      [,3]      [,4] [,5]     [1,]    4.00000    4.0000    4.0000    4.0000    4     [2,]   72.82036    4.0000    4.0000    4.0000    4     [3,] 4554.29734 1535.7284    4.0000    4.0000    4     [4,] -125.93964 -525.5418  516.2922    4.0000    4     [5,] -712.35644 -516.3164 -850.6980 -722.4494    4

this not should get.

number 4 must zero. idea please?

you can multiply `lower.tri` desired result:  matrix.loglik <- c(       0.00000,    0.0000,    0.0000,   0.0000,    0,      34.41018,    0.0000,    0.0000,   0.0000,    0,   -2275.14867, -765.8642,    0.0000,   0.0000,    0,      64.96982,  264.7709, -256.1461,   0.0000,    0,     358.17822,  260.1582,  427.3490, 363.2247,    0)  matrix.loglik <- matrix(matrix.loglik, nrow = 5, byrow = true)   aic.log <- function(x, y = 2){   -2 * x + 2 * y }   aic.log(matrix.loglik) * lower.tri(matrix.loglik)  #           [,1]      [,2]      [,3]      [,4] [,5] #[1,]    0.00000    0.0000    0.0000    0.0000    0 #[2,]  -64.82036    0.0000    0.0000    0.0000    0 #[3,] 4554.29734 1535.7284    0.0000    0.0000    0 #[4,] -125.93964 -525.5418  516.2922    0.0000    0 #[5,] -712.35644 -516.3164 -850.6980 -722.4494    0

another way modify aic.log function:

aic.log <- function(x, y = 2){   m <- -2 * x + 2 * y   lower.tri(m) * m }  aic.log(matrix.loglik)

Comments

Post a Comment

Popular posts from this blog

python Tkinter Capturing keyboard events save as one single string -

is possible store keyboard events 1 string? code below stores 1 char , prints it. card reader or bar code reader, contains collection of character/ string not 1 character @ time. goal save char pressed text variable. from tkinter import * root = tk() def key(event): text= event.char text+= event.char print ("pressed", text) def callback(event): frame.focus_set() print ("clicked at", event.x, event.y) frame = frame(root, width=100, height=100) frame.bind("<key>", key) frame.bind("<button-1>", callback) frame.pack() root.mainloop() currently, creating text variable , everytime key() function gets called, , text stores last character typed. you can define text module level variable , , use module level text inside key function - text = '' def key(event): global text text+= event.char print("pressed", text)
Read more

android - InAppBilling registering BroadcastReceiver in AndroidManifest -

i'm implementing inapp billing in android app. works well, however, i'm trying decouple broadcast receiver activity manifest. suggestion in android's trivialdrive sample: // important: dynamically register broadcast messages updated purchases. // register receiver here instead of <receiver> in manifest // because call getpurchases() @ startup, therefore can ignore // broadcasts sent while app isn't running. // note: registering listener in activity bad idea, done here // because sample. regardless, receiver must registered after // iabhelper setup, before first call getpurchases(). currently there's class extends broadcastreceiver : public class iabbroadcastreceiver extends broadcastreceiver { /** * intent action receiver should filter for. */ public static final string action = "com.android.vending.billing.purchases_updated"; private final iabbroadcastlistener mlistener; public iabbroadcastreceiver(iabbroadcastlistener listener) { m
Read more

javascript - Z-index in d3.js -

i have map , when hover effect on region of map below yellow line , circle, last ones disappear. how can solve problem? http://jsfiddle.net/ot21n9qx/ var arc = g.append("path") .style("fill", "none") .style("stroke", "yellow") .style("stroke-width", 2) .attr("d", "m" + pathorigin[0] + "," + pathorigin[1] + " q" + svgpoint[0] + "," + pathorigin[1] + " " + svgpoint[0] + "," + svgpoint[1]); var circlesize = d3.scale.linear() .domain([0, 0.5, 1]) .range([4, 10, 4]); svg not have z order draws objects in order created. need create path , circle after create highlight regions. here forked fiddle: http://jsfiddle.net/o00fzgaf/1/
Read more