This paper studies robust nonparametric regression, in which an adversarial attacker can modify the values of up to q samples from a training dataset of size N. Our initial solution is an M-estimator based on Huber loss minimization. Compared with simple kernel regression, i.e. the Nadaraya-Watson estimator, this method can significantly weaken the impact of malicious samples on the regression performance. We provide the convergence rate as well as the corresponding minimax lower bound. The result shows that, with proper bandwidth selection, supremum error is minimax optimal. The L2 error is optimal with relatively small q, but is suboptimal with larger q. The reason is that this estimator is vulnerable if there are many attacked samples concentrating in a small region. To address this issue, we propose a correction method by projecting the initial estimate to the space of Lipschitz functions. The final estimate is nearly minimax optimal for arbitrary q, up to a logarithmic factor.
Paper
Similar papers
Peer review
Summary
The paper studies robust nonparametric regression under poisoning attacks. The input is samples from some fixed distribution and the goal is to approximate some unknown function on the input space. It is assumed that the values observed in q of the samples are adversarial. In this setting, classical approaches such as k-NN estimators can fail. This is because a single adversarial sample can affect the prediction on many points. In order to avoid this issue, the paper proposes robust variants of nonparametric regression. The first result is a bound on the convergence rate of an M-estimator based on Huber loss minimization: the initial estimator has minimal optimal ell_infinity risk. When the number of adversarial samples is not too big, the ell_2 risk is optimal. It is also shown that the estimator can suffer when many of the adversarial samples are concentrated within a small region. In order to resolve this issue, the paper proposes a correction step that projects the estimator into the space of Lipschitz functions. Upper bounds on the rate of convergence of the corrected estimator are established. Numerical results show that both estimators outperform standard methods in a simple low-dimensional synthetic scenario.
Strengths
Nonparametric regression is a fundamental problem. The paper studies a natural adversarial setting where a small number of input samples are corrupted adversarially. This research direction can produce more robust useful data-analytic methods. The upper and lower bounds match up to polylogarithmic factors when the number of adversarial samples is not too large.
Weaknesses
The experimental evaluation is very limited. This is unclear whether the proposed methods can lead to significant robustness improvements in real-world scenarios. No runtime analysis is given for the corrected estimator. This is a continuous optimization problem, so it is not clear how easy it is to solve in practice, especially in high or moderately-high dimensions.
Questions
Assumption 1, c: The quantifier for x is missing. In line 253, I don't quite understand the justification for using ell_1 norm in the optimization problem. Since ell_2 imposes a larger cost in the error region, it would seem that this makes it easier to remove spikes. Definition 1 (3): One of the two floors should be a ceiling.
Rating
6: Weak Accept: Technically solid, moderate-to-high impact paper, with no major concerns with respect to evaluation, resources, reproducibility, ethical considerations.
Confidence
3: You are fairly confident in your assessment. It is possible that you did not understand some parts of the submission or that you are unfamiliar with some pieces of related work. Math/other details were not carefully checked.
Soundness
2 fair
Presentation
3 good
Contribution
2 fair
Limitations
The authors have adequately addressed the limitations of their work.
Summary
A robust version of the classic non-parametric problem is studied, when the training data is under an adversarial poisoning attack. The work is primarily theoretical, and makes assumptions on Lipschitzness and boundedness of the underlying function, boundedness of the data density, restrictions on sharpness of the data domain boundary, sub exponential noise. With some restrictions on the kernel, and appropriate parameters, upper bounds are shown on the $\ell_2$ and $\ell_\infty$ loss for the robust regression fit using Huber loss minimization. Lower bounds are obtained under the same assumptions, which match asymptotically with the upper bounds for $\ell_\infty$ and a correction is presented to match the upper bound for $\ell_2$ loss as well. Numerical simulations show usefulness of the algorithms on simulated one and two dimensional data.
Strengths
1. The authors study robust nonparametric regression under poisoning attacks for smooth underlying functions under certain assumptions on the data and domain, under which asymptotically tight rates are obtained for $\ell_\infty$ loss by suitably parameterized Huber loss. 2. A novel correction method is proposed for which tight rates are obtained for $\ell_2$ loss as well. Together with above, the optimal algorithms under the given assumptions are obtained for both losses. 3. Numerical simulations are performed to verify the theoretical results.
Weaknesses
1. The assumption needed for proving the theoretical result, specifically assumption 1(b) is too strong. Effectively, it reduces the possible test distributions to "near-uniform" distributions since there is both an upper and a lower bound on the probability density function. This assumption is too strong to make the results useful in most practical situations. 2. No insights is provided into the proofs of the theoretical results; in particular there are no proof sketches in the main body. This makes it challenging to understand the contributions, given the work is primarily theoretical. 3. Experiments only involve very specific numerical simulations. How does the approach work on real regression datasets? The effect of changing the hyperparameters is not studied and how they are selected is not described. 4. Huber loss is already known as a technique for robust nonparametric regression [1] and is not novel, the paper should mention that only the analysis under the present assumptions is new. The novel correction method is computationally intractable. [1] Maronna, R. A., Martin, R. D., Yohai, V. J., & Salibián-Barrera, M. (2019). Robust statistics: theory and methods (with R). John Wiley & Sons.
Questions
1. What is the intuitive behind needing assumption (c)? Do sharp corners make the model more susceptible to attack by a poisoning adversary? 2. Line 190-191 looks like notational overload for $h$? 3. How do T/M/h settings in the numerical simulations compare to theory?
Rating
5: Borderline accept: Technically solid paper where reasons to accept outweigh reasons to reject, e.g., limited evaluation. Please use sparingly.
Confidence
3: You are fairly confident in your assessment. It is possible that you did not understand some parts of the submission or that you are unfamiliar with some pieces of related work. Math/other details were not carefully checked.
Soundness
3 good
Presentation
2 fair
Contribution
2 fair
Limitations
The authors note several limitations of their work in a dedicated section.
Summary
The authors propose two new nonparametric regression methods that and study their resilience under poisoning attacks. The baseline for comparing the performance of the proposed methods is kernel regression (aka Nadaraya-Watson estimator). The proposed methods are two, since the initial idea that the authors have, turns out to be vulnerable to poisoning attacks that are concentrated in small regions and a large budget is spent there. Hence, the authors arrive at a `corrected' estimator which behaves better and is closer to optimal behavior. After rebuttal: For the largest part I believe that the authors have provided sufficient clarifications to various issues that were raised and moreover are willing to integrate comments and clarifications that came up during the discussion period in the final version of the manuscript. Therefore, I am increasing my score from reject to borderline accept; I am also increasing the soundness from fair to good as well as the presentation from poor to fair.
Strengths
+ New methods for regression that are resilient to poisoning attacks. + Both theoretical and practical results.
Weaknesses
+ Missing a section with preliminaries or background knowledge, where notions and notation that the authors use in the paper is well-defined. For example, in the section for preliminaries you can give information such as: - define loss functions, - discuss kernel regression and separate the presentation from your initial estimator, - discuss notions that you use when it is unclear what these things are; e.g., "bandwidth". - define functions that you use without explanations at the moment; e.g., in (18) we see Clip and med; what are the arguments and what do these functions do? + References not in alphabetical order. + In line 52-54, you indicate that your approach is similar to a combination of an $\ell_1$ and $\ell_2$ loss functions. At that point in the text, I was expecting a comparison with "elastic nets"; you may want to consider adding a small comment, or consider rephrasing and avoid such expectations from the readers. + The order by which some things are presented should probably change. For example, in lines 115-124, I think you need to give a slightly better explanation so that the reader can get better intuition, and moreover, this discussion should come up before you lay out the equations of the method that you propose. + You cannot start sentences with as in lines 167 or 168. You could add a word like "Part" or "Parts" in the beginning and make it read more naturally even if it takes a bit more space. + Larger font size in Figure 1 is expected. + It would be nice to see some commentary on the bounds and argue how they compare against each other and potentially compared to other methods that you cite in the literature. + In general, you have no comparison with parametric methods, either in theory, or in the experiments. Alternatively, provide an explanation as to why you don't have such comparisons. + It is unfortunate and normally it should not be an issue, but the authors need help on writing a paper that has fewer spelling or expression mistakes. The additional problem here is that the authors also have a big appendix, which I suspect is written similarly along the main text. So, fixing the main text is not enough. The paper needs to be proofread by someone in its entirety, including the appendix of the authors. I think that you have a very interesting story to tell, but the paper needs restructuring and better presentation of what you have accomplished.
Questions
Q1. In line 152, you define $D$ but this parameter does not show up anywhere. Q2. In line 243, what is a "linear filter"? Q3. In line 276, what do you mean by "optimal bandwidth"? Q4. In line 308, what do you mean by "number of grids"? Q5. What is $m$ in (18)? Q6. Why don't you have some comparisons with parametric methods?
Rating
5: Borderline accept: Technically solid paper where reasons to accept outweigh reasons to reject, e.g., limited evaluation. Please use sparingly.
Confidence
3: You are fairly confident in your assessment. It is possible that you did not understand some parts of the submission or that you are unfamiliar with some pieces of related work. Math/other details were not carefully checked.
Soundness
3 good
Presentation
2 fair
Contribution
3 good
Limitations
N/A
Summary
This paper studies nonparametric regression, where an adversary can corrupt $q$ samples from the training set. The paper proposes robust estimators for this problem.
Strengths
The paper studied the Huber loss minimization approach under adversarial noise, giving theoretical upper bounds and lower bounds.
Weaknesses
The organization could be improved, e.g. the last two paragraphs of section 3 could conceivably be placed elsewhere. Some aspects of this paper are not really my area, so I cannot provide many helpful comments.
Questions
Equation 2 and 3, what is expectation taken over?
Rating
6: Weak Accept: Technically solid, moderate-to-high impact paper, with no major concerns with respect to evaluation, resources, reproducibility, ethical considerations.
Confidence
2: You are willing to defend your assessment, but it is quite likely that you did not understand the central parts of the submission or that you are unfamiliar with some pieces of related work. Math/other details were not carefully checked.
Soundness
3 good
Presentation
3 good
Contribution
3 good
Limitations
The author addressed limitations.
Thank you for the response
**References not in alphabetical order** I have always interpreted the quoted text as not mixing citation styles (e.g., in some parts one says "Hasselmo et al. (1995), while in other parts says [3]). In the .tex file that is provided as an example by NeurIPS, indeed, [3] is cited first and all the references are in alphabetical order. So, alphabetical order seems to be the right approach; not the order by which you cite the papers in the main text. **Elastic net** All I am saying is that some clarification should appear in the text, which is non-existent in the current form. **Comparison with parametric methods** One could identify the cases/domains where parametric/non-parametric methods are more suitable on real-world datasets and nevertheless compare the methods to see how far off their estimates are. I think it is still interesting to understand how much the results differ in various cases and the extent to which various assumptions hold on real datasets. I appreciate all your responses. Indeed, my opinion for the paper has improved, but I still cannot recommend acceptance.
Thanks for your further response
Thanks for your response. We will definitely take your comment into our consideration and revise our paper accordingly. **References in alphebetical order** We have came across the NeurIPS 2022 accepted paper list: https://papers.nips.cc/paper_files/paper/2022 From the paper list, among the first five NeurIPS accepted papers in the previous year, three of them list references in the order of appearance, instead of in the alphabetical order: https://papers.nips.cc/paper_files/paper/2022/hash/002262941c9edfd472a79298b2ac5e17-Abstract-Conference.html https://papers.nips.cc/paper_files/paper/2022/hash/00295cede6e1600d344b5cd6d9fd4640-Abstract-Conference.html https://papers.nips.cc/paper_files/paper/2022/hash/003a96110b7134d678cb675c6aea6c7d-Abstract-Conference.html Therefore, we believe that the alphabetical order is not required. **Elastic net** Thanks. We understand that it would be better to add some clarification, in order to help readers that are not in this area to understand. Maybe we can change the statement. Instead of claiming "Huber cost is a combination of $\ell_1$ and $\ell_2$ cost", maybe it is better to change the statement to "Huber cost is $\ell_2$ cost exactly with small input, and near to $\ell_1$ cost with large input". This will help to avoid readers from thinking about elastic net. We will make corresponding revision. **Comparison with parametric methods** In lines 45-48, we have mentioned that "in many real scenarios, due to problem complexity or lack of prior knowledge, parametric models are not always available. Therefore, we hope to explore effective schemes to overcome the robustness issue of nonparametric regression." This implies that if we can provide an appropriate parametric assumption, then it would be better to use parametric model. In the cases without such prior knowledge, then nonparametric methods are better. Actually, pros and cons of parametric vs nonparametric statistics are widely discussed in literatures. We agree with the reviewer that it is helpful to add several sentences for clarification, which will be included in our revision. However, the detailed comparison is beyond the scope of our paper.
I must say that I am happy with the responses that the authors have provided to the various reviews, as well as the fact that the authors are willing to do some extra work as it is required by the reviewers for the final version of the paper. For this reason I will increase my score from reject to borderline accept.
Thank you very much!
We are grateful for the improvement of rating, as well as your valuable comments. We are making corresponding revisions. Thanks!
Reply to rebuttal
Thanks for the detailed response. In particular, thank you for answering my questions and highlighting the difference relative to prior work. Regarding the "strong density assumption", as I suspected it is a relatively _strong_ assumption to study. For example, Audibert et al. 2007 also study the "mild density assumption" and also obtain optimal learning rates under it, which is more interesting. It would be good to include a discussion of the related work (thanks for providing it!) on the "strong density assumption" and mention that going beyond this assumption (potentially using the adaptivity techniques you mention) would be interesting. Despite the relatively strong assumption, the techniques could still be interesting to researchers. However, the lack of any proof sketches makes it harder to glean the technical novelties. It would be great to include proof sketches and more experiments (e.g. including real regression data) in your revision. `We are not sure about the meaning of 'notational overload'. ` Regarding the notational overload in lines 190-191, note that $h$ is also the bandwidth parameter in equation (5).
Thank you very much for your response! We have further reply.
Thank you very much for replying to our response! In particular, we sincerely appreciate that you read the related paper of Audibert's. We would like to provide further response. # Strong Density Assumption Our work is the first attempt on adversarial attack on nonparametric regression. **From our experience in this area, compared with other initial works for some problems in nonparametric statistics, our assumption is not strong.** It is common to impose strong density assumption first. Follow-up works will relax these assumptions. In our revision, we will provide more discussions. Going beyond this assumption is indeed interesting, but this requires long and tedious analysis that distract the reader from the main focus. **Example 1.** In Audibert et al. [1], although the authors provided the matching upper bound under mild density assumption, **the purpose is to show that this bound is information-theoretic achievable, instead of providing a practical method to achieve it.** In particular, the matching upper bound is provided in Theorem 4.3 in [1], but this bound is achieved by eq (4.1). Note that (4.1) takes minimum over **a bunch of estimators $\eta\in N_{\epsilon_n}$ in which $N_{\epsilon_n}$ is very large.** Thus this method is not practical. There are several follow-ups including [2], which gives a practical matched upper bound, but **with a much larger set of unlabeled samples** (i.e. semi-supervised learning). Such unlabeled set may not always exist. Even if this unlabeled set exists, the assumption that unlabeled samples follow the same distribution as labeled samples (i.e. labeling is completely random) is not realistic. This work is then followed by [3], which does not require these unlabeled samples. **Example 2.** Another related problem is nonparametric entropy estimation. Assumption 1(c) in [4] is equivalent to the strong density assumption (the proof of equivalence is shown in Appendix G-B in [5]). [4] only analyzes the case with strong density estimation. In follow-up work [5], the assumption on the lower bound of density is removed. **Example 3.** Assumption 2 in [6]. **In our revision, we will definitely include discussions of the related work about other problems in nonparametric statistics, and mention that going beyond this assumption would be interesting.** # Proof sketches and more experiments NeurIPS submission requires 9 pages only for initial submission, thus the space is not enough. When accepted, there will be an additional content page. We will definitely include proof sketches and experiments on real data if accepted. # Notational Overload Thanks for clarification. **Throughout this paper, $h$ always refer to bandwidth**. There is no notational overload here. The estimation error is large if in $B_h(x)$, most of samples are attacked. Lines 190-191 defines $B_h(x)$ instead of $h$. In our revision, we will emphasize it to avoid misleading. Thanks again for the further response! Reference: [1] J.Audibert et al. "Fast learning rates for plug-in classifiers." Annals of Statistics 2007. [2] T. I. Cannings et al. "Local nearest neighbour classification with applications to semi-supervised learning." Annals of Statistics 2020. [3] P.Zhao et al. "Minimax rate optimal adaptive nearest neighbor classification and regression." IEEE Transactions on Information Theory 2021. [4] W.Gao et al. "Demystifying fixed $ k $-nearest neighbor information estimators." IEEE Transactions on Information Theory 2018. [5] P.Zhao et al. "Analysis of KNN information estimators for smooth distributions". IEEE Transactions on Information Theory 2019. [6] H.Jiang "Non-asymptotic uniform rates of consistency for k-nn regression." AAAI 2019.
Thanks for the additional discussion on related work on nonparametric statistics. Also thanks for clarifying the notation in $B_h(\mathbf{x})$, and yes it would be good to use a consistent notation here (e.g. avoid $B(\mathbf{x}, h)$). I have increased my score based on the above discussion.
Thank you very much!
We appreciate much for your valuable comments. In the revision, we have included related discussions in the paper. Moreover, we have corrected the notations such as $B_h(\mathbf{x})$. Thanks for the improvement of rating!
Thanks for the response. I will keep my score for now.
Decision
Reject