DDoS Attack Threat

May 26th, 2016

I recently received an email who threatens to attack a server I am running - unless I pay them 0.2 Bitcoin (about $93 US dollars). Anybody know how to handle this?

Below is the email:

—–

From: PhantomSquad@openmailbox.org
subject line: DDoS Attack Imminent - Important information

FORWARD THIS MAIL TO WHOEVER IS IMPORTANT IN YOUR COMPANY AND CAN MAKE DECISION!

We are Phantom Squad
http://lmgtfy.com/?q=Phantom+Squad

Your network will be DDoS-ed starting June 1st 2016 if you don’t pay protection fee - 0.2 Bitcoin @ 14hxG5TBZz6jCZz4H7mzfzXfgXjgaF3MqA.

If you don’t pay by June 1st 2016, attack will start, yours service going down permanently price to stop will increase to 20 BTC and will go up 10 BTC for every day of attack.

This is not a joke.

Our attacks are extremely powerful - sometimes over 1 Tbps per second. And we pass CloudFlare and others remote protections! So, no cheap protection will help.

Prevent it all with just 0.2 BTC @ 14hxG5TBZz6jCZz4H7mzfzXfgXjgaF3MqA

Do not reply, we will not read. Pay and we will know its you. AND YOU WILL NEVER AGAIN HEAR FROM US!

Bitcoin is anonymous, nobody will ever know you cooperated.

Author: Xu Cui Categories: web Tags:

Who cited my paper?

May 9th, 2016

Back in 2010 we published a paper titled “Functional Near Infrared Spectroscopy (NIRS) signal improvement based on negative correlation between oxygenated and deoxygenated hemoglobin dynamics“. It is about a really simple method which surprisingly works well in reducing motion artifact (noise) in NIRS data. As of today (2016-05-09), the paper has been cited by more than 130 times.

Who cited this paper? Well, we created this map to show where the authors are:

If you clicks the red dots you will find the detailed information about the paper, such as the author name and journal.

Then who are the top researchers who cited this paper? We count the number of citation for each author, and rank them by the number of papers. Here are the top authors who cited our 2010 paper:

Author Number of Papers
Fallgatter, AJ 9
Ehlis, AC 8
Hong, KS 8
Dresler, T 8
Herrmann, MJ 7
Scheutz, M 6
Strait, M 5
Boas, DA 5
Sato, H 4
Scholkmann, F 4
Herff, C 4
Seghouane, AK 4
Ge, SS 4
Shah, A 4
Katura, T 4
Molavi, B 4
Wolf, M 4
Author: Xu Cui Categories: brain, nirs, programming Tags:

NIRS manufactures locations

May 4th, 2016

Where are the major NIRS device manufactures? You can see from the map below.

There are NIRS manufactures in US, Japan, Europe and Korea.

If you are a NIRS manufacture and would like to add to this list, please let me know (leave a comment in this post).
If you are interested in creating a map like this, please let me know (leave a comment in this post).

Related: NIRS manufactures products http://www.alivelearn.net/?p=1335

Author: Xu Cui Categories: brain, nirs Tags:

Stork is my best research assistant (2): Grant alert

April 15th, 2016
  1. Does my boss have money?
  2. I am looking for a postdoc position; does my future boss have enough funding to support me?
  3. How much money was awarded to my field (e.g. NIRS)? And who got the money? What are they going to do with the money?

Have you ever wondered these questions? In the early years as a graduate students, I rarely asked “money” questions. It does not sound what a “true” scientist should care.  I was even puzzled when I realized my boss spent more than half of his time writing grant applications - shouldn’t he spend most of his time doing experiments and write papers?

As a postdoc I found myself spend a lot of time writing grant applications; and realized my career is critically depending on the success of securing enough funding. I also see a few colleagues had to leave academia due to lack of funding.  It would be nice if there is a tool which can notify me of the funding situation in a timely manner.

Stork is such a tool.

I entered some keywords into Stork, including “pearl chiu” (my former colleague) and “nirs brain” (my research field). Below is a letter I got from Stork:

Stork notifies me of awarded grants

Stork notifies me of awarded grants

With the information Stork provides, I know who in our field got grants and what they proposed. In fact the 3rd one is my colleague Manish who is interested in using NIRS in resting-state brain network study. I also got to know Pearl got a big gran, so I sent her a congratulation note.

Compared to journal papers alert, grants alert helps me to know the trend of my field much earlier. This is because publications are usually a few years delayed from grants.

If you also want to be the first one to know new grants in your field, why don’t you give Stork a try? I’m sure you’ll be delighted!

Author: Xu Cui Categories: brain, life, programming, web, writing Tags:

Projection of a NIRS channel on a brain surface

April 11th, 2016

When analyzing the data in our concurrent NIRS-fMRI study, we are particularly interested in how the NIRS signals were correlated to the fMRI signal. To answer the question we need to create an ROI (region of interest) in brain which is directly underneath the NIRS channel (which is on the skull). So a projection of the NIRS channel on the brain surface is necessary.

