functional programming - errors when using fpca package in R -


i using fpca package apply reduced rank model in james et al. (2000), every time try use fpca.score function, returns me error:

error in eigenfuncs.u[gridtime, 1:k] : subscript out of bounds

so pulled out source code:

function (data.m, grids.u, muhat, eigenvals, eigenfuncs, sig2hat, k)  {     temp <- table(data.m[, 1])     n <- length(temp)     m.l <- as.vector(temp)     result <- matrix(0, n, k)     n <- length(grids.u)     evalmat <- diag(eigenvals[1:k])     current <- 0     eigenfuncs.u <- t(eigenfuncs)     data.u <- matrix(as.numeric(as.vector(data.m[, -1])),nrow=nrow(data.m[,-1]), ncol = ncol(data.m[, -1]))     (i in 1:n) {         y <- as.vector(data.u[(current + 1):(current + m.l[i]), 1])         meastime <- data.u[(current + 1):(current + m.l[i]), 2]         gridtime <- ceiling(n * meastime)         muy <- muhat[gridtime]         phiy <- matrix(eigenfuncs.u[gridtime, 1:k], ncol = k)         sigy <- phiy %*% evalmat %*% t(phiy) + sig2hat * diag(m.l[i])         temp.y <- matrix(y - muy)         result[i, ] <- evalmat %*% t(phiy) %*% solve(sigy, temp.y)         current <- current + m.l[i]     }     return(result) }

data.m three-column data matrix first column id, second column measurement, third column time. input of fpca.mle function outputs include other inputs of above fpca.score function: grids, muhat, eigenvals, eigenfuncs, sig2hat, k.

apparently error occurs here 

eigenfuncs.u[gridtime, 1:k]

which comes 

n <- length(grids.u) meastime <- data.u[(current + 1):(current + m.l[i]), 2] gridtime <- ceiling(n * meastime)

is there fda guru knows how fix this? appreciated.

the issue might come mismatch between grids value , values of time in data. using newton algorithm, default value grids seq(0,1,0.002), time data needs re-scaled fall between (0,1).


Comments

Post a Comment

Popular posts from this blog

node.js - Using Node without global install -

i'm working on few node & electron (atom shell) projects , i'm curious how them work without user having install node globally. i've seen few applications including node directory in project contains node.exe binary. how go setting node project use binary? edit: thought should add, i've searched everywhere answer this, unrelated results.
Read more

How to access a php class file from PHPFox framework into javascript code written in simple HTML file? -

i've website developed using phpfox v3.0.7 . now i'm trying implement 'sse(server sent events)' 1 of php files. reference i'm giving below necessary code part of file. file titled process.class.php , present @ location /var/www/module/notification/include/service/process.class.php present on server having ip http://34.124.40.142/ . <?php /** * [phpfox_header] */ header('content-type: text/event-stream'); header('cache-control: no-cache'); defined('phpfox') or exit('no dice!'); /** * * * @copyright [phpfox_copyright] * @author raymond benc * @package phpfox_service * @version $id: service.class.php 67 2009-01-20 11:32:45z raymond_benc $ */ class notification_service_process extends phpfox_service { /** * class constructor */ public function __construct() { $this->_stable = phpfox::gett('notification'); } public function add($...
Read more

java - Null response to php query in android, even though php works properly -

Image
table i'm trying retrieve above data using 2 php files below. want pass email parameter , return bets user = email. work fine when tested advanced rest client, problem not php. however, when try through android 2 java classes below, null pointer exception (logcat below) , json object returned null. i'm not sure error is, think i'm not passing parameter php , that's why it's producing null response. logcat 05-31 15:10:50.227 26074-26074/com.example.albert.betterapp e/androidruntime﹕ fatal exception: main process: com.example.albert.betterapp, pid: 26074 java.lang.runtimeexception: unable start activity componentinfo{com.example.albert.betterapp/com.example.albert.betterapp.displayallbets}: java.lang.nullpointerexception @ android.app.activitythread.performlaunchactivity(activitythread.java:2292) @ android.app.activitythread.handlelaunchactivity(activitythread.java:2350) @ android.app.activitythread.acce...
Read more