It might be easy if the brain were a perfect sphere, or at least doable if it is smooth. But brain surface is anything but smooth. What we did is:

  1. Create a brain surface mask. This is easy using BET (http://fsl.fmrib.ox.ac.uk/fsl/fslwiki/BET)
  2. Loop over all voxels on the surface, calculate the distance between the NIRS channel (point A) and the voxel
  3. Find the voxel which is closest to the NIRS channel (It is point B in the figure below)

Projection of NIRS channel on brain surface

Projection of NIRS channel on brain surface

Author: Xu Cui Categories: brain, matlab, nirs, programming Tags:

Is measuring both oxy and deoxy-Hb useful?

March 28th, 2016

One of the advantages of NIRS is that it measure two values at a time - both oxy- and deoxy-Hb concentration, while fMRI measure only one value, called BOLD. (See a comparison between NIRS and fMRI.) However, what is the use of measuring two values if they are highly correlated (actually anti-correlated)?

This is a valid question. In our 2010 paper titled “Functional Near Infrared Spectroscopy (NIRS) signal improvement based on negative correlation between oxygenated and deoxygenated hemoglobin dynamics“, we have shown that (1) theoretically (using Balloon model) oxy and deoxy Hb should be highly correlated, close to -1; (2) with real data, they are highly correlated. See the figure below with real data (I was the subject).

oxy and deoxy-Hb are highly (anti)correlated

oxy and deoxy-Hb are highly (anti)correlated

Now go back to our question: if they are highly correlated, maybe measuring both does not bring additional advantage. If we look at NIRS publications, we see most publications only use either oxy-Hb or deoxy-Hb in their analysis, or at least use them separately. It is not very clear why we should use both.

However, it turns out it is an advantage to measure both. In a study we have tried to “predict” if the subject is doing finger tapping or not by looking at the NIRS signal only (see the paper). This is useful for brain-computer-interface purpose. In the figure below, you will see that if we use only oxy-Hb (blue) or deoxy-Hb (green), the prediction accuracy is lower than if we combine them (red).

Combined oxy- and deoxy-Hb predict better

Combined oxy- and deoxy-Hb predict better

This means while oxy- and deoxy-Hb are highly correlated, we still get more if we measure both. I think part of the reason is that there are always noise in real data. With two measurement the signals are more robust again the corruption of noise.

Author: Xu Cui Categories: brain, nirs Tags:

T value of a single subject

March 14th, 2016

I usually report a group-level T-test image in the final publications or presentations. But in the early stage of a project when only one or a few subjects have been scanned, I often need to report an activation map for an individual subject after general linear model analysis (GLM). I could show a beta image or activation image, but a T-image is sometimes desired. The question is how do I calculate the T-value of a contrast for a single subject.

I will use the following example to show how it is calculated. Assume we have two conditions, beep and flash, and the brain signal is saved in y. Based on the timing of beep and flash, we have an independent variable called x with two columns. The 1st column is for beep, the 2nd for flash.

x=rand(10,2);
y=x(:,1)*2 + x(:,2)*3 + 0.2 + rand(size(x(:,1)))/5;
[a,b,c] = glmfit(x,y);

Now you can easily find the beta value, the T value, and the p-value of each individual condition:

beta = a(2:end)'
T = c.t(2:end)'
p = c.p(2:end)'

What about the contrast between beep and flash? Well, it’s a bit more complicated. First, we need to find the covariance matrix of the betas.

covb = c.covb(2:end, 2:end);

Then we define our contrast vector, in this case it is simply contrast = [1 -1]. The variance of contrast (beep-flash) is

V = contrast  * covb * contrast';

The T-value of the contrast is

T = (beta * contrast') / sqrt(V) * sqrt(length(x));

Note: you might be tempted to use the two-sample T formula to calculate the T-value (i.e. finding the variance of beta1, variance of beta2, and then calculate). It is not the correct way. We have to consider the case where the two conditions are correlated. Using the covariance matrix above is the right way.

Author: Xu Cui Categories: brain, matlab, nirs Tags:

FDR correction matlab script

March 14th, 2016

The following script is to determine is a channel (or voxel) survives the FDR corrected threshold. Assume you have 48 channels and you already calculated the (uncorrected) p-value of each channel. You put it into an array called p. Now you want to know which channel will survive if you do a FDR correction at q=0.05. You may simply use

fdr0(p, 0.05)

The output is an array of 48 numbers, with 1s in the channel where it survives the FDR threshold.

function y = fdr0(p, q)
% y = fdr0(p, q)
%
% to calculate whether a pvalue survive FDR corrected q
%
% p: an array of p values. (e.g. p values for each channel)
% q: desired FDR threshold (typically 0.05)
% y: an array of the same size with p with only two possible values. 0
% means this position (channel) does not survive the threshold, 1 mean it
% survives
%
% Ref:
% Genovese et al. (2002). Thresholding statistical maps in functional
% neuroimaging using the false discovery rate. Neuroimage, 15:722-786.
%
% Example:
%   y = fdr0(rand(10,1),0.5);
%
% Xu Cui
% 2016/3/14
%

pvalue = p;
y = 0 * p;

[sortedpvalue, sortedposition] = sort(pvalue);
v = length(sortedposition);
for ii=1:v
    if q*ii/v >= sortedpvalue(ii)
        y(sortedposition(ii)) = 1;
    end
end

return;
Author: Xu Cui Categories: brain, nirs Tags:

Finger tapping task MatLab script

February 29th, 2016

finger tapping

finger tapping

Finger tapping is probably the most used task in brain imaging studies (fMRI and NIRS). The task is simple and elicits robust brain signal in the motor cortex. We always use it to test new devices or develop a new method.

If you want to control your NIRS device (say Hitachi ETG 4000) to start/end, please refer http://www.alivelearn.net/?p=1260

If you want to see how blood flow increases in brain motor cortex during finger tapping, check out
http://www.alivelearn.net/?p=1686

To download, please fill the form below:






Author: Xu Cui Categories: brain, matlab, nirs, programming Tags:

Hyperscanning data

February 11th, 2016

Below is sample hyperscanning data from one subject (myself). For detailed information please refer to http://www.ncbi.nlm.nih.gov/pmc/articles/PMC3254802






Author: Xu Cui Categories: brain, nirs Tags: